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ABSTRACT 


Wall  pressure  fluctuations  within  an  equilibrium  rectangular  channel  flow  geometry  were 
investigated.  The  overall  objective  of  this  investigation  was  to  obtain  experimentally  for  the 
first  time,  the  statistical  properties  of  the  wall  pressure  fluctuations  for  a  generic  turbulent  flow 
within  a  rectangular  channel  and  to  directly  compare  the  results  with  those  obtained  from  a  full 
solution  to  the  Navier-Stokes  equations  via  Direct  Numerical  Simulation  (DNS).  For  both  the 
physical  and  computational  experiments,  serious  limitations  on  the  quality  of  the  results  were 
identified  related  to  facility  generated  noise,  transducer  resolution,  Reynolds  number  effects 
and  constraints  imposed  by  the  computational  technique. 

An  extensive  data  base  was  obtained  from  the  physical  experiment.  Profiles  for  the  mean 
and  fluctuation  components  of  the  streamwise  velocity  were  obtained  at  a  number  of  spanwise 
locations  and  served  to  establish  the  existence  of  equilibrium  flow  conditions  over  a  range  of 
Reynolds  numbers.  An  array  of  flush-mounted  pressure  transducers  was  utilized  to  measure 
the  two  point  spectral  characteristics  in  the  streamwise  direction.  A  unique  signal  process¬ 
ing  scheme  was  developed  and  used  to  cancel  facility  generated  noise  which  resulted  in  new 
information  being  recorded  at  low  frequencies  for  the  Reynolds  numbers  investigated. 

Two  issues  of  fundamental  importance  were  examined.  First,  the  results  from  both  ex¬ 
periments  were  compared  directly.  The  computational  results  were  found  to  be  somewhat 
underdeveloped  for  the  higher  order  statistics  for  the  particular  data  investigated.  Resolution 
problems  were  shown  to  exist  at  both  ends  of  the  frequency  spectrum.  The  computational 
grid  spacing  was  insufficient  to  resolve  the  small  scale  structure  at  high  frequencies  and  the 
limited  size  of  the  computational  box  resulted  in  contamination  at  the  low  end  of  the  frequency 
spectrum  by  sensing  the  periodic  boundary  conditions  as  “noise”  peaks.  Secondly,  the  spectral 
results  for  rectangular  channel  flow  were  compared  directly  to  results  obtained  from  typical 
outer  flow  regimes  as  represented  by  flat  plate  boundary  layers.  Two  fundamental  results 
were  established.  Although  similarity  scaling  was  shown  to  exist,  the  turbulence  structure 
was  shown  to  decay  approximately  twice  as  fast  as  obtained  in  external  flows.  Various  scaling 
laws  were  examined  which  were  found  to  collapse  respective  spectral  content  separately,  but  no 
“universal”  law  was  found  to  exist  which  brings  both  internal  and  external  flow  regimes  into 
coincidence. 
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This  investigation  has  provided  new  information  not  previously  recorded  in  the  literature. 
An  extensive  high  quality  data  base  was  obtained  for  generic  channel  flow  at  low  Reynolds 
numbers  which  may  be  used  to  assess  the  quality  of  typical  Direct  Numerical  Simulation  (DNS) 
results.  A  careful  examination  of  the  results  will  prove  useful  in  characterizing  the  proper 
governing  parameters  necessary  for  future  DNS  initiatives. 
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PHYSICAL  AND  COMPUTATIONAL  INVESTIGATION  OF  THE 
WALL  PRESSURE  FLUCTUATIONS  IN  A  CHANNEL  FLOW 

CHAPTER  1  INTRODUCTION 

1.1  Background 


Wall  pressure  has  received  a  considerable  amount  of  attention  by  experimentalists  in  the 
last  few  years.  Most  of  this  interest  has  attempted  to  characterize  and/or  otherwise  explain 
fluid/structure  interaction  in  relation  to  the  coupling  phenomena  as  represented  by  the  gen¬ 
eration  of  vibration  and  sound.  In  a  practical  sense,  the  origins  of  these  phenomena  may  be 
expected  whenever  there  is  relative  motion  between  adjacent  bodies  of  fluid,  such  as  a  free  shear 
flow,  or  between  a  fluid  flowing  within  or  external  to  a  body.  If  the  flow  and  subsequent  body 
motions  (such  as  vibrating  panels)  become  coupled,  a  particularly  deleterious  condition  termed 
“flow'-induced  noise”  may  exist.  The  resultant  motion  of  the  body  may  approach  resonance  in 
which  case  the  motion  itself  becomes  quite  large.  In  these  circumstances,  fatigue  of  the  body’s 
surface  may  occur,  or  at  the  least,  high  intensity  sound  waves  radiate  to  the  farfield.  Thus, 
any  study  related  to  these  physical  features  of  fluid/structure  coupling  can  be  characterized  as 
the  interactions  of  the  following  basic  physical  phenomena: 

o  Unsteady  fluid  mechanics  (turbulence) 

o  Structural  response  (impedance,  energy  distribution,  damping) 
o  Radiative  acoustic  field  related  to  fluid/structure  interaction 

A  major  motivation  of  studying  such  flow  interactions  lies  in  a  desire  to  improve  the  basic 
understanding  of  the  structure  of  turbulent  flows  and  their  relation  to  flow-induced  noise,  and 
help  to  resolve  problems  of  practical  engineering  significance.  An  example  of  an  engineering 
application  may  be  found  in  the  area  of  hydroacoustics  where  the  sound  radiating  to  the  farfield 
is  attributed  to  pressure  fluctuations  at  the  wall  coupling  with  the  natural  resonances  of  the 
wall  structure.  A  related  example  of  extreme  importance  is  the  “flow  noise”  associated  with 
the  flow  of  water  over  a  sonar  platform  or  hydrophones  mounted  on  the  surfaces  of  6hips 
or  submarines.  Flow  noise  takes  on  special  significance  in  towed  array  topology,  where  the 
turbulent  flow  about  the  array  so  contaminates  the  pressure  signature  being  recorded  that  it 
becomes  difficult  to  extract  worthwhile  information. 

Since  the  early  work  of  Kraichnan  (1956a, b),  disagreement  has  persisted  on  the  identifica¬ 
tion  of  the  turbulence  source  terms  that  uniquely  contribute  to  the  wall  pressure  field  and  the 
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length  scales  and  locations  within  the  boundary  layer  in  which  this  energy  transfer  takes  place. 
It  was  originally  thought  that  the  spatial  structure  of  the  pressure  field  was  not  as  fine-grained 
as  that  of  the  velocity  field  and  received  most  of  its  energy  from  the  low  frequency  velocity 
spectrum  (Corcos,  1963).  If  one  examines  the  Poisson  equation  for  the  wall  pressure  fluctua¬ 
tions,  it  can  be  shown  that  the  pressure  generated  at  the  wall  is  coupled  to  velocity  fluctuations 
throughout  the  boundary  layer.  Furthermore,  integrating  this  equation  over  the  unbounded 
half-space  above  the  boundary  wall  shows  that  effective  contributions  decrease  with  distance 
away  from  the  wall.  Typically,  80  percent  of  the  energy  contributed  to  the  wall  pressure  ap¬ 
pears  to  occur  within  the  near  wall  buffer  region  (F+  <  50)  at  length  scales  much  smaller  than 
the  boundary  layer  thickness.  Recently,  it  has  been  argued  that  the  Burst/Sweep  events  in 
the  near-w’all  region  are  a  major  contributor  to  the  high  frequency,  high  amplitude  pressure 
fluctuations  (Her,  1986). 

Theoretical  u’ork  in  the  past  has  mainly  been  confined  to  semi-empirical  turbulence  models 
which  are  strongly  dependent  on  experimental  data  for  closure  of  the  problem.  These  models 
have  resulted  in  only  marginal  success.  Recent  numerical  solutions  to  the  full  time-dependent 
Navier-Stokes  equations  with  no  closure  assumptions  have  had  modest  success  (Handler, el  al. 
1984).  Several  characteristic  features  such  as  root-mean-square  (rms)  wall  pressure,  convec¬ 
tion  velocity  and  spectral  content  of  the  pressure  field  have  compared  favorably  to  experimental 
data.  This  numerical  work  appears  promising,  both  in  terms  of  reproducing  measured  proper¬ 
ties  of  the  wall  pressure  field  as  well  as  giving  insight  into  physical  features  of  this  field  which 
are  not  amenable  to  experimental  measurements. 

A  major  difficulty  encountered  in  evaluating  recent  numerical  methods  is  the  absence  of 
detailed  experimental  studies  of  the  wall  pressure  field  for  a  generic  flow.  The  small  body  of 
existing  wall  pressure  data  is  for  high  Reynolds  number  external  flows  which  are  typically  two 
orders  of  magnitude  higher  than  can  be  achieved  in  the  present  calculations.  In  addition,  these 
previous  experiments  are  poorly  documented  and  contain  measurements  of  only  some  of  the 
relevant  quantities. 

The  purpose  of  thi6  investigation  is  to  conduct  a  comparison/analysis  of  two  wall  pressure 
fluctuation  experiments;  one  computational  simulation  experiment  and  one  physical  labora¬ 
tory  experiment  that  are  of  the  same  order  magnitude  in  Reynolds  number.  For  the  first  time, 
detailed  measurements  of  the  wall  pressure  field  in  a  low-noise,  low  Reynolds  number,  fully  de¬ 
veloped  turbulent  channel  flow  will  be  compared  to  the  results  obtained  from  Direct  Numerical 
Simulations  (DNS)  of  the  Navier-Stokes  equations. 
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In  this  introduction,  an  attempt  will  be  made  to  review  the  more  salient  characteristics  and 
important  results  of  past  experiments  on  wall  pressure.  In  particular,  efforts  will  be  made  to 
point  to  those  classical  studies  which  are  considered  as  standards  in  the  literature.  In  addition, 
a  brief  synopsis  of  theoretical  and  computational  work  done  to  date  will  be  given,  including 
Direct  Numerical  Simulations  (DNS)  as  it  exists  at  this  time.  Finally,  the  specific  objectives 
of  this  study  will  be  put  forth. 


1.2  Review  of  Wall  Pressure  Studies 


Kraichnan  (1956  a,b)  performed  probably  the  earliest  complete  analysis,  in  a  qualitative 
sense,  of  the  important  features  of  the  wall  pressure  fluctuation  structure  within  boundary 
layer  flow  over  a  smooth  fiat  plate.  His  reasoning  and  dimensional  arguments  for  scales  and 
orders  of  magnitude  of  various  parameters  has  set  the  foundation  for  much  of  the  present  work. 
Combining  the  equation  of  continuity  with  the  divergence  of  the  incompressible  Navier-Stokes 
equations,  he  derived  a  Poisson  equation  for  the  unsteady  pressure  in  a  turbulent  field  which 
was  a  function  of  the  velocity  fluctuations  throughout  that  field.  At  this  point,  Kraichnan 
resorted  to  an  analysis  based  on  similarity  arguments  and  orders  of  magnitude,  since  detailed 
measurements  of  the  fluctuating  velocity  field  and  related  gradients  were  not,  at  the  time, 
well  documented.  He  showed  that  if  the  equation  for  pressure  was  written  in  integral  form,  a 
relation  for  wall  pressure  was  obtained  which  showed  the  fluctuating  pressure  to  be  coupled  to 
velocity  fluctuations  throughout  the  flow  region  being  investigated,  with  those  effects  closest 
to  the  wall  contibuting  the  major  share. 

Kraichnan ’s  work  on  the  theoretical  estimate  of  the  rms  value  of  the  wall  pressure  in  re¬ 
lation  to  the  velocity  fluctuations  via  similarity  arguments  was  extremely  innovative.  Further¬ 
more,  he  was  the  first  to  obtain  values  for  these  ratios  and  speculate  that  they  were  Reynolds 
number  dependent.  He  determined  that  the  ratio  of  rms  pressure  to  the  wall  6hear  stress  should 
be  of  the  same  order  of  magnitude,  approximately  equal  to  six,  and  that  for  the  same  class  of 
similarity  flows,  this  ratio  should  be  a  constant.  One  of  the  more  interesting  observations  of 
Kraichnan ’s  analysis  was  that  the  integration  of  the  correlation  function  over  the  wall  for  all 
time  delay  should  vanish,  and  therefore,  the  total  mean  square  force  per  unit  area  exerted  on 
the  wall  should  be  zero.  This  observation  is  still  open  to  discussion  within  the  literature. 

Willmarth  (1959)  performed  one  of  the  first  experiments  on  wall  pressure  within  a  flat 
plate  boundary  layer.  He  speculated  that  there  were  two  basic  mechanisms  for  the  generation 
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of  noise  by  turbulent  boundary  layers  in  general.  In  relation  to  rigid  walls,  he  determined 
that  a  rapid  change  of  fluid  properties  by  the  turbulent  flow  field  caused  noise  generated  by 
turbulence  only.  On  the  other  hand,  for  more  flexible  walls,  the  wall  acted  as  a  diaphragm  or 
speaker  and  hence  coupled  with  the  low  wave  number  harmonics  of  turbulence  which  amplified 
this  sound  and  radiated  to  the  farfield.  The  subsequent  generation  of  sound  may  be  a  resonance 
phenomena  involving  moving  waves  and  ripples  in  the  wall  which  are  produced  and  amplified 
by  a  convected  pattern  of  pressure  fluctuations  within  the  boundary  layer.  Willmarth  found 
that  the  high  frequency  end  of  the  spectrum  was  attenuated  due  to  the  spatial  limitations  of  the 
large  transducers  available  at  that  time  ( dt/6 “  =  1.1  to  4.6).  He  found  Prmt/Qoo  —  .00425  but 
as  dt/6*  decreased  to  the  order  of  1.1,  he  obtained  a  value  approaching  .0055.  Willmarth  further 
proposed  that  in  the  limit  as  the  transducer  diameter  approached  zero,  this  ratio  should  reach 
a  constant  value  of  0.006  and  be  independent  of  Reynolds  number,  contrary  to  Kraichnan's 
findings.  Willmarth  further  suggested  that  the  turbulent  pressure  fluctuations  lost  half  their 
identity  in  traveling  two  boundary  layer  thicknesses  and  almost  all  of  it  by  10 6.  In  agreement 
with  Kraichnan’s  work,  Willmarth  suggested  that  one  can  not  compare  pressure  fluctuations 
to  those  of  velocity  in  terms  of  amplitudes  or  frequency  content,  since  the  pressure  fluctuations 
depend  on  an  integral  relationship  with  the  velocity  throughout  the  boundary  layer. 

As  interest  began  to  develop  in  the  wall  pressure  field,  it  became  apparent  that  at  this 
early  stage  in  experimentation,  state-of-the-art  pressure  transducers  lacked  resolution  capability 
due  to  their  large  diameters.  In  addition,  as  mentioned  earlier,  turbulence  acted  as  a  source 
mechanism  of  unwanted  acoustic  sound  which,  in  most  cases,  contaminated  the  low  end  of  the 
frequency  spectrum.  In  1962  Willmarth  and  Wooldridge  conducted  one  of  the  first  extensive 
investigations  of  wall  pressure  on  a  flat  plate,  wherein  care  was  taken  to  investigate  transducer 
size  effects  and  minimize  background  contaminating  acoustic  noise.  Their  ratio  of  boundary 
layer  thickness  to  transducer  diameter  was  6/dt  -  30,  one  of  the  smallest  transducers  used  until 
that  time.  They  obtained  ratios  of  Prmt/Twti  *nd  prmt/Qoo  of  2.64  and  .00561  respectively. 
They  found  only  a  slight  Reynolds  number  effect  in  the  ratio  of  pr m,  to  rwau  which  decreased 
as  the  Reynolds  number  increased.  Of  course,  with  increasing  Reynolds  number,  the  frequency 
content  of  the  wall  pressure  would  extend  to  a  higher  frequency  range.  Since  the  transducer 
diameter  was  fixed,  more  attenuation  would  be  apparent  at  the  higher  frequencies,  so  this  result 
was  not  totally  unexpected. 

Willmarth  and  Wooldridge  addressed  the  issue  of  low  frequency  acoustic  contamination  by- 
cutting  off  the  spectrum  at  that  frequency  where  the  contamination  began,  and  assumed  a  con¬ 
stant  spectrum  amplitude  for  all  lower  frequencies.  They  made  combined  spatial  and  temporal 
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cross-correlation  measurements  which  showed  that  pressure  fluctuations  were  convected  down¬ 
stream.  These  measurements  were  hi-pass  filtered,  hence  low  frequency,  turbulence-related 
pressure  data  w*as  eliminated  from  their  results.  In  any  case,  they  determined  that  the  con¬ 
vection  speed  of  a  pressure  producing  eddy  varied  from  approximately  0.56  to  0.83  and  that 
eddies  of  a  certain  size  decayed  on  the  order  of  four  to  six  boundary  layer  thicknesses.  The  large 
convection  velocity  suggested  that  significant  contributions  came  from  the  large  scale  motion 
in  the  outer  region.  Finally,  Willmarth  and  Wooldridge  found  that  longitudinal  and  trans¬ 
verse  macro  scales  were  of  the  same  order  of  magnitude  and  suggested  that  for  future  studies, 
the  transducer  should  be  of  the  order  of  dt  =  Uoo/2Fmax  (where  Fmax  equals  the  maximum 
frequency  for  which  turbulence  information  is  present  in  the  flow)  in  order  to  prevent  high 
frequency  attenuation. 

Corcos  (1962)  presented  a  quasi-theoretical  scheme  which  addressed  the  issue  of  finite 
resolution  of  the  high  frequency  content  of  the  wall  pressure  field.  Using  the  results  of  Willmarth 
and  Wooldridge  (1962),  Corcos  derived  functional  relationships  which  allowed  all  statistical 
averages  related  to  the  6pace-time  covariance  of  the  wall  pressure  field  to  be  systematically 
corrected  for  finite  size  transducers.  These  corrections  were  found  to  be  a  very  sensitive  function 
of  the  convection  velocity,  generally  non-linear  and  to  a  large  extent,  over  corrected  in  those 
regions  which  had  been  highly  attenuated. 

Bakewell,  et  al.(  1963)  conducted  one  of  the  first  comprehensive  experiments  in  a  fully 
developed  pipe  flow.  They  applied  Corcos ’s  (1962)  corrections  to  all  data  but  found  minimal 
effect  even  though  their  ratio  of  pipe  radius  to  transducer  diameter  was  only  25  : 1,  which  was 
slightly  smaller  than  Willmarth  and  Woodridge’s.  Over  a  broad  band  analysis,  Bakewell  found 
the  convection  velocity  to  be  Ue/Uoo  —  0.7  and  appeared  to  be  somewhat  Reynolds  number 
dependent.  On  the  other  hand,  for  a  narrow  frequency  band,  e.g.,  low-passed  at  500 Hz,  only 
large  structures  were  present  moving  at  close  to  the  maximum  velocity  hence,  Ue/Ue 0  =  0.98. 
For  a  high  frequency  band  represented  by  hi-passing  at  7000Hz ,  only  the  small  scale  (high 
frequency)  structure  was  present,  and  a  much  lower  convection  speed  was  obtained  on  the 
order  of  0.63.  The  rms  value  of  the  wall  pressure  was  computed  by  averaging  the  time  history 
data  from  sets  of  octave  bands  and  found  Prm»/Q oo  —  0.0048  to  0.0063  which  was  somewhat 
Reynolds  number  dependent. 

At  that  time,  an  interest  was  beginning  to  develop  in  the  actual  structure  of  the  wall 
pressure  field  and  whether  it  scaled  with  definable  similarity  variables.  Corcos  (1963)  reviewed 
a  large  compendium  of  work  done  to  that  date,  with  some  attempt  at  explaining  similarities 
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and  differences  in  the  results.  He  substantiated  that  prmt/Twaii  was  weakly  Reynolds  number 
dependent  and  that  length  scales  within  the  inner  region  of  a  boundary  layer  seemed  to  be 
impressed  by  the  outer  flow.  Corcos  hypothesized  that  the  pressure  field  was  not  as  fine-grained 
as  the  velocity  field  and  that  the  spectrum  should  peak  at  a  value  of  u6m/U<X)  =  0.2  to  0.3.  He 
found  that  mixed  inner /outer  variables  for  the  magnitude  and  outer  variables  for  the  frequency 
appeared  to  collapse  all  Reynolds  number  flows  into  coincidence.  However,  he  further  theorized 
that  this  was  most  likely  due  to  the  cancellation  effects  of  finite  size  transducers.  Corcos  felt 
that  these  differences  at  high  frequencies  were  not  resolvable  and  that  similar  differences  at 
the  low  frequency  end  of  the  spectrum  were  due  to  secondary  flows  in  the  free  stream  and/or 
radiated  noise. 

Although  other  authors  had  made  cross-spectral  density  measurements,  Corcos  pointed 
to  the  importance  of  this  function  as  indicative  of  the  true  convective  nature  of  the  flow  and 
that  measurement  of  the  cross-spectrum  phase  quantified  a  narrow-band  measure  of  convection 
velocity  t/e(w,£).  Corcos  finally  concluded  that  the  lack  of  similarity  between  various  Reynolds 
number  flows  was  due  to  the  fact  that  as  pressure  patterns  convect  downstream,  evidence  seems 
to  suggest  that  instead  of  becoming  incoherent,  the  dispersive  nature  of  the  source  terms  causes 
the  pattern  to  spread  out  and  become  dispersive  instead. 

Willmarth  and  Roos  (1965)  made  the  assumption  that  if  the  transducer  responds  linearly  to 
the  pressure,  then  a  complete  and  accurate  set  of  wall-pressure  measurements  can  be  corrected 
for  transducer  size  effects.  They  tested  four  different  6ize  transducers  in  the  range  0.122  < 
dt/6m  <  0.442.  They  found  that  as  the  transducer  diameter  approached  zero,  the  ratio  of  prm, 
to  Twau  approached  the  limiting  value  of  prm»/Twaii  =  2.66  for  Rt  =  38,000.  They  also  found 
that  the  similarity  corrections  of  Corcos  (1962)  were  not  valid  when  u >6*/Uoo  >  3  for  spatial 
separations  such  that  Ax  <  0.16,  and  that  there  existed  departures  even  as  low  as  u>6m/Uoo  —  1. 

Schloemer  (1966)  investigated  pressure  gradient  effects  on  wall  pressure  fluctuations  along 
a  flat  plate.  His  conclusion  was  that  convection  velocities  were  changed  due  to  the  distortion 
of  the  mean  velocity  profiles,  which  were  caused  by  imposed  pressure  gradients.  With  respect 
to  zero  pressure  gradient  flows,  UeIUoo  was  higher  in  the  favorable  gradient  and  lower  in  the 
adverse  case.  Schloemer  determined  that  loss  of  coherence  or  decay  in  the  flow  direction  was 
more  rapid  in  the  adverse  than  in  zero  pressure  gradient  case  and,  as  expected,  the  favorable 
pressure  gradient  demonstrated  slower  rates  of  decay.  The  most  important  result  of  his  work 
was  that  he  was  one  of  the  first  to  establish  or  hypothesize  that  convection  velocity  was  a 
function  of  some  non-dimensional  frequency  based  on  boundary  layer  thickness.  He  also  argued 
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that  smaller  eddies  at  higher  frequencies  must  “physically”  exist  closer  to  the  wall  in  lower  speed 
regions,  therefore,  Ue/Uoc  should  drop  to  zero  as  the  velocity  goes  to  zero. 

Bull  (1967)  concluded  that  the  wall  pressure  field  has  a  structure  produced  by  contributions 
from  turbulent  sources  in  the  boundary  layer  with  a  wide  range  of  convection  velocities  and 
therefore  comprise  two  families  of  convected  wave  number  components.  The  first,  at  high  wave 
number,  is  associated  with  turbulent  motion  in  the  constant  stress  layer  and  is  coherent  for 
times  proportional  to  the  time  to  be  convected  a  distance  6  —  A.  The  second  family,  composed 
of  A’s  greater  than  twice  the  boundary  layer  thickness  lose  coherence  as  a  group  more  or  less 
independent  of  wavelength  and  are  associated  with  large  6cale  eddy  motion  outside  the  constant 
stress  region.  By  referring  to  these  definitions,  Bull  could  address  the  issue  of  where  in  the  flow 
field  the  pressure  source  terms  were  likely  to  be  located. 

Both  Bakewell  (1968)  and  Clinch  (1969)  investigated  wall  pressure  fluctuations  in  water 
flow.  Bake  well’s  experiment  was  on  the  exterior  surface  of  a  cylinder  for  which  in  the  flow 
direction,  agreement  with  flat  plate  data  was  obtained.  Only  in  the  spanwise  statistical  prop¬ 
erties  did  curvature  effects  become  dominant.  Clinch  investigated  fully  developed  pipe  flow  and 
found  Prmi/Qoo  to  be  independent  of  Reynolds  number,  on  the  order  of  0.007.  Clinch  found 
macro  and  micro  length  scales  in  the  longitudinal  direction  to  exceed  those  in  the  transverse 
direction  by  6  and  1.5,  respectively. 

Blake  (1970)  made  wall  pressure  measurements  on  a  flat  plate  with  both  smooth  and  rough 
walls.  In  order  to  address  the  issue  of  transducer  resolution,  his  experiment  was  the  first  to 
use  small  pin-hole  microphones  which  provided  resolution  approximately  three  times  better, 
relative  to  displacement  thickness,  than  had  been  previously  accomplished.  He  determined  that 
the  6pace-time  decay  rate  was  considerably  higher  than  previously  reported.  Blake  found,  as 
Schloemer(1966)  first  noted,  that  convection  velocity  for  small  separations  dropped  off  sharply 
at  low  frequencies.  Blake  gave  three  definitions  of  convection  velocity.  If  one  draws  a  curve 
tangent  to  the  space-time  correlation  curves  then: 

1.  Ue(r\)  =  dri/dr  —  instantaneous  broad  band  convection  velocity, 

and 

2.  Ue(r i)  »  fj/r  —  average  broad  band  convection  velocity. 
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Blake’s  third  definition  was  a  narrow  band  convection  velocity  defined  from  the  cross-spectral 
density  function,  where  Ve  =  -rjw/o,  which  was  in  agreement  with  Corcos’s  definition  (1963). 
Blake  further  found  that  the  low  frequency  pressure  spectra  scaled  on  variables  mixed  from 
the  inner  as  well  as  outer  flow,  but  that  the  smaller  scale,  high  frequency  components  scaled 
only  on  the  inner  variables  /  a  v/u*  and  t>  «  tt*.  Blake’s  results  were  probably  the  most 
definitive  work  done  in  which  great  care  was  taken  to  minimize  flow  disturbances  and  in  which 
an  attempt  was  made  to  improve  transducer  resolution  by  an  order  of  magnitude.  Indeed,  even 
to  this  day,  Blake’s  results  for  flat  plate  flow  are  still  used  as  baseline  data  for  comparative 
purposes. 

Willmarth  and  Yang  (1970)  conducted  a  comparison  experiment  between  the  flow  on  a 
flat  plate  to  that  external  to  a  cylinder.  They  found  that  longitudinal  and  transverse  scales 
were  less  for  the  cylinder  than  a  plane  boundary  layer  which  they  interpreted  as  the  transverse 
curvature  reducing  the  size  of  the  pressure  producing  eddies.  In  addition,  they  determined 
Prm»/Twaii  to  approach  2.5  for  both  experiments. 

Wills  (1970)  made  measurements  of  the  wave-number /phase  velocity  spectrum.  Up  to  this 
point  in  time,  most  experimentalists  had  evaluated  pressure  data  by  computing  the  space- time, 
correlation  function  and  had  integrated  (Fourier  transformation)  to  obtain  the  pressure 
spectrum.  Wills  chose  to  integrate  twice  (ie.,  with  multiple  transducers)  to  obtain,  for  the  first 
time,  the  full  (wave  number,  frequency)  spectrum.  Thus  one  could  look  at  a  contour  map  for 
which  a  definition  of  Ue(k)  was  depicted  explicitly.  Wills  found  that  convection  velocity  at 
low  wave  number  (low  frequency)  was  as  high  as  0.9Uoo  and  asymptoted  to  0.551/qo  at  the 
largest  wave  number  components.  His  data  also  showed  that  at  lower  wave  numbers,  there  was 
a  “decreasing”  trend  of  the  convection  velocity  indicated,  thereby  collaborating  Schloemer  and 
Blake’s  earlier  work  with  pin-hole  transducers.  Wills’s  method  provided  a  means  for  making 
corrections  both  for  transducer  effects  as  well  as  low  frequency  contaminating  noise. 

Willmarth  (1975)  gave  a  comprehensive  review  of  the  extensive  effort  towards  investigating 
wall  pressure  fluctuations.  He  outlined  that  it  was  not  possible  to  measure  pressure  fluctuations 
within  the  flow  field  due  to  the  intrusive  nature  of  the  probe  and  that  at  the  wall,  resolution 
problems  were  still  the  rule  rather  than  the  exception.  He  further  showed  that  the  pressure  at 
any  one  point  is  produced  by  momentum  fluctuations  at  many  otheT  points  and  that  theoretical 
work,  up  to  that  point,  mostly  employed  turbulence  models  based  on  experimental  work. 
Summarizing  the  efforts  of  many  experimentalists,  Willmarth  suggested  that  the  sound  field 
at  low  frequencies  contaminates  most  data  and  is  difficult,  if  not  impossible,  to  eliminate  and 
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that  although  pin-hole  effects  were  still  not  well  known  or  understood,  he  commended  Blake's 
(1970)  results  as  then  the  best  resolvable.  Willmarth  pointed  out  that  other  than  Schloemer. 
Blake  and  Wills,  experimenters  for  the  most  part  were  calculating  convection  velocity  from 
correlation  curves  rather  than  from  the  cross-spectral  density  functions;  ».e.,  a  broad-band 
analysis  only. 

From  the  review  by  Willmarth,  it  became  apparent  that  the  two  overriding  important 
themes  for  making  high  quality  measurements  of  the  wall  pressure  held  was  first  to  eliminate 
the  low  frequency  contaminating  acoustic  sound  and  secondly,  placing  a  more  careful  emphasis 
on  obtaining  smaller  transducers.  During  this  past  decade,  technological  advances  in  piezo¬ 
electronics  has  provided  much  smaller  transducers  and  a  number  of  new  noise  cancellation 
techniques  were  combined  with  facility  redesign  to  reduce  or  eliminate  the  effects  of  the  facility¬ 
generated  noise. 

Willmarth  and  Yang  (1970)  made  one  of  the  earliest  attempts  at  reducing  facility-generated 
noise  from  interfering  with  low  frequency  wall  pressure  content.  Willmarth  showed  that  a  finite 
correlation  was  obtained  for  negative  time  delay  between  flush-mounted  transducers  separated 
in  the  streamwise  direction.  By  combining  this  temporal  quantity  with  the  spatial  separation 
between  the  two  sensors,  a  velocity  was  derived  equal  to  the  sr  and  speed  (relative  to  the 
flow  velocity)  in  the  fluid.  Willmarth  interpreted  this  quantity  as  flow-generated  acoustic 
noise  and  subtracted  it  from  the  original  output  to  eliminate  its  contribution  to  the  overall 
signal.  However,  this  technique  is  somewhat  subjective,  and  in  reality,  can  only  be  applied  to 
correlation  measurements  when  it  is  known  a  priori  that  the  low  frequency  noise  is  composed 
of  a  single  frequency  component. 

A  new  technique,  which  employs  a  signal  processing  scheme,  was  introduced  first  by  Wamb- 
sganss  and  Zaleski  (1970),  and  later  refined  by  Wilson  et  al. (1979)  and  Horne  and  Hansen 
(1981).  Utilizing  three  pressure  transducers  located  in  a  plane  perpendicular  to  the  flow  direc¬ 
tion,  simple  subtraction  of  the  time-dependent  pressure  signals  provided  a  means  of  separating 
the  individual  signals  due  to  turbulent  wall  pressure,  background  acoustic,  and  electronic  noise, 
respectively.  The  method  may  be  applied  to  wider  band  noise  but  relies  on  the  belief  that  the 
low  frequency  noise  only  propagates  as  a  plane  wave.  Wilson  et  al.  investigated  air  flow  in 
an  annular  pipe  flow  configuration  at  a  Reynolds  number  of  300, 000.  The  turbulent  pressure 
signal  was  some  20 dB  less  than  their  overall  signal  from  a  single  transducer,  resulting  in  the 
subtraction  of  very  large  numbers  to  obtain  a  much  smaller  one.  Their  values  for  Prm$/Q<» 
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ranged  from  .005  to  .0062  also  indicating  a  Reynolds  number  dependence.  Although  an  innova¬ 
tive  technique,  this  method  requires  the  noise  to  propagate  as  a  plane  wave  which  is  somewhat 
restrictive.  Furthermore,  as  will  be  shown  later  in  this  study,  the  subtraction  scheme  involves 
a  major  source  of  potential  error  (although  subtle)  which  involves  “all”  signal  processing  sub¬ 
traction  techniques. 

Panton  et  a/.(1980)  investigated  the  low-frequency  region  of  the  wall  pressure  spectrum 
on  the  external  fuselage  of  a  sailplane.  The  zero-pressure  gradient  flow  had  low  free-stream 
turbulence  intensity  and  showed  the  pressure  spectra  to  decrease  with  decreasing  frequency  in 
agreement  with  calculations  relative  to  a  turbulence-mean  shear  interaction  mechanism.  The 
measurements  showed  virtually  no  acoustic  contamination  at  low  frequency  and  were  among  the 
first  to  show  an  approximate  “u>2”  dependence  with  rising  frequency.  An  important  point  made 
by  Panton,  was  how  the  pressure  fluctuations  act  as  a  foot-print  of  the  velocity  fluctuations 
throughout  the  boundary  layer,  thereby  providing  a  convenient  mechanism  by  which  the  flow 
structure  in  the  outer  region  can  influence  events  going  on  within  the  inner  layer  close  to  the 
wall. 


Schewe  (1983)  investigated  the  resolution  problem  inherent  with  high  frequency  wall  pres¬ 
sure  measurements  by  testing  small  (non-pinhole)  flush  mounted  transducers  ranging  from 
25.4mm  to  1mm  in  diameter.  These  sizes  correspond  to  viscous  diameters  (rf+  =  dtum  fv)  of 
d+  =  19  to  333  and  are  considered  to  be  the  most  definitive  work  on  transducer  resolution  to 
date.  Schewe  showed  that  as  the  transducer  diameter  decreased,  the  probablility  distribution 
function  exhibited  considerable  expansion.  In  other  words,  not  only  does  the  larger  diameter 
transducer  filter  off  the  higher  frequency  (smaller  structure)  components  of  the  wall  pressure, 
but  it  filters  out  very  large  pressure  amplitudes  which  contribute  significantly  to  the  rms  level 
of  pressure  but  have  short  life  times.  For  the  smallest  transducer,  (rf+  =  19),  Schewe  found 
Prm»/Qoo  —  0.0102.  However,  he  used  Corcos’6  (1962)  corrections  to  adjust  all  data  with  a 
five  percent  error  in  the  smallest  transducer.  If  one  compares  his  spectral  data  (scaled  on 
inner  variables)  with  the  work  of  Blake  (1970)  and  Emmerling  (1973),  it  appears  that  the  low 
frequency  part  of  the  spectrum  is  somewhat  low.  In  other  words,  all  the  transducers  appear 
to  be  somewhat  attenuated  at  low  frequency  which  could  be  indicative  of  a  slight  calibration 
problem.  Nevertheless,  Schewe’s  resolution  of  the  high  end  of  the  spectrum  was  pursued  with 
great  care  and  may  be  regarded  as  a  standard  for  zero-pressure  gradient  wall  bounded  flows. 

Farabee  (1986),  investigated  wall  pressure  fluctuations  beneath  a  non-equilibrium  turbu¬ 
lent  flow.  He  conducted  experiments  with  pinhole  transducers  similar  to  those  used  by  Blake, 
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in  a  flow  passing  over  a  backward-facing  step.  The  objective  of  such  a  study  was  to  identify 
those  turbulent  structures  in  the  flow  that  contribute  to  the  source  terms  of  the  wall  pressure 
field.  Data  was  taken  for  both  the  perturbed  and  non-perturbed  states  for  comparison  pur¬ 
poses.  The  high  frequency  wall  pressures  were  found  to  scale  on  inner  variables  suggesting  that 
the  sources  for  this  frequency  range  can  be  found  within  the  inner  layer  of  the  boundary  layer 
close  to  the  wall.  The  maximum  in  the  non-perturbed  boundary  layer  was  found  to  occur  at 
uS’/Uoo  =  0.3.  The  extensive  two-point  statistics  were  in  excellent  agreement  with  the  work 
of  Blake,  Bull  and  Emmerling,  and  provided  6ome  guidance  for  the  present  work. 


1.3  Theoretical  Studies  and  Numerical  Simulation 


For  the  most  part,  analytical  studies  have  attempted  to  understand  the  relationships  ex¬ 
isting  between  the  different  scales  of  turbulence  and  how  they  induce  and/or  affect  pressure 
fluctuations.  As  illustrated  above,  wall  pressure  fluctuation  experiments  have  been  more  or  less 
limited  to  examining  the  statistical  properties  of  the  pressure  field  and  by  conjecture,  relate 
them  to  turbulence  structures  known  to  exist.  This  level  of  understanding  has  for  the  most 
part,  been  limited  to  equilibrium  type  boundary  layers  for  high  Reynolds  numbers  and  very 
little  evidence  exists  tying  shear  layers  and  internal /external  wall  flows  together  in  a  balanced 
theory. 

Following  the  pioneering  work  of  Kraichnan  (1956),  Landahl  (1967)  attempted  to  numeri¬ 
cally  solve  the  nonhomogeneous  Orr-Sommerfeld  equation  for  a  turbulent  boundary  layer  and 
assumed  that  the  least  damped  eigenvalues  of  the  problem  would  be  the  dominant  sources  of 
the  wall  pressure  field.  His  results  for  the  pressure  cross- spectral  density  for  a  flat  plate  bound¬ 
ary  layer  agreed  favorably  with  Willmarth  and  Wooldridge’s  experimental  results.  Landahl’s 
calculated  streamwise  decay  rates  agreed  well  with  the  experimental  results  of  Bull  (1967). 
However,  the  calculated  convection  velocities  were  found  to  be  too  low  by  30  percent. 

Panton  and  Linebarger  (1974)  calculated  the  one-dimensional  wavenumber  (streamwise 
wavenumber)  spectrum  of  pressure  fluctuations  beneath  equilibrium  turbulent  boundary  layers. 
They  numerically  solved  the  Poisson  equation  for  wall  pressure,  following  the  development  of 
Kraichnan,  by  using  empirical  models  for  the  turbulence  velocity  source  terms  in  the  equation. 
An  interesting  feature  of  their  work  is  that  they  ho  wed  the  contributions  to  the  total  pressure, 
as  a  function  of  wavenumber,  that  come  from  various  regions  of  the  boundary  layer.  The  inner 
portion  of  the  boundary  layer  dominates  the  contributions  at  the  highest  wavenumbers.  When 


11 


1 


contributions  come  entirely  from  the  wall  layer  and  log  region,  a  universal  overlap  region  in 
the  pressure  spectrum  was  shown  to  exist.  The  presence  of  an  overlap  region  can  be  argued 
following  the  same  dimensional  arguments  that  are  used  to  show  the  presence  of  an  overlap 
region  for  the  mean  flow  velocity  profile.  Within  this  overlap  region,  the  spectrum  should  scale 
as  u>-1.  Interestingly,  the  contributions  to  the  low  wavenumber  components  come  from  the 
wake  region,  yet  the  contribution  to  the  lowest  wavenumber  components  come  from  the  middle 
region  of  the  boundary  layer. 

Chase  (1980)  '•('"ducted  an  extensive  theoretical  work  attempting  to  model  the  wall- 
pressure  field  util./'  _  'nown  facts  extracted  from  detailed  experimental  work.  His  objective 
was  to  provide  a  desc  ,;tian  °f  the  pressure  spectrum  that  would  be  valid  for  both  the  convec¬ 
tive  (high  wavenumoer)  components  as  well  as  for  low  wavenumber  components.  Chase  used 
the  Poisson  equation  formulation  given  by  Kraichnan  and  modeled  the  fluctuation  velocities  as 
propagating  waves,  as  suggested  by  Morrison  and  Kronauer  (1969).  An  analytical  expression 
for  the  pressure  spectrum  was  obtained  containing  adjustable  constants  that  were  evaluated 
from  experimental  results.  Chase's  model  agreed  well  with  experiments,  predicted  the  V'2r 
dependence  of  the  pressure  spectrum  at  low  frequencies,  and  deduced  that  the  phase  convection 
velocity  should  tend  to  zero  at  low  frequencies. 

Ffowcs- Williams  (1982)  used  the  similarity  model  of  Corcos  (1962)  to  model  the  turbulent 
sources  of  pressure,  not  the  pressure  itself,  so  that  a  Corcos-like  description  of  the  pressure 
spectrum  could  be  derived  that  was  consistent  with  the  constraints  imposed  by  the  governing 
equations.  The  relation  which  was  derived,  resulted  in  a  fairly  compact  specification  of  the 
pressure  spectrum  with  a  set  of  constants  similar  to  Chase’s  (1980)  model,  which  are  dependent 
on  experimental  results.  He  concluded  that  the  existing  body  of  literature  formed  an  inadequate 
base  for  determining  these  constants  explicitly  and  that  new  experiments  were  needed.  Ffowcs- 
William6  showed  that  the  surface  spectrum  was  proportional  to  the  square  of  frequency  at  low 
wavenumbers  with  subsonic  phase  speed,  Uc  =  u>/k.  He  further  showed  that  the  acoustics  of  the 
problem  determines  the  strength  of  all  those  low-wavenumber  elements  with  supersonic  phase 
velocities.  The  sound  radiated  by  turbulence  and  the  ability  to  discriminate  against  self-noise 
for  incoming  acoustic  waves  is  therefore  set  entirely  by  spectral  elements  whose  magnitude 
can  not  be  determined  experimentally  unless  the  experiment  is  conducted  in  a  facility  that  is 
acoustically  anechoic. 

Since  the  equation  for  pressure  at  the  wall  is  a  Poisson  type,  it  is  natural  to  expect 
the  evolution  of  various  attempts  to  solve  this  equation  to  take  the  path  as  indicated  above. 
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However,  all  such  methods  resulted  in  undefined  constants  that  must  rely  on  experimental 
measurements  to  explicitly  predict  the  pressure  spectra. 

Direct  numerical  simulations  (DNS)  of  turbulent  shear  flows  have  recently  provided  re¬ 
searchers  with  a  new  tool  that  has  been  used  to  study  the  dynamics  of  turbulence  in  greater  de¬ 
tail.  DNS  represents  three-dimensional,  time- dependent  numerical  solutions  of  the  full  Navier- 
Stokes  equations  that  attempt  to  compute  all  significant  scales  of  motion  without  using  any 
turbulence  models  for  closure.  In  some  cases,  these  calculations  have  provided  information  that 
would  be  difficult  or  impossible  to  obtain  from  experiments  and  in  other  cases,  provided  new 
results  later  confirmed  by  experiments.  As  a  review,  a  brief  synopsis  follows  of  the  principal 
efforts  that  have  been  undertaken  by  numerous  investigators  in  the  direct  numerical  simulation 
of  turbulence,  with  special  emphasis  on  a  simulation  of  wall  bounded  turbulence.  The  reader 
may  refer  to  reviews  by  Rogallo  and  Moin  (1984)  and  Schumann  and  Friedrich  (1987)  for  a 
more  exhaustive  review  of  Direct  Numerical  Simulation  (DNS)  techniques  and  results. 

Orszag  and  Patterson  (1972)  performed  the  first  simulation  of  homogeneous  isotropic  tur¬ 
bulence  which  required  “323”  grid  points.  A  spectral  Galerkin  approach  was  used  based  upon 
trigonometric  functions.  Orszag  and  Kells  (1980)  and  Orszag  and  Patera  (1980)  investigated 
turbulent  channel  flow  in  which  one  of  the  homogenous  coordinates  used  Chebyshev  polynomi¬ 
als  in  place  of  trigonometric  functions.  Orszag  and  Kells  investigated  the  problem  of  transition 
to  turbulence  and  Orszag  and  Patera  extended  this  work  to  a  full  turbulence  calculation. 

Marcus  (1984)  investigated  the  rotating  cylinder  problem  represented  as  Taylor- Couette 
flow  and  satisfactorily  predicted  wave  speeds.  Marcus’s  use  of  Green  functions,  or  so-called 
influence  matrix  methods,  were  an  improvement  to  DNS  by  a  more  acceptable  implementation 
of  the  necessary  boundary  conditions. 

Handler  et  al. (1984, 1988)  and  Hansen  et  al. (1987)  performed  the  first  analysis  of  the  wall 
pressure  field  from  a  DNS  calculation.  The  single  and  two-point  spectra  were  obtained  for  a 
channel  flow  Reynolds  number  of  R*  =  5000  with  a  32s  calculation.  The  dominant  features 
to  include  the  general  convective  nature  of  decaying  turbulence  was  successfully  demonstrated. 
However,  the  smallest  scales  of  motion  were  unresolved,  resulting  in  a  much  slower,  broadband 
decay  rate  then  experiments  have  shown.  They  speculated  that  increased  spatial  resolution 
would  improve  the  results  at  high  frequency.  They  were  the  first  to  suggest  that  aliasing 
errors  were  possible  due  to  spatial  filtering  when  transferring  back  and  forth  between  the 
Fourier  computational  space  and  the  “physical”  domains  as  is  typically  carried  out  in  direct 
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numerical  simulation  calculations.  They  further  speculated  that  both  direct  as  well  as  large 
eddy  simulations  (LES)  of  the  Navier-Stokes  equations  could  provide  the  means  of  studying 
and  interpreting  the  underlying  physical  processes  taking  place  within  the  turbulent  velocity 
field.  Thus,  evaluating  results  from  such  calculations  might  give  information  about  the  physics 
of  wall  pressure  generation  not  easily  obtained  within  the  laboratory. 

For  the  most  part,  investigations  involving  DNS  for  the  incompressible  case  result  in  a 
number  of  time  point  realizations  for  which  all  three  components  of  velocity  and  the  pressure  are 
solved.  Many  of  the  statistical  functions  of  turbulence  involving  velocity  have  been  investigated 
with  the  results  from  DNS.  However,  other  than  the  work  of  Handler  et  a/.,  no  detailed  analysis 
of  the  wall  pressure  field  has  evolved  from  any  of  the  calculations  to  date. 


1.4  Objectives  of  Research 


The  review  of  the  literature  reveals  that  a  number  of  questions  about  the  wall  pressure  field 
remain  unresolved.  It  has  been  shown  that  wall  pressure  fluctuations  appear  to  be  generated  by 
velocity  fluctuations  throughout  the  boundary  layer.  To  this  date,  it  has  not  been  substantiated 
how  events  within  the  flow  field,  such  as  burst-sweep  events  affect  these  fluctuations.  It  also 
has  not  been  resolved  whether  pressure  convected  patterns  are  similar  or  directly  related  to 
velocity  convected  patterns.  Taking  into  account  all  of  the  past  experimental  and  theoretical 
work,  some  specific  questions  could  be  asked.  Do  axial  and  transverse  length  scales  for  wall 
pressure  have  limiting  values  for  higher  Reynolds  number  flows?  Do  in  fact,  unique  length, 
velocity  and/or  time  scales  exist  for  both  internal  and  external  flows  such  that  all  wall  bounded 
flows  have  universal  features?  For  example,  it  has  been  shown  both  via  experiments  and  to 
a  certain  extent,  with  theoretical  modeling  that  eddies  of  a  specific  size  decay  over  a  distance 
proportional  to  their  size.  Do  external  and  internal  flows  exhibit  the  same  characteristics,  and 
what  are  the  Reynolds  number  effects? 

Our  current  understanding  of  these  fluctuations  is  specific  to  high  Reynolds  number,  exter¬ 
nal  wall-bounded  flows  and  free  shear  layers.  Anomalies  that  may  exist  at  lower  speeds  can  not 
be  addressed  realistically  by  examining  the  past  literature.  In  light  of  new  analytical  techniques 
(DNS),  studies  of  the  wall  pressure  field  and  turbulence/turbulence  and  turbulence/pressure 
interaction  are  now  being  pursued.  A  full  experimental  study  close  to  the  Reynolds  number 
range  of  a  typical  DNS  calculation  would  prove  invaluable  for  code  validation  as  well  as  for 
documentation.  Therefore,  a  generic  wall  pressure  experiment  would  prove  useful  for  critiquing 
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DNS  as  well  as  for  comparing  future  DNS  applications  to  wall  roughness  and  non-Newtonian 
applications. 


The  specific  objectives  of  this  research  are: 

o  to  extend  the  present  data  base  by  comparing  for  the  first  time,  two  wall  pressure  experi¬ 
ments;  one  physical  experiment  in  the  laboratory  and  one  computational  simulation 
on  a  high  speed  computer; 

o  to  document  detailed  measurements  of  the  wall  pressure  field  at  a  range  of  Reynolds 
numbers  (based  on  channel  half- width)  from  10,000  to  about  25,000  in  a  turbulent 
channel  flow; 

o  to  document  spectral  statistics  in  the  wavenumber  and  frequency  domains  for  both 
experiments,  where  possible; 

o  to  utilize  novel  signal  processing  techniques  in  order  to  examine  the  low  frequency 
range  of  the  pressure  spectrum  thus  contributing  to  our  fundamental  understanding 
of  the  mechanism  of  flow-induced  noise; 

o  to  conduct  a  direct  comparison  between  experiments  at  similar  Reynolds  numbers  that 
will  prove  useful  in  assessing  both  qualitatively  and  quantitatively  the  validity  of  the 
recent  theoretical  work; 

o  to  examine  scaling  laws  for  both  flows  to  address  their  range  and  validity,  and  to 
determine  whether  similarity  laws  have  common  features  for  both  internal  and 
external  flows. 
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CHAPTER  2 


EQUATIONS  OF  TURBULENCE 


2.1  Turbulence  Definitions  and  Corresponding  Equations 


Laminar  and  turbulent  flow  may  be  classified  as  one,  two  or  three  dimensional,  depending 
on  the  number  of  spatial  coordinates  needed  to  describe  the  flow.  In  the  steady  laminar  case, 
the  fluid  streamlines  exhibit  a  constancy  in  time  with  no  superimposed  fluctuations  or  random 
interplay  of  any  one  streamline  with  another.  For  unsteady  flows,  the  streamlines  become  time 
dependent,  but  still  remain  relatively  smooth  with  no  apparent  fluctuations  superimposed  on 
the  overall  flow  pattern. 

By  contrast,  turbulent  flow  is  depicted  by  disorder,  randomness  and  unsteady  fluctuations, 
and  this  added  complexity  precludes  any  simple  analysis  as  in  the  laminar  case.  Turbulent 
flows  are  predominant  over  external  surfaces  such  as  airfoils,  turbine  and  compressor  blades, 
ship  hulls  and  missiles;  and  through  internal  configurations  such  as  ducts,  pipes,  nozzels  and 
diffusers.  Beyond  certain  critical  Reynolds  numbers,  most  laminar  flows,  whether  external 
or  internal,  become  unstable  when  subjected  to  small  flow  disturbances  and  move  through  a 
transition  phenomenon  towards  the  more  prevalent  turbulent  flow  characterization. 

There  are  two  alternatives  in  proceeding  with  any  turbulent  flow  analysis.  The  first  involves 
a  practical  semi-empirical  analysis  which  attempts  to  predict  the  macro  properties  of  the  flow 
such  as  mean  velocity  and  temperature  profiles,  boundary  layer  thicknesses,  Bkin  friction  and 
heat  transfer  coefficients.  Once  this  information  is  extracted,  the  more  useful  properties  of 
engineering  significance  such  as  lift,  drag,  flow  separation  and  heat  transfer  may  be  determined. 
The  second  alternative  involves  a  statistical  theory  of  turbulent  correlation  functions  which  seek 
to  analyze  these  random  flow  fluctuations.  This  mode  of  solution  endeavors  to  discover  the 
structure  of  turbulence  in  a  scientific  knowledge- accumulating  fashion.  There  are  many  flow 
investigations  for  which  both  these  theories  of  solution  overlap  but  almost  always  at  a  rather 
advanced  frontier-research  level  of  fluid  mechanics.  As  a  general  rule,  engineering  studies 
involve  the  former,  for  the  method  of  analysis  lends  itself  nicely  in  producing  the  parameters 
of  engineering  significance,  whereas,  the  latter  method  is  the  predominant  form  of  analysis  in 
a  more  functional  investigation  of  the  origins  of  turbulence  to  include  its  relative  scales  and 
energy  content.  However,  one  of  the  most  complicated  and  important  aspects  of  experimental 
turbulence  is  that  quantities  that  are  usually  obtained  in  the  laboratory  are  measured  because 
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they  are  ’'easy’*  to  get  rather  than  of  fundamental  significance.  In  other  words,  the  important 
quantities  which  may  further  our  understanding  are  not  necessarily  easily  obtained. 

In  molecular  motion,  molecules  at  most,  only  affect  the  motion  of  their  immediate  neigh¬ 
bors.  On  the  other  hand,  within  turbulence,  the  motion  at  any  point  within  the  flow  affects 
motion  at  other  points  within  the  flow  via  the  pressure  field  and  therefore  is  only  “quasi- 
random”.  No  two  realizations  of  turbulent  flow  parameters  (i.e.,  velocity  or  pressure)  are  the 
same.  In  this  case,  statistical  descriptions  of  the  distributions  of  mean  values  are  all  that  make 
sense.  Finally,  one  hears  turbulence  characterized  as  “homogeneous”  and  “isotropic”.  These 
terms  do  not  have  the. same  or  similar  meanings.  The  first  means  a  particular  turbulent  flow  is 
invariant  or  independent  of  one  or  more  coordinates  in  space  and/or  time.  The  second  means 
that  turbulent  "scales”  are  invariant  in  space.  Various  size  scales  with  quasi-organized  motion 
or  convection  depending  on  size  may  exist  within  homogeneous  turbulence,  but  only  one  size 
scale  may  exist  in  isotropic  turbulence. 

The  general  relationships  for  a  generic  flow  will  now  be  presented,  and  when  appropriate, 
reduced  to  appropriate  forms  for  a  2D-channel  flow  environment. 


2.1.1  General  Equations  for  Turbulent  Flows 

In  the  solution  of  a  typical  turbulent  flow  problem,  the  basic  equations  of  mass,  motion  and 
energy  for  a  viscous  fluid  are  assumed  valid  at  any  instant  of  time.  These  must  be  transformed 
into  alternate  forms  by  averaging  over  a  time  interval  consistent  with  our  ability  to  study  such 
flows  experimentally.  This  time  averaging  is  really  a  necessary  consequence  of  our  inability  to 
track  the  instantaneous  fluctuations  of  the  real  turbulent  flow.  In  general,  let  any  flow  variable, 
G,  be  subdivided  into  mean  and  fluctuating  components,  G,  and  g ,  respectively  such  that  the 
following  definitions  hold  true. 

G  =  G  +  g  (2.1) 


G  =  lim  ^ 
T- oo  T 

jT G  dt 

Jo 

(2.2) 

1 

fT 

0  =  Urn  - 

T-»oo  1 

gdt 

Jo 

(2.3) 

1 

fT-  , 

0  *  lim  - 
T—  oo  T  j 

/  Gg  dt 
>0 

(2.4) 
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(2.5) 


G‘  +  g2  =  lim  ~  [*&  dt 

T-»oo  2  2o 

The  period  T  over  which  the  integration  is  performed  must  be  chosen  significantly  larger 
than  the  largest  period  of  the  fluctuating  component,  g ,  and  the  integration  refers  to  taking 
the  time  average  of  the  variable  in  question.  The  instantaneous  equations  of  motion  for  an 
incompressible,  Newtonian  fluid,  with  constant  viscosity,  can  be  written  in  Cartesian  tensor 
notation  as, 


W  W  _  _1  dP 
dt  *  dx  j  pdii  ^  dx2j 


(2.6) 


where  V  =  V(V1=U,V2=V,V3=W)  represents  the  velocity  vector,  P,  pressure,  p,  density,  and 
v,  the  kinematic  viscosity.  Equation  (2.6),  representing  the  Navier-Stokes  equations,  is  time 
dependent,  incompressible  with  constant  viscosity  and  consists  of  three  equations  for  the  four 
independent  variables  l\V,W  and  P.  If  we  take  the  divergence  of  equation  (2.6),  we  obtain  a 
Poisson  equation  for  the  pressure  as, 


V2P  = 


dVj  dVj 
dii  dij 


(2.7) 


One  may  note,  that  as  written,  equations  (2.6)  and  (2.7)  involve  the  instantaneous  variables 
V ,  V,  IP  and  P  which  themselves  are  implicit  functions  of  time,  t  ,  and  the  spatial  coordinates, 
x,y,  and  z.  We  will  leave  until  section  3.2,  a  description  of  the  difficulties  and  complexities 
involved  with  any  attempt  to  solve  the  above  two  equations.  Suffice  it  to  say,  that  only  in  the 
last  decade  with  the  advent  of  high  speed  computers  have  these  two  equations  been  amenable 
to  computational  solutions  directly,  and  only  at  extremely  low  Reynolds  numbers  for  very 
restrictive  flow  geometries. 


As  depicted  above,  all  four  variables  represent  pseudo  random  functions  of  both  space  and 
time,  and  the  equations  defining  them  usually  require  some  manipulation  in  order  to  make  their 
solution  more  tractable.  Following  Reynolds  (1895),  we  decompose  both  variables,  velocity  V' 
and  pressure  P,  according  to  the  definitions  of  a  mean  and  fluctuating  component  as  outlined 
in  equation  (2.1).  Paying  close  attention  to  the  rules  specified  in  equations  (2.2)  to  (2.5)  we 
substitute  V  =  V  +  v  and  P  =  P  +  p  into  equations  (2.6)  and  (2.7),  take  the  time  average,  and 
obtain  the  following  equations  for  the  velocity  and  pressure  respectively. 


dt  + 


l  dP_ 

p  dxi  V  dx2j 


1  d  . _ , 


V2P  = 


dVj  dVj  dv:  dv,  ~ 

^  dxidij  dxidxj 


(2.8) 

(2.9) 
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As  indicated,  we  have  introduced  nine  new  variables  (as  represented  by  the  last  terms  in  both 
equations)  into  an  already  complicated  set  of  equations  which  prevents  any  simple  solution, 
since  they  require  additional  knowledge  of  the  turbulence  structure.  There  are  no  a  priori 
analytical  representations  for  these  terms,  and  hence,  they  are  usually  modeled  empirically. 
These  new  variables  cause  what  has  been  termed  the  "closure”  problem  of  solving  the  time 
mean  Navier-Stokes  equations.  Two  relatively  simple  and  similar  empirical  models  are  based 
on  the  eddy  viscosity  and  mixing  length  hypotheses  based  on  a  1st  order  closure  method.  The 
reader  may  refer  to  any  text  on  advanced  fluid  mechanics  for  a  presentation  of  this  and  higher 
order  models  of  closure  for  various  geometric  configurations. 


2.1.2  Channel  Flow 


The  derivation  of  a  solvable  set  of  relations  for  the  internal  flow  through  a  channel  is 
relatively  simple  due  primarily  to  the  fact  that  geometry  prohibits  a  continuous  growth  of  a 
characteristic  thickness.  We  assume  that  x  is  in  the  direction  of  flow,  y,  normal  to  the  surface, 
and  z  lies  in  the  plane  of  the  channel  walls  which  extend  to  infinity.  If  this  channel  is  long 
enough,  there  is  no  further  growth  in  the  boundary  layer  as  it  stretches  to  the  centerline,  and 
we  would  expect  all  variables  and  derivatives  thereof,  to  be  independent  of  x  except  for  the 
mean  pressure.  The  only  component  of  mean  velocity  we  expect  is  U  and  the  flow  lies  entirely 
in  the  xy-plane,  with  no  variation  in  a  direction  normal  to  that  plane.  Recognizing  that  V  and 
W  are  zero,  it  can  be  shown  that  equation  (2.8)  reduces  to  the  following  set  of  equations  for 
steady  flow. 


0 


1  dP  d  , ,  vdPU 

p  dx  dy  Ut)  dy 2 


(2.10) 


\dP_ 
dy 


Integrating  equation  (2.11)  with  respect  to  y,  results  in, 


(2.11) 


P{*>V)  ,  77  _  Po{*) 

-  +  V*  St  - 

P  P 


(2.12) 


Since  v2  is  independent  of  x,  differentiating  equation  (2.12)  with  respect  to  x  shows  to  be 
equal  to  which  is  constant  in  magnitude,  otherwise  the  flow  would  accelerate.  Using  this 
information,  equation  (2.10)  may  be  integrated  once  with  respect  to  y  to  obtain, 


0  = 


V  dPo 
p  dx 


uv  + 


vdV 

dy 


(2.13) 
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Due  to  symmetry,  the  total  shear  -ut>  +  must  equal  zero  at  the  channel  centerline,  hence 
equation  (2.13)  reduces  to  the  well  known  relation  for  shear  velocity  and  shear-stress  at  the 
wall  as,  __ 

T„alt  =  pul  =  -h^  (2.14) 

In  the  inertial  sublayer,  it  can  be  shown  that  U(y)  must  satisfy  both  a  logarithmic  law  related 
to  inner  variables,  and  a  velocity  defect  law  in  the  core  region,  similar  to  the  flow  along  a  flat 
plate.  The  resultant  relations  are, 

Log  -  Layer  ^ill)  —  1.  /n  0^-)  +  a  (2.15) 

ti»  A  u 

Core  -  region  - ~ — —  =  —  In  )  +  b  (2.16) 

\im  A  n 

where  U0  is  the  centerline  velocity,  k  is  von  Karman’s  constant,  and  the  constants  ‘a’  and 
‘6’  are  finite,  yet  independent  of  Reynolds  number.  The  reader  may  refer  to  the  works  of 
Schlichting  (1955),  Townsend  (1956),  or  Tennekes  and  Lumley  (1972)  for  a  detailed  derivation 
of  these  results.  In  the  literature,  there  is  some  scatter  in  the  constants  V  and  ‘5’  which 
may  be  due  to  a  shortage  of  experiments  being  conducted  at  large  enough  Reynolds  number 
to  ensure  their  independence.  A  second  reason  for  the  scatter,  is  that  it  becomes  difficult  to 
ideally  match  parameters  between  one  facility  and  another,  even  to  the  extent  of  matching  day 
to  day  conditions  for  the  same  facility. 

An  attempt  may  be  made  to  solve  equation  (2.13)  provided  one  assumes  some  eddy  vis¬ 
cosity  or  mixing  length  model  for  the  UtJ  term  to  effect  closure.  Using  finite  differencing  or 
finite  element  techniques,  solutions  to  this  equation  have  been  shown  to  be  relatively  simple 
and  straightforward,  provided  an  accurate  set  of  initial  and  pressure  gradient  conditions  can 
be  found.  In  the  next  chapter,  we  will  present  the  methodolgy  for  solving  equations  (2.8)  and 
(2.9)  by  the  method  of  DNS  with  no  "closure"  assumptions.  We  proceed  now  to  review  the 
definitions  of  higher  order  statistics  of  a  time  dependent  function. 


2.1.3  Higher  Moments 

As  discussed  previously,  the  details  of  the  turbulent  velocity  field  can  not  be  reproduced 
instantaneously.  It’s  only  in  the  mean,  or  on  average,  that  the  mean  values  within  stationary 
time  independent  flows  can  be  obtained.  The  mean  value  (over  time)  of  velocity  is  the  first 
moment,  but  higher  moments  for  variance  <r2,  skewness  5,  and  the  kurtosis  or  flatness  factor 
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A',  can  also  be  defined.  Since  any  variable  has  been  defined  to  equal  the  sum  of  its  mean  value 
and  a  fluctuating  component  about  the  mean,  the  higher  moments  of  interest,  are  those  due 
to  the  fluctuations  themselves.  These  moments  serve  to  describe  the  amount  of  disorder  or 
randomness  present  and  characterize  the  turbulence  more  definitively.  All  the  moments  can  be 
defined  either  through  time  averages  or  the  probability  density  function  B( u),  as  illustrated 
below. 


romance 


(skewness)  S  -  =  lim  4  / 

<r3  T—oc  T  J o 


i  dt  =  f  u  B(u)du 

J—oo 

(2.17) 

=  uju7  =  u2 

(2.18a) 

=  r  ti2  B(u)du 

J  —  OO 

(2.186) 

r  •v'S? 

(2.19) 

,T  ri'dt  1  3_, 

—  =  —  /  t i3B(u)du 

)  a * 

(2.20) 

7  u*dt  1  f°°  4Df 

—  =  ^  j.J B(u>du 

(2.21) 

The  probability  density  function  B(u)  represents  the  number  of  u’s  found  over  an  interval 
ttj  and  Uj  +  A «<,  weighted  by  the  width  of  the  interval,  and  is  a  function  of  all  u,’s  tested.  The 
function  thu6  defined,  has  unit  area  under  the  curve  as  depicted  below. 

I  B(u)du  =  1  (2.22) 

J—OO 

As  shown,  o  is  a  measure  of  the  rms  value  of  the  fluctuation,  and  the  skewness  5,  gives  a 
representation  of  the  symmetry  of  u.  Positive  skewness  means  the  left  side  of  the  probability 
curve  is  steeper  than  the  right  side  and  vice  versa.  5  =  0,  is  a  symmetric  prob&blility  density  as 
represented  by  a  truly  Gaussian-like  random  function  with  zero  skewness.  K,  otherwise  known 
as  the  Kurtosis  or  flatness  factor,  is  a  measure  of  the  curvature  of  B(u)  (or  how  rounded  it  is) 
at  its  maximum  value.  We  expect  K  for  the  velocity  to  have  high  curvature  (very  sharp  profile) 
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dose  to  the  centerline  and  very  low  curvature  (less  sharp  or  very  rounded)  as  one  approaches 
the  wall. 

Only  a  complete  set  of  moments  of  all  orders,  contain  all  the  turbulent  information.  How¬ 
ever,  some  distributions  can  be  described  by  a  Unite  set  of  moments.  A  Gaussian  distribution 
is  completely  described  by  the  2nd  moment  a.  AH  odd  moments  are  equal  to  zero  and  all 
even  ^pments  n  >  2  are  equal  to  a  constant  times  the  2nd  moment.  For  instance,  5  =  0  and 
K  =  3.  The  departure  of  turbulence  from  Gaussian  in  nature,  is  a  manifestation  of  the  highly 
non-linear  behavior  of  turbulence.  Note,  the  2-parameter  probability  density  functions  at  a 
point  also  exist  as  shown  by, 

Ru,u,  =  uTuT  =  lim  i  [  u iUjdt  =  f  f  UiUjB(ui,Uj)duiduj  (2.23) 
T-oo  T  Jo  J-oc  J-eo 

Here,  Ru,u,  equals  the  normal  pressure  stresses  due  to  velocity  for  i  =  j  and  the  Reynolds  stress 
terms  when  i  ft  j  .  The  6ame  definitions  hold  true  for  the  pressure  P.  Schewe  (1983)  showed 
that  the  probability  density  function  B(p)  can  be  a  very  important  indicator  of  transducer 
cancellation  effects  related  to  the  size  of  the  transducer.  B(p)  shows  that  the  pressure  signature 
contains  amplitudes  as  high  as  15tr  which  have  very  short  lifetimes  and  contribute  significantly 
to  the  value  of  o  calculated,  depending  on  transducer  size. 

2.2  Statistical  Relationships  of  Turbulence 

As  illustrated  above,  the  derived  averages  are  all  first-order  time  mean  values  with  the 
end  result  equal  to  some  constant,  depending  on  the  parameter  being  measured  and  on  some 
independent  variable,  e.g.  time,  or  a  spatial  coordinate.  However,  when  conducting  a  detailed 
analysis  of  a  turbulent  flow,  we  are  interested  in  the  various  scales  of  turbulence,  their  energy 
content  and  convective  properties,  and  how  these  various  scales  interact  one  with  another.  It 
becomes  necessary  to  investigate  the  wave-like  behavior  using  second  order  statistical  functions 
obtained  with  fourier  transformations.  We  start  by  characterizing  a  turbulent  flow  in  6uch  a 
way  as  to  give  meaning  to  the  following  statistical  relationships. 

The  turbulence  quantities  of  interest  are  random  functions  of  both  space  and  time  and 
are  considered  to  be  stochastic  processes.  A  stochastic  process  refers  to  dealing  with  functions 
(such  as  pressure  or  velocity  components)  which  are  functions  of  an  independent  variable  or  a 
specified  combination  of  independent  variables  such  that, 

Y  -  Y  («1,S2,S3 . s„)  (2.24) 
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It  may  also  mean  that  a  series  of  functions  Y,  exist  with  independent  variables  s,  such  that  an 
ensemble  average  can  be  characterized  through  appropriate  statistical  properties.  This  means 
that  we  are  dealing  with  some  statistical  comparison  between  dependent  variables  such  as  Vj 
and  yj  which  are  functionally  dependent  on  the  same  prescribed  set  of  independent  variables  s, 
throughout  the  comparison  process.  Since  “statistics”  implies  a  time  average,  it’s  not  surpising 
that  one  of  the  independent  variables  of  most  interest,  is  the  variable  depicting  time  such  as  t 
or  r  =  t  +  At.  We  further  restrict  our  averaging  process  to  be  stationary  in  time.  This  means 
a  process  for  which  the  defining  statistics  remain  the  same  by  shifting  the  independent  variable 
by  At.  For  instance,  suppose  we  have  a  dependent  variable  Y  =  F(t)  for  which  we  collect  a 
series  of  snapshots  over  a  time  period  T  in  sequential  time  slots  and  then  form  an  ensemble 
average  by  summing  all  the  snapshots  as  follows, 


Yt 

—  I  (to,  #i,  1 2, .. 

...,  <0, ) 

y2 

-  Y(ta„ . 

••••? 

y3 

-  v(f.# . 

Fn-1  =  1  (ta._  . »to,_i) 

In  =  . o«) 

n 

V(t)  =  <  Yi  >  =  £  Yi  .  (2.25) 

i 

The  braces  <  ...  >  represent  the  ensemble  averaging  process.  As  noted,  this  particular  process 
began  at  time  t  =  to-  We  can  form  another  ensemble  average  begining  at  time  t  =  to  +  At 
such  that  yj(to)  =  y^to  +  At).  If  these  two  ensembles  are  equal  and  in  fact,  are  equal  to 
any  other  ensemble  created,  we  call  y(t)  a  stationary  process.  Finally,  if  each  record  of  the 
ensemble  is  statistically  equivalent  to  every  other  record  of  the  ensemble  such  that  we  can 
replace  an  ensemble  averaged  process  with  a  time  average  over  a  longer  time  interval  of  a  single 
representative  record,  this  process  is  considered  “ergodic”.  An  ergodic  process  is  one  in  which 
ensemble  averages  may  be  replaced  with  time  averages  and  vice-versa.  In  summary,  this  study- 
concerns  the  examination  of  a  fully  developed  channel  flow  in  which  the  defining  statistics 
concern  stochastic  processes  which  may  be  considered  as  stationary  and  ergodic. 
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The  most  general  statistical  function  of  a  variable  l'(x,t)  can  be  represented  as  the  space- 
time  correlation  covariance  as  represented  by, 

C«.l,i(x,x,,f,f')  =  <  tij(x,t)uj(x',t')  >  (2.26) 


where  the  braces  <  ...  >  again,  represent  either  a  time  average  or  a  suitable  ensemble  average 
according  to  the  above  definitions.  The  auto-correlation  is  defined  for  «  =  j  and  the  cross- 
correlation  for  t  j.  If  CUiUj  is  a  stationary  function  of  time,  and  homogeneous  in  space,  then 
Cu,u,  is  only  a  function  of  the  time  difference,  r  =  t'  -  t  and  the  separation  vector  {  =  x'  -  x 
such  that, 

Cu.ujiC, 7)  =  <  u<(x,t)ii>(x  +  ^,t  +  r)>  (2.27) 

To  give  meaning  to  the  averaging  process,  we  may  use  the  definitions  of  “stationarity”  and 
“ergodicitv”  to  define  Cu,Uj  by, 


T 

C„|U)(£,r)  =  lim  i  /  u i(x,t)uj(x  +  (,t  +  r)dt  (2.28) 

T— 00  ■»  Jo 

It’s  very  difficult  to  measure  the  full  space- time  correlation  covariance,  since  equation  (2.27)  or 
(2.28)  indicate  an  average  over  an  infinite  spatial  domain.  If  one  had  an  array  of  transducers 
separated  in  space,  one  may  average  each  one  in  time,  store  on  some  medium  such  as  a  hard 
disc  or  tape,  then  perform  the  spatial  integration.  However,  for  velocity  and  pressure,  the 
integration  is  three  and  two  dimensional  respectively,  a  process  not  altogether  easily  accom¬ 
plished.  What  is  usually  done  in  the  laboratory,  is  to  compute  the  time  covariance  only,  at  a 
fixed  spatial  separation  distance  {  =  s,  such  that, 


Cu.Uj-^ii  t) 


lim 

T~oo 


f  Jo  Ui^UJ^t  +  T^dt 


(2.29) 


where  the  cross-correlation  is  for  t  j  and  the  auto  correlation  i6  for  i  =  j.  Often,  it  is  more 
convenient  to  examine  the  space-time  correlation  functions  with  their  Fourier  transform  pair 
variables  (k,w)  by  the  following  transformation  of  CU,U/ . 


-  ////: 


Ctijtt,  (£,  r)exp[-m(k  •  {  +  ur)]drdi 


(2.30) 


m  =  y/—l 
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Equation  (2.30)  is  referred  to  as  the  full  auto/cross-spectra  density  function,  depending  on 
i  —  j  and  i  ^  j  respectively.  Finally,  the  Fourier  transform  pair  is  complete  by  inverse  Fourier 
transforming  equation  (2.30)  to  obtain  the  covariance  once  again. 

=  ^7  J  J  J  J  (k,w)eip[m(k  •  £  +  uT)}<L)dk  (2.31) 

One  of  the  parameters  of  most  interest  in  this  study  is  the  rms  pressure  of  the  fluctuations  at 
the  wall.  The  auto-covariance  of  the  pressure  may  be  defined  by, 

7* 

cw(T)  =  J1™  \  t  P(*M*  +  T)  dt  (2.32) 

T—oo  1  JQ 

and  if  r  =  0,  we  have, 

C„( 0)  =  Jim  i  /  p(t)2  dt  (2.33) 

00  1  Jo 

For  no  spatial  separation,  the  fourier  transform  of  the  covariance  of  pressure  as  defined  via 
equation  (2.30)  is, 

*vv(u)  =  /  Cpp(T)exp[-m{uT)]dT  (2.34) 

J  —  OO 

and  it’s  Fourier  transform  partner  is, 

cpp{t )  =  *PP(u)exp[+m{uT)}(Lj  (2.35) 

If  we  set  r  =  0,  we  obtain  the  additional  relation  for  rms  pressure  as, 

C„(r  =  0)  =  P}m,  =  (Sj <2'36> 

One  may  note,  that  there  are  three  methods  of  obtaining  the  rms  pressure  squared;  by  a  long 
time  average  of  the  square  of  the  pressure  signal,  evaluation  of  the  auto-covariance  function  at 
r  =  0,  or  by  integrating  the  auto-spectrum  $(w)  over  all  frequencies.  Note,  all  three  general 
functions  are  mathematically  related  however,  laboratory  hardware  usually  determines  which 
function  is  more  easily  obtained.  In  the  case  of  buried  noise  within  the  pressure  signal  itself, 
one  or  more  methods  may  be  eliminated  due  to  the  contamination.  Assuming  that  p(t)  contains 
noise,  a  simple  time  average  is  contaminated  and  thus  equation  (2.36)  may  be  the  only  way 
to  obtain  p*m#,  provided  an  accurate  measurement  of  $(u?)  is  obtained.  This  issue  will  be 
discussed  further  in  chapter  three. 

Of  particular  interest  is  the  cross-spectral  function  of  the  wall  pressure  9PiJ>.  between  two 
transducers  separated  over  a  spatial  distance  (  in  the  flow  direction.  The  streamwise  coherence 
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function,  72,  is  derived  by  normalizing  the  magnitude  squared  of  with  the  square  of  the 
auto  spectrum  as  demonstrated  below. 


Bull  (1967)  found  this  function  to  decay  exponentially  when  plotted  versus  the  cross  spectrum 
phase,  and  appears  to  universally  represent  similarity  scaling. 


For  a  propagating  (or  convective)  field,  the  phase  of  the  cross-spectrum  has  been  shown 
by  Blake  (1984)  and  Farabee  (1986)  to  be  proportional  to  the  convective  speed,  Uc(v,().  This 
relationship  can  be  expressed  by, 

«“■<)  -  -sj&5  (2'38) 

where  a(u>,£)  is  the  narrow  band  phase  which,  due  to  the  dispersive  nature  of  turbulence,  is  not 
necessarily  a  linear  function  of  frequency.  The  convection  velocity  Uc  for  the  pressure  producing 
eddies,  as  defined  by  equation  (2.38),  represents  the  convective  speed  of  all  wavenumber  scale 
components  contained  within  a  narrow  frequency  band  centered  at  the  frequency,  u>. 


Prior  to  the  last  decade,  a  single  convection  speed  for  the  turbulence  was  determined 
for  all  scales  and  frequencies  by  analyzing  the  results  from  the  cross-correlation  covariance  of 
two  hot  wire  sensors  separated  in  space  or  time.  Consider  equation  (2.28)  for  the  fluctuating 
component  of  velocity  in  the  flow  direction  over  space  and  time.  Following  the  previously 
established  conventions,  this  is  the  auto  correlation  (t  =  j)  over  one  spatial  coordinate  for  all 
time  and  can  be  represented  by, 

=  lim  ^  I  u(x,t)u(x  +  {,t  +  r)dt  (2.39) 

T—oo  i  Jo 

This  total  function  is  impossible  to  measure  instantaneously.  However,  suppose  the  following 
process  is  performed.  First,  various  values  for  time  delay  are  chosen,  e.g.  r  =  nr0  and  the 
distance  {  continuously  varied.  Secondly,  the  similar  but  opposite  operation  is  performed,  that 
of  chocwing  fixed  (  =  but  evaluated  for  all  time.  Therefore,  the  following  relations  for  the 
■pace  correlation  at  successively  fixed  tJs  and  time  correlation  for  succesively  fixed  {<>'«  are 
represented  respectively  as, 

{fixed  t0)  =  Cuu((,To)  =  lim  ^  I  u(x,i)u(x  +  (,t  +  rQ)dt  (2.40) 

T—  oo  i  Jo 
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(2.41) 


i  r 

{fixed  to)  =  Cutt(£0,T)  =  lim  -  /  u(x,t)t«(x  +  £0,<  +  r)dt 

7— oc  1  J0 

Equations  (2.40)  and  (2.41)  for  various  £,’s  and  r,’s  are  illustrated  in  figure  2.1.  As  shown,  both 
covariances  illustrate  the  convective  nature  of  turbulence.  A  convection  velocity  UC(  obtained 
from  figure  2.1a  by  the  value  of  £  =  £mar  for  the  maximum  in  each  hump  for  each  r0  is  defined 
as, 


dCuu(Z,nr0) 

di 


=  0 


Uc{  = 


(2.42) 


'  on 


Similarly,  a  related  convection  velocity  UCr  may  be  obtained  from  a  similar  analysis  of  the 
respective  humps  in  figure  2.1b  by, 


dCuu(n(0,r)  _ 
dr 


(2.43) 


Equation  (2.39)  is  represented  in  figure  2.2  and  may  be  obtained  from  either  equation 
(2.40)  or  (2.41)  by  summing  over  all  t0's  or  £0's  respectively.  As  indicated  in  figure  2.2,  the 
two  definitions  for  convection  velocity  are  not  the  same.  As  shown,  the  contours  as  represented 
in  this  figure  (and  in  real  turbulence)  are  elongated  and  somewhat  curved  in  the  r^-plane. 
Only  in  isotropic  turbulence  (invariant  scales)  within  which  Taylor’s  hypothesis  holds  (frozen 
pattern  with  no  temporal  rate  of  change  in  the  pattern),  are  these  two  definitions  equal. 

As  a  final  remark,  equation  (2.40),  and  hence  the  convection  velocity  UC(  as  defined  by 
equation  (2.42)  is  difficult  if  not  impossible  to  measure.  What  is  more  easily  obtained  and 
reported  in  the  literature  is  to  evaluate  equations  (2.41)  and  (2.43)  for  UCr  as  illustrated  in 
figure  2.1b.  Referring  to  the  curves  in  this  figure,  one  may  note  that  the  peak  for  r  =  0,  besides 
being  of  higher  intensity,  has  a  much  narrower  peak  and  sharper  roll  off  than  for  larger  values  of 
spatial  separation.  This  broadening  of  the  so-called  “hump”  with  increasing  spatial  separation 
6hows  that  the  correlation  is  due  to  larger  structures  extending  over  a  wider  region  of  the 
flow  field.  Thus  the  eddy  field  is  convected  downstream  with  the  larger  structure  remaining 
coherent,  while  the  higher  frequency  smaller  scale  structure  loses  its  identity. 


2.3  Macro/Micro  Length  and  Time  Scales 

As  mentioned  previously,  characterizing  the  scales  of  turbulence  over  both  space  and  time 
is  of  extreme  importance  in  furthering  our  knowledge  of  this  random  process.  The  correlation 


function  is  defined  by  dividing  the  covariance  function  by  the  mean  square  value  o1 .  Taking 
the  time  correlation  of  velocity  as  defined  by  equation  (2.41)  for  zero  separation  results  in, 


*..<T)  =  (2.44) 

which  is  illustrated  in  figure  2.3.  This  function  represents  a  measure  over  time  of  how  long  the 
velocity  is  correlated  with  itself  for  all  structures  within  the  flow.  A  macro  time  scale,  Tmai< 
can  be  defined  by  integrating  equation  (2.44). 


(2.45) 


Tmax  is  a  measure  of  the  maximum  amount  of  time  over  which  the  eddie-  remember  their  path. 
We  may  also  define  a  minimum  or  micro  time  scale,  rm,„,  by  fitting  a  parabola  to  the  function 
in  equation  (2.44)  starting  at  time  r  =  0  and  evaluating  the  2nd  derivative. 


(0 

Ruu{r)  = 

'min 

(«) 

f  SR' 

-2 

Ur2. 

T  =  0 

T2 

min 

Similar  relations  exist  for  macro  and  micro  spatial  scales,  as 
defined  by, 


(2.46a) 

(2.466) 

Lx  and  lx  respectively  and  are 

(2.47) 

(2.48) 


The  macro-spatial  scale  Lx  is  a  relative  measure  of  the  size  of  the  eddies  containing  most  of 
the  energy  and/or  the  longest  distance  over  which  eddies  of  a  given  size  are  correlated,  and 
lx  is  an  approximate  measure  of  the  minimum  scale  at  which  dissipation  occurs.  Kraichnan 
(1956)  argued  that  integration  of  the  pressure  correlation  function  Rvv(t)  over  a  flat  plate  (and 
channel)  should  be  zero,  hence  for  pressure,  no  macro  time  scale  can  be  computed.  Hansen, 
et  al. (1987)  showed  that  this  might  not  be  true,  that  channel  flows  appeared  to  be  able  to 
support  pressure  scales  of  infinite  length  at  fixed  frequency.  However,  neither  hypothesis  gives 
an  alternative  definition  via  pressure  statistical  functions.  One  can  only  define  a  qseudo-macro 
scale  T0,  defined  as  that  value  of  time  r  where  the  time  correlation  function  Rpv{t)  =  0  at  its 
first  zero  crossing.  From  Taylor’s  hypothesis,  one  can  define  a  macro  length  scale  as, 


Lx  =  U0Te 


(2.49) 
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where  U0  is  either  the  convection  velocity  or  maximum  velocity  in  the  flow.  For  the  micro 
scales  of  length  and  time,  a  similar  analysis  as  outlined  in  equation  (2.46)  and  figure  2.3 
may  be  attempted  for  the  pressure  correlation.  However,  in  practice,  the  effort  expended  is 
somewhat  subjective  and  requires  creativity  on  the  part  of  the  researcher,  because  there  are 
usually  only  a  few  points  defining  Rpp(r)  close  to  zero.  Therefore,  using  an  analysis  first  put 
forward  by  Taylor  (1938),  a  relation  for  the  time  micro-scale  rmtn  can  be  defined  as, 

4-  =  K  /”/’*(/)<«■  (2-50) 

and  using  equation  (2.49),  a  micro-length  scale  lx  can  be  defined  as, 

lx  —  U0  Tmtn  .  (2.51) 

These  macro  and  micro  scales  are  computed  and  further  discussed  in  chapter  5. 


2.4  Relevant  Measured  Time-Dependent  and  Statistical  Quantities 

In  order  to  characterize  the  kinematics  of  the  flow  for  the  physical  experiment,  the  mean 
value  of  velocity  in  the  flow  direction  U(x,z)  and  the  fluctuating  component  ti(x,z,i)  were 
computed  from  the  instanteous  data  6et  U{x,y,z,t).  This  data  set  was  obtained  via  laser 
Doppler  velocimetry  and  is  explained  in  detail  in  the  next  chapter.  The  longitudinal  pressure 
gradient  dP0(x)/dx ,  along  the  channel  centerline,  was  measured  directly  during  each  Reynolds 
number  experimental  run  from  which  the  wall  shear  stress  and  shear  velocity  could  be  calculated 
via  equation  (2.14). 

A  five  element  array  of  flush  mounted  pressure  transducers  were  recorded  simultaneously  in 
time,  at  various  flow  velocities.  These  transducers  only  respond  to  the  fluctuating  component  as 
represented  by  p(x,  z,  t),  and  will  be  discussed  in  more  detail  in  the  next  chapter.  A  number  of 
first-order  and  second-order  multi-dimensional  time/statistical  functions  have  been  discussed. 
Intuitively,  one  may  hope  that  most  relations  presented  can  feasibly  be  measured.  In  a  realistic 
sense,  all  that  is  measured  is  the  time  dependent  pressure  at  the  wall  and  the  time  dependent 
velocity  at  various  points  within  the  flow  field.  These  data  are  either  processed  immediately 
with  dual  channel  analyzers,  or  digitized,  stored  on  some  suitable  medium  and  processed  later 
in  an  appropriate  manner.  In  this  experimental  study,  the  wall  pressure  p(x,z,t)  was  recorded 
simultaneously  for  five  channels  of  data  for  each  Reynolds  number.  This  data  was  stored  on 
digital  tape,  and  processed  later  by  software  developed  by  the  author.  The  highest  order 
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statistical  function  obtained  in  this  study  is  the  second  order  2-dimensional  auto/cross  spectral 
density  function  $pp(£r»/)  of  the  wall  pressure  fluctuations.  All  other  2nd  order  statistical 
functions,  such  as  the  space- time  covariance,  as  defined  by  equation  (2.41),  are  derived  from 
this  function  via  Fourier  transformation,  and  all  first  order  relations,  such  as  the  probability 
density  function,  are  obtained  directly  from  p(z,z,i). 

Below  is  a  summary  of  all  values,  either  measured  directly  or  computed,  that  are  used  in 
this  study.  The  details  of  measuring  or  computing  these  parameters  will  be  given  in  the  next 
chapter. 


U(x,  y,  z,t) 
u(z,y,z,t) 
U(z,y,z) 

®rm» 

&Po(x,t) 

*pp(f) 

Prms 

°p.pi  (£♦  /) 

VAU) 

Lx  »  It 

To  , 


streamwise  velocity 
streamwise  velocity  fluctuation 
streamwise  mean  velocity 

root-mean-square  of  streamwise  velocity  fluctuations 

static  pressure  gradient 

pressure  frequency  spectrum 

root-mean-square  of  the  pressure  fluctuations 

pressure  cross-spectrum 

cross-spectrum  phase 

narrow-band  convection  velocity 

Macro,  micro  length  scales 

Macro,  micro  time  scales 


2.5  Statistical  Error  Definitions 

The  full  space-time  power  spectral  density  functions  may  be  computed  by  Fourier  trans¬ 
forming  correlation  covariance  functions  through  equation  (2.30)  or  by  direct  Fourier  transfor¬ 
mation  of  the  time  history  records.  Since  the  introduction  of  finite  Fourier  transform  algorithms 
around  1965,  the  latter  approach  has  dominated  due  to  the  faster  processing  times.  For  in¬ 
stance,  the  time  to  compute  the  correlation  covariance  through  equation  (2.29),  then  Fourier 
transforming  into  spectral  space,  typically  takes  approximately  N2  calculations  where  N  is 
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the  number  of  ensemble  averages.  The  advent  of  the  Cooley-Tukey  algorithm  in  1965  for  the 
calculation  of  a  complex  Fourier  series  requires  approximately  Nln(N)  calculations.  Special 
firmware  has  been  designed  with  programmed  instructions  for  fast-Fourier-transformation  and 
are  the  primary  instructional  set  within  dual  channel  analyzers  used  in  the  laboratory. 

However,  these  computations  are  only  estimates  and  not  the  actual  or  true  values  of 
the  complete  functions,  as  one  always  deals  with  finite  time  records.  Estimates  of  the  errors 
in  computing  such  functions  have  been  derived  by  Bendat  and  Piersol  (1978,1980).  These 
errors  form  two  classes  described  as  bias  and  random.  Bias  errors  are  systematic  errors  that 
appear  with  the  same  magnitude  and  at  the  same  frequencies  from  one  analysis  to  another. 
For  instance  background  pump  noise  or  vibrational  energy  transferred  through  walls  to  the 
transducer  contribute  to  bias  errors.  Random  errors,  on  the  other  hand,  are  considered  a 
haphazard  scatter  in  the  results  and  are  due  to  averaging  operations  conducted  over  a  finite 
time  window'  or  number  of  ensembles.  Bendat  and  Piersol  have  shown  biasing  errors  to  be 
second  order,  thereby  contributing  less  overall  decay  to  spectral  estimates.  For  the  broad  band 
spectrums  (typical  of  velocity  and  wall  pressure  spectral  functions),  random  errors  have  more 
significance.  The  non-dimensional  error  estimates  c,  for  the  spectral  functions  computed  in  this 
experimental  study,  are  listed  below'  with  jV  equal  to  the  total  number  of  ensembles  or  records 
and  Np  equal  to  the  total  number  of  points  per  record. 


auto  —  spectrum . e[$PJ>(u;)]  —  -j= 

vAT 

(2.52) 

cross  —  3pecfrtim....f[$pp({,/)]  =  -  ■= 

7  iZJWN 

(2.53) 

. ■‘b’K./M  = 

(2.54) 

*“*■ . = 

(2.55) 

fluctuations . <fp^J  = 

(2.56) 

All  of  the  quantities  listed  on  the  righthand  side  of  the  above  equations  are  considered  to  be 
the  true  functions.  Their  values  are  not  known  a  priori,  so  the  common  practice  is  to  insert  the 
calculated  values  for  these  functions.  As  mentioned,  the  finite  time  records  available  in  most 
analyses  are  the  primary  contributor  to  the  random  error,  hence  all  the  error  definitions  above 
are  shown  to  be  proportional  to  N ,  the  number  of  records  or  ensembles.  While  the  error  in  the 
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auto  spectrum  is  proportional  only  to  N  and  that  for  the  rms  of  the  fluctuations  proportional 
to  JVjVp,  all  the  other  error  estimates  are  functions  of  the  true  coherence  between  two  signals 
and  are  also  functions  of  frequency.  In  either  case,  to  decrease  the  error,  hence  improving  the 
confidence  level  in  the  accuracy  of  obtaining  all  the  functions,  it  is  apparent  from  equations 
(2.52)  to  (2.56)  that  the  larger  the  value  of  A\  the  more  accurate  are  the  computed  estimates. 
In  the  present  study,  the  computational  and  physical  experiments  had  values  of  N  equal  to 
2048  and  1750  respectively.  This  results  in  error  estimates  of  0.024  and  0.022  for  the  power 
spectral  density  function  for  the  physical  channel  flow  and  numerical  simulation  experiments 
respectively. 
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FIGURE  2.3.  AUTO  TIME-CORRELATION  COEFFICIENT  FOR 
VELOCITY  (zero  spatial  separation) 


CHAPTER  3 


DESCRIPTION  OF  PHYSICAL  AND  COMPUTATIONS 
EXPERIMENTS 


3.1  Physical  Experiment 

Upon  the  conception  of  this  experimental  endeavor,  a  number  of  critical  elements  had  to 
be  considered.  The  first  was  to  choose  a  facility  that  would  ensure  the  presence  of  a  fully 
developed  two-dimensional  turbulent  flow  with  low  centerline  turbulence  intensity  level,  such 
that  boundary  layer  variables  (both  dynamic  as  well  as  kinematic)  could  accurately  be  mea¬ 
sured.  Secondly,  due  to  the  small  intensity  levels  of  the  wall  pressure  being  made,  this  facility 
should  have  low  noise  characteristics,  in  terms  of  facility  related  vibration  and  propagating 
acoustic  noise  which  typically  contaminate  wall  pressure  measurements  in  liquid  flow  envi¬ 
ronments.  Thirdly,  the  facility  should  accomodate  test  plates  of  various  configurations  and 
retain  an  uncontaminated  flow  characterization,  not  withstanding  numerous  insertions  of  such 
devices.  Finally,  the  facility  should  have  a  wide  range  of  operating  conditions. 

While  the  aforementioned  specifications  may  appear  formidable,  they  represent  the  basis 
for  performing  high  quality  measurements,  particularly  of  a  generic  nature,  to  facilitate  accu¬ 
rate  comparisons  with  similar  experiments.  In  particular,  it  was  mentioned  previously  that 
these  experimental  results  would  serve  to  validate  the  correctness  of  more  recent  numerical 
simulations  which  solve  the  full  Navier-Stokes  equations  in  the  same  geometrical  configuration. 
We  can  not  therefore,  overemphasize  the  importance  of  ensuring  that  the  above  designated 
criteria  be  strictly  enforced. 


3.1.1  Channel  Description 

The  chosen  facility  was  originally  constructed  in  1978,  at  the  Naval  Research  Laboratory, 
under  the  auspices  of  the  Office  of  Naval  Research  (ONR)  compliant  coating  research  program 
of  study.  A  unique  feature  characteristic  of  this  facility  is  that  it  operates  in  a  blow-down 
configuration  to  eliminate  contaminating  noise  due  to  structure  vibration  and/or  vibrationally 
induced  acoustic  phenomena  which  are  typical  of  a  recirculating  system.  The  flow  through  the 
channel  can  be  driven  by  a  gravity  head  alone  or  by  a  combination  of  increasing  air  pressure 
and  gravity  head  which  serve  to  produce  a  constant  discharge  rate  or  driving  pressure. 
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Figure  3.1  shows  a  schematic  of  the  NRL  water  channel  in  its  present  configuration.  The 
major  features  of  the  facility  are  a  pressure  tank,  upstream  liquid  reservoir,  turbulence  man¬ 
agement  section,  entrance  section  for  the  developing  turbulent  flow,  the  test  section  and  a 
downstream  reservoir.  All  of  the  features,  other  than  the  pressure  tank  and  upstream  reser¬ 
voir,  are  enclosed  in  a  temperature  controlled,  acoustically  isolated  laboratory.  The  pressure 
tank  has  a  volume  of  14.16m3  with  a  working  pressure  of  552kpa  (SQpsig).  It  is  equipped  with 
an  inlet,  outlet,  relief  valve,  drain,  pressure  gage  and  an  inspection  manhole.  The  pressure 
tank  is  continuously  charged  to  its  working  pressure  with  flow  control  valves  at  the  outlet 
for  regulating  the  discharge  pressure.  The  upstream  liquid  reservoir  is  a  stainless  steel  tank 
approximately  10m  high,  with  a  diameter  of  3m  and  a  volume  of  28.4m3  (7000paJ).  There 
are  five  ports  of  differing  sizes  to  accomodate  a  primary  inlet  and  outlet  and  for  various  other 
instrumentation  purposes.  The  tank  is  filled  through  a  5.08cm  diameter  pipe  leading  to  the 
bottom  of  the  reservoir.  The  outlet,  which  couples  at  the  same  location  as  the  fill  pipe,  is 
a  10.15cm  copper  pipe  which  is  routed  inside  the  laboratory  through  a  flow  control  valve  to 
the  turbulence  management  section.  At  the  downstream  end  of  the  facility  is  a  reservoir  for 
collection  of  the  working  fluid.  This  is  an  open  tank  of  welded  steel  construction  with  a  total 
capacity  of  28.4m3  (7000 gal). 

In  figure  3.2,  a  more  detailed  schematic  indicates  the  interior  components  of  the  channel. 
The  turbulence  management  section  consists  of  a  number  of  sectional  screens  and  honeycombs 
with  an  entrance  component  having  a  contraction  ratio  of  approximately  7  to  1.  The  de¬ 
veloping  turbulent  boundary  layer  proceeds  through  a  constant  cross-sectional  area  member, 
approximately  4  meters  in  length,  which  has  an  aspect  ratio  of  18  to  1  with  a  channel  height 
of  2h  =  2.54 cm  (1  inch).  The  test  section  contains  a  number  of  large  plates  within  which 
smaller  test  windows  or  transducer  plates  may  be  arranged  in  many  different  configurations. 
These  windows/plates  allow  for  the  placement  of  small  pressure  probes  or  the  non-intrusive 
access  of  the  beams  of  a  laser  Doppler  anemometer.  The  ratio  of  downstream  distance  of  the 
measurement  point  to  the  half-height  of  the  channel  x/h ,  is  320.  This  value  of  x/h  is  well 
within  the  region  of  fully  developed  turbulent  channel  flow  (Hussain  and  Reynolds  1975). 

The  side  walls  of  the  channel  are  constructed  from  2.54cm  high  by  10.16cm  (1*4  inches ) 
wide  continuous  pieces  of  clear  plexiglass.  For  single  point  velocity  measurements,  an  LDV 
system  may  enjoy  access  through  these  side  panels  if  necessary.  The  entire  channel  system  is 
located  approximately  1.22m  (Aft)  off  the  floor  on  massive  aluminum  I-beams.  The  massive 
supports  ensured  that  a  minimal  amount  of  vibrational  energy  would  be  transferred  through 
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to  the  channel  inner  walls  and  subsequently,  the  working  fluid.  The  test  section  is  described  in 
more  detail  later  in  this  chapter. 

The  above  components  coupled  with  a  full  reservoir  and  the  charged  pressure  tank,  provide 
a  range  of  velocities  at  the  test  section  of  0.22  to  2.44  m/sec  (0.72  to  6.0  ft/ sec)  corresponding 
to  Reynolds  numbers  based  on  channel  half  width  (Rh.  —  Uavgh/ v)  of  3000  to  approximately 
30,000. 


3.1.2  Operational  Characteristics 

The  basic  principle  of  operation  is  to  control  the  bulk  flow  rate  in  the  test  section  by 
maintaining  a  constant  total  dynamic  head  or  driving  pressure  from  the  inlet  to  outlet.  This 
has  typically  been  characterized  as  a  “blow-down”  facility.  By  monitoring  the  pressure  at  some 
point  in  the  upstream  reservoir  and  subtracting  the  pressure  at  the  channel  exit,  a  differential 
(dynamic)  driving  pressure,  characteristic  of  a  specific  bulk  or  average  velocity  can  be  contin¬ 
uously  monitored.  As  mentioned  previously,  there  is  a  large  upstream  liquid  reservoir  which 
provides  the  working  fluid.  However,  the  quantity  of  fluid  is  finite  and,  as  the  flow  proceeds 
through  the  system,  this  measured  pressure  decreases  as  the  fluid  level  drops  in  the  reservoir. 
In  order  to  maintain  a  constant  driving  pressure,  air  is  added  above  the  liquid  at  such  a  rate 
as  to  compensate  for  the  liquid  deficiency.  From  the  equation  of  continuity  (for  incompressible 
flow),  we  may  define  the  following  relations. 


flow  rate  =  Q[t)  -  ArVr(t)  =  AeVe(t)  =  Constant  (3.1) 


Vr(t)  =  - 


dim 

dt 


(3.2) 


Here,  Q  is  the  discharge  rate,  Ar  and  Vr(t)  are  the  cross-sectional  area  and  average  velocity  in 
the  upstream  reservoir  and  Ac,  Ve(t)  are  the  corresponding  values  in  the  channel.  H(t)  is  the 
height  of  the  liquid  as  a  function  of  time.  Combining  equations  (3.1)  and  (3.2),  results  in, 


(3.3) 


By  monitoring  the  height  of  the  liquid  fluid  as  a  function  of  time  (hence  the  reservoir’s  velocity), 
equation  (3.3)  shows  a  simple  relation  for  determining  the  channel’s  velocity.  For  a  fixed  starting 
level  of  liquid,  one  would  obtain  the  same  velocity  for  all  runs.  Therefore,  as  shown  in  figure 
3.1,  a  flow  control  valve  is  inserted  in  front  of  the  turbulence  management  section.  By  choosing 
various  openings  of  the  valve,  one  may  obtain  a  number  of  velocities  as  indicated  previously. 
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A  Celesco  model  455  differential  pressure  transducer  is  used  to  effect  this  monitoring 
function.  This  transducer  requires  a  specific  stainless  steel  diaphragm  for  each  pressure  range 
of  operation  and  can  be  statically  calibrated  for  each  run  if  necessary.  For  the  blowdown 
procedure,  a  1  Ops  ip  diaphram  is  inserted  within  the  transducer  and  subsequently  calibrated 
as  the  reservoir  is  being  charged  with  working  fluid  (A  calibration  curve  for  this  transducer  is 
shown  in  the  next  section.).  The  output  is  0  to  10  voltsDC  with  a  frequency  response  good 
to  4000#*.  In  order  to  facilitate  the  monitoring  of  this  pressure  and  subsequent  calculation 
of  average  velocities,  an  HP-200  series  lab  scientific  computer  with  a  multi-channel  721  multi¬ 
processor  is  utilized.  The  Celesco  transducer  ouput  is  sampled  in  time  at  a  1000/7 z  with  500 
pts  taken  over  half-second  periods  and  averaged  such  that  at  time  tj,dH/dt  may  be  calculated. 

Figure  3.3  shows  typical  results  of  the  monitoring  instrumentation  for  four  different  con¬ 
stant  velocities  (Reynolds  numbers).  The  height  of  the  liquid  reservoir,  monitored  over  time, 
is  recorded  with  the  respective  lines  for  each  velocity  representing  a  least  squares  best  fit  to 
the  data.  As  the  equation  of  continuity  suggests,  the  data  is  linear  to  a  great  deal  of  accuracy; 
better  than  1%  repeatability  for  the  same  velocity.  With  respect  to  the  velocity  range  indicated 
previously,  the  total  amount  of  fluid  available  resulted  in  run  times  from  15mtn  to  over  an  hour 
for  the  highest  to  lowest  velocities  respectively. 


3.1.3  Test  Section 

The  channel  contains  large  removable  sections  that  provide  access  within  the  channel  as 
well  as  a  platform  for  the  insertion  of  smaller  test  plates.  Thus  the  experimenter  has  the 
ability  to  design  arbitrarily  any  desirable  test  plate/probe  combination  to  insert  within  the 
flow  system.  For  this  experiment,  a  7.97cm  (ID)  hole  was  drilled  within  the  center  of  one  of 
the  removable  sections  with  8  symmetric  bolt  holes  for  fastening  purposes.  A  small  test  plate 
7.97cm  (OD)  in  diameter  was  designed  for  insertion  of  the  small  flush-mounted  transducers. 
Figure  3.4  indicates  the  location  and  orientation  of  this  tefct  plate  with  respect  to  the  large 
removable  section,  and  figure  3.5  presents  a  detailed  illustration  of  the  placement  of  holes  for 
the  transducers.  As  shown,  there  were  20  holes  drilled  0.17  cm  in  diameter  in  a  linear  array 
approximately  0.20cm  apart.  This  linear  array  was  offset  from  the  centerline  by  2.54cm  with 
a  similar  hole  aligned  with  the  first  hole  in  the  array,  symmetrically  placed  the  other  side  of 
the  centerline.  This  spanwise  transducer  would  be  used  for  noise  cancellation  purposes.  Only 
5  transducers  were  available  for  this  experimental  study,  with  four  inserted  within  the  linear 
array. 
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The  eight  symmetric  bolt  holes  provide  a  means  to  rotate  the  plate  to  differing  angles  of 
theta  (with  respect  to  the  streamwise  direction)  to  extract  more  detailed  wavenumber  depen¬ 
dence  if  desired.  For  this  experiment,  data  was  taken  at  the  theta  angle  of  0  degrees  only  and 
will  be  characterized  as  inline  or  streamwise  data. 

Since  only  four  transducers  were  placed  in  the  flow  direction,  the  additional  holes  were 
filled  with  rtv,  smoothed  even  with  the  surface  and  lightly  sanded  with  emery  cloth  to  exact 
a  flush  surface.  The  entire  plate  was  anodized  to  prevent  rust  and/or  additional  corrosive 
deposits  from  adhering  to  the  surface. 

As  indicated  in  figure  3.4,  the  large  removable  plate  was  4.445cm  thick  and  also  indicates 
that  the  test  plate  was  of  a  similar  thickness.  In  order  to  ensure  a  flush  mating  of  both  plates, 
the  test  plate  after  fabrication,  was  inserted  within  the  large  removable  section  and  sanded 
down  to  a  tolerance  of  24.5 pm  (microns).  An  O-ring  maintained  a  water-tight  seal  once  the 
two  plates  were  inserted  within  the  flow  channel. 


3.1.4  Wall  Shear  Stress  Measurements 

The  geometry  being  investigated  with  rectangular  cross-section,  is  known  as  plane  Poisuille 
flow  when  the  aspect  ratio  is  large,  and  the  flow  is  steady  and  incompressible.  If  one  chooses 
a  differential  volume  located  at  that  point  where  the  flow  is  considered  fully  developed  and 
analyzes  the  applicable  surface  forces,  one  may  obtain  the  well  known  relation  for  wall  shear 
stress  as  represented  by, 

It  =  (3-5) 

Here,  h  is  the  channel  half-width,  P0  is  the  static  pressure  and  i  is  the  spatial  coordinate  in 
the  flow  direction.  Thus  by  an  appropriate  measurement  of  the  pressure  gradient  within  the 
channel  at  each  respective  Reynolds  number,  the  wall  shear  stress  and  hence  the  shear  velocity, 
u'  =  y/rwau/p  may  be  determined.  A  number  of  1/4  inch  taps  are  placed  along  the  side  wall 
of  the  channel,  approximately  0.3048m  apart.  Each  tap  was  connected  in  parallel  to  a  Celesco 
model  451  differential  transducer  containing  a  0.1  psig  stainless  6teel  diaphragm.  This  type 
transducer  is  simple  to  calibrate.  A  manometer  was  stragecially  placed  at  the  channel  test 
section  and  coupled  to  the  pressure  taps  via  appropriate  tubing.  Distilled  water  was  added 
to  the  manometer  at  various  heights  above  the  channel  centerline  which  provided  a  simple 
means  for  calibration.  Figure  3.6  shows  a  typical  calibration  curve  for  this  transducer  as  well 
as  the  previously  mentioned  model  used  for  monitoring  the  height  of  water  in  the  upstream 
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reservoir.  A  least  squares  curve  was  fitted  to  the  data  in  each  case,  and  as  shown,  the  linear 
characteristics  of  the  transducers  were  validated. 


During  a  particular  channel  run  at  constant  velocity,  the  pressure  gradient  transducer 
was  systematically  switched  through  each  of  the  12  pressure  taps  and  a  time  history  of  the 
transducer  output  was  recorded.  An  average  for  each  record  was  then  obtained  and  plotted 
versus  channel  location.  Figure  3.7  indicates  the  linearity  throughout  all  velocities  tested  and 
further  substantiates  that  the  channel’s  velocity  is  being  maintained  constant.  Since  the  static 
pressure  is  linear,  the  slope  of  each  line  (a  constant)  may  be  used  in  equation  (3.5)  to  determine 
the  wall  shear  stress  and  hence  shear  velocity,  for  each  Reynolds  number  investigated.  Jones 
(1976)  performed  an  extensive  study  demonstrating  that  high  aspect  ratio  plays  an  extremely 
important  role  in  determining  the  friction  factor  ‘A’.  The  friction  factor  is  defined  by, 

A  =  (3.6) 

Q 

Jones  showed  that  in  order  for  data  from  differing  aspect  ratio  ducts  to  demonstrate  similarity 
for  turbulent  flows,  the  Reynolds  number  based  on  hydraulic  diameter  must  be  modified.  Using 
the  following  definitions, 

D,  =  Dhydraulie  =  ,3.7) 

=  E^£i  (3.8, 

the  friction  factor  A,  can  be  plotted  vs.  Rd  wnich  is  defined  by  Jones  as, 

R‘dh  =  *(At)  RDh.  (3.9) 


V(A,)  represents  a  function  of  aspect  ratio  A„  which  Jones  was  able  to  determine  for  various 
sized  rectangular  ducts.  Jones  further  showed  that  to  second  order,  thi6  function  could  be 
represented  by  the  following. 


(3.10) 


In  figure  3.8,  typical  data  for  the  Reynolds  numbers  investigated  have  been  plotted  together 
with  the  well  known  relations  of  Prandtl  and  Blassius. 
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3.1.5  Measurement  and  Data  Acquisition  of  the  Mean  Flow  and  Turbulence 

Before  making  any  pressure  measurements  of  a  generic  nature,  a  characterization  of  the 
channel’s  velocity  was  desirable.  The  purpose  being  to  ensure  the  two-dimensionality  of  the 
flow  and  that  no  extraneous  large  scale  structure  from  the  side  walls  interferred  with  the  core 
flow. 

Velocity  was  measured  with  an  argon  TSI  laser  Doppler  velocimeter  (LDV)  model  No. 
9100-10  which  was  a  3-component  (5  beam  )  system  operating  in  backscatter  mode.  Orig¬ 
inally,  it  was  desired  to  pass  the  beams  through  a  window  located  in  one  of  the  removable 
plates  opposite  the  pressure  transducer  test  plate.  However,  when  the  first  measurements  were 
attempted,  it  was  found  that  the  closest  position  that  the  probe  volume  could  be  placed  to  the 
wall  was  on  the  order  of  0.025cm  which  was  half  the  length  of  the  ellipsoid  characterizing  the 
shape  of  the  probe  volume.  This  corresponded  to  a  Y+  value  of  approximately  100,  which  was 
in  the  outer  region  of  the  boundary  layer  completely  ignoring  the  viscous  sub-layer  and  buffer 
layer.  The  optics  which  were  available  prevented  reducing  this  characteristic  dimension  unless 
first  passing  through  a  beam  expander.  Doing  so  however,  resulted  in  a  reduction  in  the  focal 
length,  thus  preventing  going  through  the  window.  It  was  decided  that  the  clear  glass  making 
up  the  side  walls  of  the  channel,  provided  suitable  access  for  two  beams  lying  in  a  horizontal 
plane.  By  passing  the  beams  through  the  3.75  beam  expander  module,  one  could  reduce  the 
ellipsoid  diameter  to  approximately  0.015mm,  resulting  in  a  V+  value  on  the  order  of  5  to  15 
depending  on  the  velocity.  This  did  not  allow  for  a  complete  characterization  of  the  viscous 
sub-layer,  but  the  buffer  layer  was  well  indicated. 

An  optical  table  lm  by  2m  provided  a  flat  bed  upon  which  the  entire  optical  train  as  well 
as  the  laser  itself  was  placed.  The  optical  table  rested  on  the  two  dimensional  segment  of  a 
nulling  machine  which  provided  accurate  measurements  to  within  0.025mm  in  two  directions  of 
a  horizontal  plane.  The  entire  apparatus  then  rested  on  an  aluminum  frame  connected  at  three 
points  with  stepping  motors.  This  provided  a  simple  method  of  aligning  the  two  beams  in  a 
horizontal  plane  plus  providing  profile  capability  in  the  vertical  direction  to  within  0.015mm. 

The  two  beams  which  passed  through  the  side  walls  were  taken  from  the  blue  beam  of  an 
argon  laser  having  a  wavelength  of  488nm.  Both  beams  lay  in  a  horizontal  plane  and  one  was 
shifted  with  a  bragg  cell  by  2  M Hz  such  that  the  fringes  would  appear  optically  to  move  in  the 
opposite  direction  to  the  mean  velocity.  After  exiting  from  the  laser,  both  beams  pass  through 
the  optical  train  approximately  22 mm  apart  then  pass  through  the  beam  expander  module. 
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At  the  exit  from  the  final  lens,  both  beams  are  80mm  apart,  have  an  included  angle  of  1 8deg 
with  a  final  focal  lenth  of  455mm.  At  the  intersecting  point  of  the  two  beams,  an  ellipsoid  is 
formed  with  a  characteristic  length  and  diameter  of  0.15mm  and  0.03mm  respectively.  This 
resulted  in  being  able  to  traverse  to  the  wall  as  close  as  0.015mm  in  the  vertical  direction. 

In  backscatter  mode,  the  reflected  light  from  the  probe  volume  was  collected  by  the  return 
optics,  focused  on  a  photo-multiplier  tube  with  an  appropriate  analog  signal  sent  to  a  model 
9898  TSI  counter.  This  counter  included  all  appropriate  instrumentation  for  post  processing  of 
the  frequency  shifted  beams  and  scattered  light  to  include  validation  of  the  particle,  counting 
of  fringe  crossings  and  tabulation  of  the  time  between  validated  bursts.  This  information  was 
then  sent  in  real  time  to  an  HP-600  computer  and  stored  for  further  analysis  at  a  later  time. 
At  each  station  chosen  in  a  velocity  profile  for  each  respective  run,  12,800  validated  points 
were  collected,  stored  on  hard  disc  and  the  probe  volume  traversed  to  the  next  profile  station. 
These  data  could  then  be  processed  for  mean  and  fluctuating  components  of  velocity,  V  and 
ti  respectively.  Figure  3.9  shows  a  schematic  of  the  test  section  with  the  location  of  the  two 
coincident  beams.  Three  stations  were  chosen  to  investigate  and  are  indicated  by  the  points 
A,B  and  C  in  the  diagram.  As  mentioned  earlier,  the  channel  is  18  inches  wide,  and  these 
station  points  were  chosen  to  bracket  the  inner  core  of  flow  as  far  as  the  beam  focal  length 
would  reach.  Station  A,B  and  C  correspond  to  5,  8  and  11  inches  from  the  inside  of  the  near 
wall  and  bracket  the  area  of  the  small  wall  pressure  test  plate.  Three  profiles  were  taken  at  each 

location  to  check  repeatability.  Figure  3.9  also  indicates  a  traverse  labeled  “a - a" 

which  corresponds  to  the  path  made  by  the  focal  point  of  the  two  beams,  from  the  near  side  wall 
into  the  flow  at  various  distances  from  the  bottom  or  top  wall.  These  traverses  were  inserted 
to  substantiate  that  flow  was  indeed  two  dimensional  and  that  no  side  wall  effects  interfered 
with  the  core  flow.  For  more  detailed  information  on  the  general  theory  of  LDV  operation  and 
that  more  specific  to  the  model  9100-10,  the  reader  is  referred  to  Durst  et  al. (1981)  and  TSI 
(1985). 

3.1.6  Measurement  and  Data  Acquisition  of  the  Wall  Pressure  Fluctuations 

Five  small  transducers  were  embedded  flush  at  the  wail  within  a  small  test  plate  in  the 
test  section.  These  transducers  are  model  8514-10  piezo-electric  sensors  manufactored  by  EN- 
DEVCO.  Their  nominal  sensitivities  as  reported  by  the  manufacturer,  were  -227 dB  referenced 
to  with  a  flat  frequency  response  out  to  140 kHz.  The  active  area  of  the  transducer,  as 
reported  by  Galib  and  Zandia  (1984),  is  approximately  0.5mm.  This  dimension  gives  a  ratio  of 


h/dt  equal  to  25/1,  and  a  viscous  scale,  d+  =  (dtu’)f v  in  the  range  20  to  40.  The  signals  from 
each  transducer  were  amplified  to  a  suitable  level,  low-pass  filtered  to  prevent  aliasing,  and 
then  simultaneously  sampled  in  time,  digitized  and  stored  on  tape  for  post  processing.  Hence, 
for  each  velocity  run,  simultaneous  records  in  time  of  all  five  transducers  were  recorded. 

The  data/acquisition  was  accomplished  with  a  Quad  Systems,  Inc.(QSI)  16  channel  model 
721  data  processor.  This  unit  includes  a  12-bit  A/D  converter,  16  sample  holds,  a  variable 
quartz  clock  providing  sampling  frequencies  from  0.5  to  500 fc-ffz  resulting  in  an  overall  through¬ 
put  storage  capability  of  23,000  samples/sec.  The  storage  medium  was  a  cipher  tape  drive 
utilizing  IBM  compatible  phase  encoded  (PE)  formating  and  stored  the  data  as  two-byte  in¬ 
teger  words  in  2’s  complement  binary  on  a  10  inch  reel.  The  versatility  of  the  system  was 
realized  by  its  large  storage  capacity  (approximately  25  mega-words),  and  the  capability  of 
the  tape  medium  to  be  read  by  almost  any  tape  reader  that  has  the  ability  to  read  variable 
length  records  on  1600  bpi  10  track  1/2  inch  tape.  Figure  3.10  illustrates  the  overall  pressure 
measuring  system.  For  the  Reynolds  numbers  investigated,  a  sampling  frequency  of  4 kHz  was 
found  sufficient  to  contain  all  relevant  information.  The  overall  throughput  was  therefore  5 
x  4,000  =  20,000  samples/sec  which  was  effectively  handled  by  the  tape  drive’s  throughput 
processor. 

The  manufacturer’s  reported  specifications  on  these  transducers  are  generally  not  realiz¬ 
able  in  the  laboratory  due  to  the  additional  impedances  contributed  by  the  various  pieces  of 
instrumentation  used  for  amplification  and/or  filtering.  In  addition,  as  will  be  addressed  in 
more  detail  later,  due  to  the  nature  of  some  of  the  statistical  functions  necessary  to  be  cal¬ 
culated,  it  was  imperative  to  have  an  accurate  estimate  of  the  sensitivity  of  each  transducer 
during  each  respective  run.  In  order  to  satisfy  any  curiosity  about  bandwidth  or  amplitude 
limitations,  all  transducers  were  calibrated  in  a  small  J-ll  piston  actuated  standtube  calibra¬ 
tor  filled  with  water.  The  calibrator  could  be  operated  over  a  range  from  40 Hz  to  10kHz  and 
each  transducer  was  connected  to  all  expected  instrumentation  with  appropriate  amplification 
and  filter  bandwidths  preset.  These  small  transducers  are  very  fragile  hence,  susceptible  to 
destruction  by  frequent  removal  from  the  test  plate  for  subsequent  calibrations.  Therefore, 
since  the  flow  within  the  channel  is  considered  homogeneous  in  the  plane  of  the  wall,  it  was 
only  necessary  to  carefully  remove  a  single  transducer  for  calibration.  During  subsequent  runs, 
comparison  with  the  other  transducers  on  a  continuous  basis  afforded  updated  calibrations. 
Figure  3.11  shows  the  individual  calibrations  of  each  transducer  over  a  bandwidth  of  500Hz, 
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and  figure  3.12  shows  the  power  spectral  density  (PSD),  $(/)  of  three  representative  trans¬ 
ducers  during  typical  Reynolds  number  runs  of  Rh  =  10,000  and  15,000.  As  indicated,  other 
than  spurious  electronic  signals  seen  by  each  transducer,  virtually  identical  results  are  obtained 
throughout  the  examined  bandwidth. 


3.1.7  Software  Approach  for  Spectral  Computations 

All  transducer  data  was  recorded  simultaneously  in  time  and  transferred  to  a  VAX  780 
computer.  At  the  Naval  Research  Laboratory,  these  computers  act  as  front  end  processors  to  a 
CRAY  X-MP/24  super  computer.  The  data  were  transferred  to  the  CRAY  through  the  VAX 
front  ends,  and  software  was  written  to  analyze  the  time  dependent  data.  For  each  velocity  run, 
the  data  were  read  into  appropriate  arrays  and  the  respective  tranducer  sensitivities  applied 
immediately.  The  software  computed  the  power  and  cross  spectral  densities  and  their  Fourier 
transform  pair,  auto  and  cross  correlation  functions,  together  with  all  coherence  and  phase 
information  where  appropriate.  As  an  example  of  the  speed  of  the  CRAY  processing,  the 
software  upon  execution  could  process  all  data  from  one  Reynolds  number  run,  providing 
all  possible  combinations  of  frequency  spectra  and  time  correlation  functions  (as  indicated  in 
section  2.3)  for  more  than  12  transducers  in  approximately  2  1/2  minutes.  More  than  100  files 
could  be  generated  and  stored  simultaneously  for  each  velocity  run. 

The  choice  of  using  the  software  approach  was  predicated  on  the  ability  of  securing  almost 
immediate  turn-around  on  processed  data.  This  allowed  for  various  combinations  of  transform 
6ize  and  number  of  ensemle  averages  to  be  investigated  to  extract  as  much  relevant  information 
as  possible.  One  of  the  goals  of  this  experimental  6tudy  was  to  compare  a  physical  experiment 
with  a  computational  experiment.  This  software  method  of  analysis  provided  compatability 
with  the  processing  of  both  6ets  of  data.  As  a  final  comment,  the  cancellation  methodology 
outlined  in  the  next  section,  was  more  easily  applied  to  the  data  of  the  physical  experiment  by 
software  processing. 


3.1.8  Cancellation  Methodology 

Wambsganss  and  Zaleski  (1970)  first  introduced  a  temporal  subtraction  scheme  in  an  at¬ 
tempt  to  cancel  low  frequency  contaminating  noise  within  fully  developed  annular  flow.  Refine¬ 
ments  were  made  by  Wilson  et  al.( 1979)  and  the  present  author  in  1981.  If  the  contaminating 
noise  was  acoustic  in  origin  and  could  be  assumed  to  propagate  as  a  plane  wave,  this  temporal 
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subtraction  scheme  worked  quite  well,  provided  pressure  signals  were  dealt  with.  However, 
one  very  important  assumption  implicit  within  these  methods  but  “not”  explicitly  pointed  out 
in  this  past  work,  is  the  fact  that  “any”  cancellation  scheme  involving  the  addition  and/or 
subtraction  of  transducer  output  signals,  a  priori  assumes  that  the  individual  transducer  sen¬ 
sitivities  must  be  equal.  This  seems  trivial,  but  as  shown  below,  there  has  been  a  very  subtle 
misconception  in  the  past  work  related  to  transducer  output  and  sensitivity. 

First,  assume  a  simple  2  transducer  arrangement  as  depicted  in  figure  3.13  within  a  rect¬ 
angular  channel  or  along  a  flat  plate.  V  represents  the  velocity  in  the  streamwise  direction  and 
transducers  1  and  2  are  separated  transverse  to  the  flow  velocity.  The  distance  between  the 
transducers  indicated  as  /  is  very  large  (l  >>>  0)  such  that  the  signal  due  to  turbulence  at 
transducer  1  is  not  correlated  with  the  turbulence  at  transducer  2.  The  assumption  is  further 
made  that  the  contaminating  noise  whether  narrow  band  acoustic  in  origin,  or  of  a  wide  band 
random  phenomena,  is  fully  correlated  at  each  location  with  itself  but  not  with  the  signal  due 
to  turbulence.  Assume  that  the  pressure  signal  at  either  location  is  the  sum  of  the  signals  due 
to  noise.  ”n”  and  turbulence  "t”  as  represented  by  the  following  relation. 

Pj  -  b  +  tij  =  pt)  +  pn>  (3.11) 

Note,  the  key  point  being  made  is  that  equation  (3.11)  refers  to  "actual”  pressure  signals  and 
"not”  voltage  signals.  In  most  analyses,  the  temptation  is  strong  to  assume  that  sensitivity  is 
just  a  constant,  that  may  be  applied  at  any  point  in  the  processing  of  data  that  is  convenient 
for  the  researcher.  When  the  desired  output  deals  with  transformed  quantities  such  as  obtained 
from  Fourier  transforms,  this  is  not  necessarily  true.  Using  bj  and  ej  as  the  sensitivity  and 
voltage  respectively,  the  following  relations  can  be  formed. 


( total  pressure)  pj  =  bj  ej 

(3.12a) 

(pressure  due  to  turbulence)  pti  —  bj  et> 

(3.126) 

( pressure  due  to  noise)  pnj  =  bj  en> 

(3.12c) 

and, 

Pj  =  bj  ej  =  bj  et>  +  bj  en> . 

(3.13) 

In  the  laboratory,  the  voltage  signals  Cj,  etj ,  and  en.  axe  the  components  usually  considered 

rather  than  the  pressures  themselves.  In  fact  the  whole  field  of  transducer  technology  is  geared 
to  transforming  the  physical  entities  of  pressure,  force,  length  and  time  (to  name  a  few)  to  an 
easily  measured  voltage  or  current  signal  with  appropriate  sensitivity. 
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The  following  conventions  are  used.  Lowercase  letters  will  represent  the  signals  in  temporal 
space,  i.e.,  as  a  function  of  time,  and  capitalized  letters  will  refer  to  their  Fourier  transforms 

into  spectral  space.  The  braces  <  . >  represent  the  expected  value  over  all  ensembles  and 

the  superscript  •  represents  the  complex  conjugate  of  a  Fourier  transform.  ipp  and  will 
represent  the  auto  or  cross  spectrums  of  the  total  pressure  and  voltage  signals  respectively.  As 
an  example,  the  auto  spectrum  of  the  pressure  at  a  transducer  ‘j’  can  be  represented  as, 

*p,p,  =  <  p;  Pi  >  «  <  (Tj  +  Nj)m  (Tj  +  Nj)  > 

«  <  (T-Tj)  +  (n;Nj)  > 

~  ^p»,p>,  ^p»t p«,  (3.14 ) 

Assuming  the  transducer  signals  are  due  to  voltages  only,  then  the  total  sensor  signal  is, 


€i  =  £t,  +  en> 

and  the  Fourier  transform  in  spectral  space  results  in  the  following. 


(3.15) 


*e>e,  =  <  E]  Ej  >  =  <  (Et,  +  En,)m  (. Et ,  +  £„,)  > 


Furthermore,  it  can  be  shown  that, 


*  <  (£•£«,)  +  (£;,£„,)  > 
s  ^ et)  etj  "h 


(3.15) 


* 


and, 


hence, 


n,n> 


*PiPi  -  ®  [(*«,i«,i)  + 


(3.16) 


The  temporal  subtraction  scheme  based  on  the  assumptions  given,  consists  of  performing  a 
temporal  subtraction  between  the  signals  from  transducers  1  and  2  and  then  correlating  the 
resulting  subtracted  signal  with  either  transducer.  In  other  words,  the  power  spectrum  desired 
can  be  represented  as,  ♦(p,-pa)p,.  Bear  in  mind  that  the  assumptions  state  that  the  turbulence 
at  locations  1  and  2  are  not  correlated  with  each  other,  and  the  noise  is  correlated  with  itself 
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and  at  either  location,  but  not  with  the  turbulence.  Performing  the  required  algebra  results  in 
the  following. 

^(Pi-Pa)Pi  =  <  (ft  ~  ft)  ft  > 

=  <  (ft  +  N!  -  ft  -  N2)m  (ft  +  Ni)  > 

=  <  r; ft  +  n;n2  -  n;nx  > 

The  last  two  terms  are  assumed  equal  and  therefore  cancel  resulting  in, 

^(Pi-Pa)Pi  =  ^P<jPf,  •  (3-1 «) 

If  the  same  analysis  is  performed  on  the  voltage  signals,  with  proper  substitution  of  the  sensi¬ 
tivities,  the  following  equation  is  obtained. 

^(Pj-Pa)Pi  =  ^P«iP«i  +  (l  —  ^p»p»  (3.18) 

As  indicated,  equations  (3.17)  and  (3.18)  are  the  same,  “if  and  only  if”  the  sensitivities  bx 
and  i>2  are  equal.  It  makes  no  difference  whether  the  subtraction  is  performed  in  "temporal" 
or  “spectral”  space.  The  problem  arises  due  to  the  researcher  conducting  analyses  in  which 
the  real  signals  (with  proper  units)  are  being  dealt  with  on  paper,  whereas  in  the  laboratory, 
voltage  signals  are  the  entities  being  recorded. 

The  above  monolith  seems  simple,  but  in  reality  poses  enormous  possiblities  of  error  when 
dealing  with  multi-sensor  arrays.  For  this  experimental  study,  the  sensitivities  of  all  five  trans¬ 
ducers  would  have  to  be  maintained  constant  and  equal  throughout  every  experimental  run. 
This  was  virtually  impossible  due  to  the  finite  time  allowed  for  each  experimental  run  as  well 
as  the  limited  instrumentation  available.  However,  if  the  transducers  have  been  consistently 
calibrated,  such  that  precise  knowledge  of  their  sensitivity  is  known  at  all  times,  then  the 
software  approach  to  processing  the  data  as  presented  above  provides  a  convenient  means  for 
applying  the  sensitivities  to  each  channel  of  data  as  it  is  read  back  into  the  analysis  program. 
In  these  circumstances,  the  actual  “pressures”  are  being  manipulated  without  introducing  any 
errors. 

To  backtrack  slightly,  the  temporal  subtraction  scheme  outlined  above  works  quite  well, 
provided  the  contaminating  noise  is  acoustic  in  origin  and  propagates  as  a  plane  wave  across 
any  cross-sectional  area  of  the  channel.  This  assumption  is  quite  restrictive  and  in  particular, 
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doesn’t  provide  a  means  of  cancelling  random  vibrational  energy  that  may  interfere  with  the 
transducer  signal.  Horne  and  Handler  (1989)  have  presented  a  unique  cancellation  scheme 
utilizing  the  concept  of  a  “Wiener”  filter  (Wiener,  1949).  This  filter  was  designed  originally  to 
cancel  wide  band  frequency  signals  in  electrocardiography  and  the  broad-band  interference  in 
the  sidelobes  of  an  antenna  array  and  assumes  that  a  signal  free  reference  is  available  which 
represents  the  contaminating  noise.  In  applying  the  filter  technology  to  the  measurement  of 
wall  pressure  fluctuations,  certain  assumptions  provided  the  author  with  a  unique  derivation 
of  error  relationships  which  precisely  dictated  the  maximum  error  possible  in  the  calculation 
of  auto  and/or  cross  spectral  properties. 

To  briefly  summarize,  transducer  1  is  considered  to  be  the  signal  free  reference  in  this 
experimental  study.  It  is  assumed  to  be  separated  from  the  other  transducers  in  the  linear 
array,  a  distance  sufficient  to  insure  zero  correlation  of  turbulence.  In  addition  it  is  assumed 
that  the  channel  represents  truly  homogeneous  turbulent  flow,  such  that  the  auto  spectral 
densities  at  each  location  are  identical.  As  a  final  assumption,  the  noise  at  any  location  is 
assumed  to  be  uncorrelated  with  the  turbulence  at  neighboring  locations.  Note,  no  stipulation 
has  been  put  upon  the  propagation  characteristics  of  the  noise,  whether  it  is  acoustic  in  origin, 
single  frequency  or  random  by  nature.  The  only  assumption  made,  is  that  the  noise  represented 
at  each  location  has  “finite”  coherence  with  the  noise  located  at  transducer  1,  and  that  this 
coherence  is  measurable.  The  results  of  the  cancellation  scheme  are  presented  below  for  the 
auto  and  cross  spectral  densities  $p<Pl  and  $p,iPlj  of  the  turbulence  respectively.  The  reader 
may  refer  to  Horne  and  Handler  (1989)  for  a  more  detailed  presentation  of  the  method. 

Auto  Spectrum  — *  *p<iPli  =  $p>Pi  (1  -  7?i)  (3.19) 

t  =  2, 3, 4, 5 

Cross  Spectrum  — ►  *p,.p  =  $PlP,  -  (3.20) 

vpip\ 

i  =  2  only,  and  j  =  3,4,5  only. 

The  power  spectral  densities  in  equations  (3.19)  and  (3.20)  are  computed  from  the  sensor 
signals  themselves  and  the  order  of  the  indices  is  imperative.  Note,  other  than  the  denominator 
in  the  second  term  of  equation  (3.20),  the  quantities  computed  are  complex  variables  with  $*lPi 
representing  the  complex  conjugate  of  this  cross  spectrum  respectively. 
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The  above  two  relations  were  used  extensively  in  this  study  to  obtain  true  measurements 
of  the  power  spectral  densities,  coherence  and  phase  for  the  physical  experimental  data.  The 
results  thus  obtained  provide  new  information  not  yet  published  and  confirm  existing  theories 
on  the  ‘u>2’  dependence  at  low  frequencies  for  the  wall  pressure  fluctuation  amplitude. 


3.2  Computational  Experiment 

The  wavenumber-frequency  spectrum  of  turbulence-induced  wall  pressure  fluctuations  is 
of  both  fundamental  as  well  as  practical  interest.  As  mentioned  previously,  this  spectrum  may 
provide  insight  into  the  fundamental  structure  of  turbulent  flow  in  the  near-wall  region.  From 
a  practical  point  of  view,  it  may  allow  prediction  of  the  small  amplitude  vibration  response  of  a 
structure  to  turbulence-induced  wall  pressure  fluctuations.  In  the  past,  however,  this  spectrum 
has  proven  somewhat  elusive  due  to  its  difficulty  in  being  measured.  The  numerical  simulation 
experiments  offer  significant  promise  in  this  area. 

3.2.1  Methodology  of  Solution 

The  governing  equations  for  the  current  DNS  work  at  NRL  begin  with  equation  (2.6) 
written  in  vector  notation  as, 

H  +  u-VU  =  -V-  +  g  +  J/V2U  (3.21) 

ot  p 

By  using  the  vector  identity, 

U  VU  =  -  Uxfi  (3.22) 

the  time  dependent  Navier-Stokes  equations  in  the  so-called  rotational  form  are  obtained  as 
given  by, 

~  =  Uxfl  -  +  g  +  x/V*U  (3.23) 

From  continuity,  the  following  incompressibility  relation  is  obtained, 

VU  =  0.  (3.24) 

Equations  (3.23)  and  (3.24)  form  a  non-linear  set  of  4  equations  for  the  four  unknown  variables, 
ti,  v,  w  and  p.  In  these  two  equations,  v  equals  the  kinematic  viscosity,  ft  =  V  x  u  equals 
the  vorticity,  and  g  represents  a  known  body  force.  A  pseudo-spectral  solution  technique  is 
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used  to  solve  for  the  velocity  and  pressure  fields.  Velocity  components  are  expanded  in  Fourier 
series  in  the  streamwise  and  spanwise  directions  and  Chebychev  polynomials  in  the  wall-normal 
direction.  All  spatial  derivatives  are  evaluated  in  the  transform  domain,  after  which  inverse  fast 
Fourier  transforms  are  employed  to  transform  the  solution  for  one  time  step  to  the  physical 
domain.  The  solution  is  moved  forward  to  the  next  time  step  using  the  Adams  Bashforth 
time  stepping  algorithm.  Fast  Fourier  transforms  are  then  used  to  move  this  updated  solution 
back  into  the  transform  domain,  where  spatial  derivatives  are  again  evaluated  spectrally.  In 
carrying  out  the  solution,  the  flow  was  assumed  to  be  initially  laminar.  After  time  t  =  0,  a 
large  unstable  disturbance  was  added  to  the  laminar  velocity  field.  The  establishment  of  fully- 
developed  conditions  required  approximately  4000  time  steps.  The  pressure  at  both  upper 
and  lower  walls  for  all  points  in  the  computational  grid  was  obtained  from  equation  (2.7), 
and  presented  to  the  author  by  his  colleagues  for  analysis.  The  reader  is  referred  to  Handler, 
et  al.(  1989)  for  a  more  definitive  description  of  the  methodology.  In  any  event,  this  investigation 
seeks  only  to  compare  and/or  critique  the  results  from  a  typical  calculation  (particularly  the 
wall  pressure  fluctuations),  rather  than  comment  on  the  accuracy  of  the  technique  used. 


3.2.2  Description  of  Physical  Model 

For  the  data  to  be  presented,  computation  continued  until  a  fully  developed  channel  flow 
was  obtained.  The  computational  volume  included  16  points  within  the  flow  direction  x,  33 
points  in  the  wall  normal  direction  y,  and  64  points  in  the  transverse  direction  2,  lying  in  the 
wall.  Hence,  on  both  the  upper  and  lower  walls  of  the  channel,  there  were  a  total  of  2  x  (16  x  64 ) 
=  2048  separate  realizations  for  the  wall  pressure.  At  each  of  these  2048  computational  points, 
512  full  realizations  in  time,  of  velocity  and  pressure  were  calculated  and  stored  on  tape.  These 
given  parameters  somewhat  limited  the  spectral  analysis.  Since  only  512  time  points  were 
computed,  this  provided  only  257  points  in  the  frequency  domain  to  be  resolved.  In  addition, 
each  computation  point  was  treated  as  an  ensemble,  with  all  points  added  together  to  give 
2048  averages.  In  a  statistical  evaluation,  this  does  not  translate  into  a  true  2048  ensemble 
average.  Many  of  the  points  used  are  well  within  the  correlation  length  for  the  wall  pressure  at 
the  wall  and  therefore  are  not  true  statistically  independent  records,  as  described  in  chapter  2. 
The  calculation  was  performed  at  a  Reynolds  number  based  on  shear  velocity  and  half-width 
of  the  channel  as, 

Rm  *=  =  125.  (3.25) 
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The  physical  experiment  had  an  R‘  value  ranging  from  approximatly  500  to  1000,  or  4  to  8 
times  larger.  In  contrast,  most  of  the  experimental  work  to  date,  has  values  two  orders  of 
magnitude  or  greater,  on  the  order  of  Rm  greater  than  or  equal  to  5,000. 
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FIGURE  3.4.  TEST  SECTION  CONFIGURATION 
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IE  3.6a.  TRANSDUCER  CALIBRATION  FOR  BULK  VELOCITY 
CONTROL 


FIGURE  3.6b.  TRANSDUCER  CALIBRATION  FOR  CHANNEL 
PRESSURE  GRADIENT 
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FIGURE  3.7.  CHANNEL  PRESSURE  GRADIENT  USING  MULTIPLE 
TAPS  FOR  A  RANGE  OF  BULK  VELOCITIES 


REYNOLDS  NUMBER  (Rdh) 

FIGURE  3.8.  FRICTION  COEFFICIENT  FOR  RECTANGULAR  CHANNEL 
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FIGURE  3.9.  SCHEMATIC  OF  LDV  ORIENTATION  w/RESPECT  TO 
TUNNEL  SIDEWALLS 
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FIGURE  3.10.  BLOCK  DIAGRAM  OF  WALL  PRESSURE  DATA  ACQUISITION 


40.0 


FIGURE  3.12.  POWER  SPECTRAL  DENSITY  OF  THREE  TRANSDUCERS  AT 
TWO  REPRESENTATIVE  REYNOLDS  NUMBERS 
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|  U  (FLOW  DIRECTION) 

| - ► 

TRANSDUCER  No.  1  TRANSDUCER  No.  2 

FIGURE  3.13.  TRANSDUCER  ORIENTATION  FOR  CANCELLATION 


CHAPTER  4  EXPERIMENTAL  RESULTS  FOR  VELOCITY 

4.1  Introduction 


This  chapter  presents  the  results  for  the  streamwise  mean  velocity  and  it’s  fluctuating  rms 
value  U  and  ti  respectively,  for  the  physical  and  computational  experiments.  The  objectives  of 
these  measurements  were  to  establish  that  both  experiments  represented  a  fully  developed  tur¬ 
bulent  channel  flow  for  the  Reynolds  numbers  investigated.  Furthermore,  it  was  also  important 
to  establish  that  the  large  aspect  ratio  of  the  channel  for  the  physical  experiment  ensured  that 
no  transverse  variations  in  the  turbulent  structure  interfered  with  the  flow’s  two-dimensionality. 
As  pointed  out  in  chapter  3,  the  LDV  system  was  limited  to  a  single  component  in  the  flow 
direction  since  the  side  wall  plexiglass  permitted  access  to  only  two  beam.  Although  limited, 
it  was  felt  that  these  measurements  were  sufficient  to  establish  that  the  physical  experiment 
demonstrated  fully  developed  turbulent  flow. 

The  data  obtained  for  the  computational  experiment  included  profiles  of  the  mean  and 
fluctuating  components  of  velocity  as  a  function  of  the  distance  normal  to  the  wall.  These  data 
were  averaged  over  the  time  period  of  the  computational  run  as  well  as  the  entire  computational 
box.  The  justification  for  this  averaging  technique  lies  in  the  assumption  that  if  the  flow  is 
considered  fully  developed,  then  horizontal  planes  represent  statistically  independent  samples 
of  homogeneous  turbulent  flow.  Due  to  the  finite  correlation  lengths  which  exist  in  both 
streamwise  and  transverse  directions  for  turbulent  flow,  statistical  independence  is  not  strictly 
met  thereby  introducing  a  source  of  random  error.  These  profiles  were  represented  as  t/(y)/u* 
and  u(y)/u*  respectively,  as  a  function  of  Y+  =  yu‘  jv.  As  mentioned,  the  governing  Reynolds 
number  for  the  computational  experiment  was  the  value  Rm  =  125,  which  was  based  on  the 
shear  velocity  «*,  channel  half- height  h,  and  kinematic  viscosity  v,  respectively.  These  three 
parameters  provide  all  that  is  necessary  to  derive  other  profile  characterizations  which  are 
outlined  in  table  4.1. 

Profiles  were  taken  in  the  rectangular  channel  of  the  physical  experiment  at  three  trans¬ 
verse  locations  labeled  a s  A,  B  and  C.  These  profiles  were  taken  a  number  of  times  to  insure 
their  repeatability.  Position  B  was  approximately  located  on  the  channel’s  center  line,  at  8 
inches  from  the  near  side  wall  and  were  those  measurements  compared  with  the  computational 
experiment.  Four  Reynolds  numbers  were  investigated  at  each  station  and  were  character¬ 
ized  a £  Rh  =  Uavgh/u  =  5000,10000,15000  and  20000  respectively.  Pertinent  characteristics 
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describing  representative  profiles  at  each  location  are  presented  in  table  4.2.  It  should  be 
noted  that  the  Reynolds  numbers  of  the  physical  experiment  were  significantly  higher  than  the 
computational  experiment. 


Table  4.1  Parameters  for  Computational  Experiment 


Given  : 

U+  =  U (y)/ U*.  N  vklum 

Y+  —  y\lm /v,  N  values 

ttt  =  u(y)/um,  N  vftluM 

Rm  =  vTh/u 

Derived : 

Vo/*  =  K 

u(y)/U(y)  =  «;/r; 

u(y)/u0  =  vfiv% 

vjh  =  y;/r‘ 

10. -«»)]/«■  »  v$-vf 
«0r)/V.  =  «+/PJ 

Rk  =  V%  R‘ 

It  should  be  noted  that  the  velocity  measurements  and  wall  pressure  measurements  for  the 
physical  experiment  were  not  taken  at  the  same  time.  Some  of  the  same  instrumentation  was 
used  in  the  characterization  of  velocity  and  pressure  and  precluded  simultaneous  measurement 
in  some  cases.  In  addition,  the  complexity  of  obtaining  both  types  of  data,  which  involved 
constant  monitoring  of  a  number  of  controlling  parameters  (indicated  by  various  readout  de¬ 
vices),  permitted  only  one  type  of  data  (velocity  or  wall  pressure)  to  be  obtained  at  any  one 
time.  Therefore,  the  velocity  data  was  taken  first  at  the  Reynolds  numbers  specified.  Later, 
when  the  wall  pressure  data  was  investigated,  it  was  found  that  the  lower  velocity  runs  pro¬ 
vided  pressure  intensities  of  such  low  magnitude  that  the  pressure  sensors  could  not  resolve  the 
pressure  signature  adequately,  due  to  their  limited  dynamic  range.  It  will  be  shown  in  chapter 
5,  that  only  the  larger  three  Reynolds  numbers  could  be  investigated  along  with  an  additional 
higher  Reynolds  number  at  =  25000. 
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Table  4.2  Governing  Parameters  for  Physical  Experiment 


Rh 

station 

u’/V0 

Rm 

Uavg(m/sec) 

5,000 

A 

.0509 

373 

0.457 

5,000 

B 

.0511 

375 

0.448 

5,000 

C 

.0503 

369 

0.448 

10,000 

A 

.0482 

560 

0.719 

10,000 

B 

.0480 

563 

0.719 

10,000 

C 

.0477 

559 

0.722 

15,000 

A 

.0460 

839 

1.158 

15,000 

B 

.0451 

836 

1.125 

15,000 

C 

.0454 

834 

1.131 

20,000 

A 

.0446 

1015 

1.417 

20,000 

B 

.0447 

1015 

1.396 

20,000 

C 

.0443 

1010 

1.406 

4.2  Mean  Velocity  Profiles 

The  mean  velocity  profiles  U(y)/U0,  as  a  function  of  distance  from  the  wall  y/h  repre¬ 
senting  all  four  Reynolds  numbers,  are  shown  in  figures  4.1  to  4.4  for  the  physical  experiment. 
For  comparison  purposes,  a  laminar  flow  profile  for  the  same  flow  configuration  is  also  shown. 
In  comparison  with  laminar  flow,  all  profiles  for  each  Reynolds  number  demonstrate  the  fuller 
shape  which  is  characteristic  of  typical  fully  developed  turbulent  flows. 

Observation  of  the  data  show  that  no  variation  across  the  channel  is  apparent  with  excellent 
agreement  among  the  three  stations.  An  interesting  anomaly  may  be  noted  in  relation  to 
resolution  at  the  wall.  As  one  examines  figures  4.1  thru  4.4,  resolution  with  respect  to  the 
wall  normal  coordinate  ‘y’  appears  to  get  worse  as  the  Reynolds  number  increases.  As  the 
velocity  increases  in  the  channel,  so  does  the  data  rate  of  detected  particles  passing  through 
the  probe  volume.  However,  larger  pan  tides  which  obliterate  the  signal  from  the  probe  volume 
are  also  more  numerous  and  create  what  is  termed  “shot  noise”  to  the  photo  multipliers,  thus 
redudng  the  signal  to  noise  ratio.  As  the  flow  moves  from  low  to  high  Reynolds  number,  one 
obtains  first  less  noise  but  low  data  rate  and  Anally  high  data  rate,  but  too  much  “shot  noise” 
or  flare  up;  hence  moving  away  from  the  wall  is  required.  This  phenomenon  is  indicated  in 
flgure  4.4  by  showing  the  dosest  distance  to  the  wall  available  as  y/h  =  .04  in  comparison  with 
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a  value  of  .01  in  figure  4.2.  Secondly,  by  examining  each  station,  .4,2?  or  C,  one  may  note  B 
appears  to  give  the  best  resolution  regardless  of  Reynolds  number.  This  can  be  explained  by  the 
following  argument.  Station  A  is  close  to  the  side  wall  and  allows  reflections  from  the  plexiglass 
to  interfere  with  the  signal.  As  the  beams  are  traversed  towards  station  C,  a  greater  length 
of  beam  travels  through  the  flow  environment  and  requires  the  photo  multipliers  and  related 
instrumentation  to  be  more  discriminating  in  validating  particles.  Both  situations  reduce  the 
signal  to  noise  ratio  of  the  LDV  system  resulting  in  poorer  resolution.  By  examining  each  of 
the  figures,  station  B  appears  to  represent  the  best  compromise  between  both  anomalies. 

In  order  to  establish  more  resolutely  the  two-dimensional  nature  of  the  flow,  the  laser 
probe  volume  was  positioned  at  the  two  distances  from  the  wall  of  y/h  =  0.2  and  1.0.  At  each 
position,  the  probe  volume  w'as  traversed  across  the  channel  from  approximately  2  thru  12  to 
14  inches  from  the  near  wall  for  the  two  middle  Reynolds  numbers  of  Rh.  =  10000  and  15000. 
These  profiles  are  presented  in  figures  4.5  and  4.6  with  the  location  of  stations  A,B  and  C 
also  indicated.  No  gross  interference  from  the  side  wall  can  be  detected  until  the  probe  volume 
approached  approximately  4  inches  to  the  wall  as  indicated  by  the  results  at  Rh  =  10000 
in  both  figures.  The  slight  variations  noticed  are  due  to  the  slight  movement  of  the  probe 
volume  as  the  beams  move  through  the  side  wall.  As  the  beams  pass  through  the  plexiglass  at 
different  locations,  small  imperfections  cause  the  beam  to  diffract  and  change  its  path.  This 
results  in  the  probe  volume  moving  to  a  slightly  different  location  from  that  expected.  This 
effect  becomes  more  pronounced  the  closer  the  beams  are  moved  towards  the  wall  due  to  the 
presence  of  a  higher  shear  gradient.  This  is  indicated  by  the  slight  waviness  in  the  profiles 
which  is  somewhat  more  noticible  in  figure  4.5  than  figure  4.6,  as  expected. 

Figure  4.7  takes  the  results  from  station  B  at  all  Reynolds  numbers  and  compares  them 
with  those  from  the  computational  experiment.  In  this  figure,  the  measurements  of  Hussain  and 
Reynolds  (1975)  and  Laufer  (1951)  are  indicated  where  the  Reynolds  numbers  are  appropriate. 
Two  conclusions  may  be  drawn  from  the  comparison.  First,  the  data  of  the  physical  experiment 
are  in  excellent  agreement  with  other  reported  investigations  at  similar  Reynolds  numbers.  The 
data  of  Laufer  and  Hussain  were  also  taken  in  rectangular  channels  with  aspect  ratios  of  12  and 
18,  respectively.  Secondly,  the  Reynolds  number  dependence  corresponding  to  outer  variable 
scaling  i6  well  indicated. 

The  data  taken  along  the  channel  centerline  (station  B)  is  compared  with  the  compu¬ 
tational  results  in  figure  4.8  using  inner  variable  scaling  for  the  classical  "Law  of  the  Wall" 
profiles.  The  solid  line  represents  the  following  equation  derived  by  Spalding  (1961), 
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y+  =  U+  +  exp{-KB)  \exp(*B)  -  1  -  kU+  -  -l-p~  -  ^Up  ]  (4.1) 

with  the  constants  k  and  B  equal  to  0.40  and  5.5  respectively.  Also  6hown  in  figure  4.8  are 
the  different  regions  of  a  typical  boundary  layer  governed  by  a  wall-bounded  turbulent  flow. 
The  linear  sublayer  (U+  =  V+)  covers  the  range  0  <  Y+  <  5,  the  buffer  layer  covers  the 
range  approximately  5  <  Y+  <  40  and  the  log-law  region  from  40  <  y+  to  y  <  0.2h.  These 
three  regions  cover  the  inner  layer  of  the  boundary  layer.  The  outer  or  "velocity  defect"  region 
exists  for  y/h  >  0.2.  All  the  profiles  in  figure  4.8  collapse  on  the  solid  line  and  appear  to 
merge  smoothly  from  the  log-law  region  to  the  inner  sublayer  region.  As  shown,  the  limited 
spatial  resolution  due  to  the  instrumentation,  only  gave  a  few  points  in  the  buffer  layer  for  the 
mid-range  Reynolds  numbers,  but  qualitatively  the  overall  behavior  is  in  excellent  agreement 
with  the  classical  results.  At  larger  values  of  y+,  there  exist  differing  values  depicting  the 
onset  of  the  wake  or  “velocity  defect"  region  for  each  profile.  The  computational  data  appears 
to  terminate  at  y+  =  90,  with  V+  values  of  200,400,600  and  800  for  the  physical  Reynolds 
number  flows  of  5000, 10000, 15000  and  20000,  respectively. 

To  better  illustrate  this  defect  region,  the  profiles  of  all  five  flows  are  plotted  with  outer 
variables  in  figure  4.9  as  [U0  -  U{y)]/u*  versus  y/h.  The  defect  region  is  well  indicated, 
beginning  at  y/h  approximately  equal  to  0.2,  with  the  profiles  becoming  more  concave  and 
fuller  with  decreasing  Reynolds  number  as  the  wall  is  approached. 


4.3  Streamwise  Velocity  Fluctuations 

The  turbulence  intensity  profiles  at  station  B  for  all  Reynolds  numbers  of  the  physical 
experiment  are  shown  in  figures  4.10  and  4.11  along  with  the  computational  results.  These 
data  are  plotted  as  a  function  of  the  centerline  velocity  U0  in  figure  4.10,  and  the  local  velocity 
U(y)  in  figure  4.11.  A  definite  Reynolds  number  dependence  is  demonstrated,  with  the  peak 
value  increasing  and  moving  away  from  the  wall  with  decreasing  Reynolds  number.  The  data 
presented  in  both  figures  show  good  agreement  with  similar  experiments  conducted  in  other 
large  aspect  ratio  channels. 

The  rms  values  of  fluctuating  velocity  are  plotted  with  inner  variable  scaling  in  figure  4.12 
as  17+  =  V / ti*  versus  Y+  =  yum/v.  The  data  of  Kreplin  and  Eckelmann  (1979)  and  Clark 
(1968)  aic  oL~>vn  for  comparison.  The  overall  behavior  exhibits  trends  similar  to  other  classic 
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studies  in  that  the  results  increase  monotonically  away  from  the  wall,  peak  at  Y+  values  of 
between  13  to  18,  and  then  slowly  decrease  to  the  value  measured  at  the  channel  centerline. 
The  rate  at  which  the  results  decrease  is  Reynolds  number  dependent  as  shown  in  the  literature. 
It  should  be  noted  that  the  peak  value  measured  is  limited  by  the  transducer  used.  In  most 
cases,  hot  films  or  hot  wires  are  used  in  which  case  the  length  of  the  wire  (or  film)  limits  the 
scale  of  fluctuations  capable  of  being  measured.  Kreplin  and  Eckelmann  used  hot-film  sensors 
in  oil  for  which  the  effective  sensor  length  was  as  low  as  two  viscous  units.  In  the  present 
study,  the  LDV  represents  a  non-intrusive  probe  (producing  no  obstruction  to  the  flow),  60  it 
can  be  argued  that  the  probe  volume  diameter  is  the  limiting  dimension.  The  beam  diameter 
of  .005  inches  corresponds  to  approximately  10  viscous  units  for  the  larger  Reynolds  number 
and  should  not  be  a  limiting  factor. 

The  peak  value  of  about  2.8  to  3.0  obtained  in  the  physical  experiment  is  well  within  the 
range  measured  by  other  investigators  and  occurs  at  a  y+  value  of  about  15.  The  computational 
results  peak  in  the  same  approximate  location  from  the  wall  but  the  peak  is  somewhat  lower, 
reaching  a  value  of  only  2.5.  Kreplin ’s  data  had  an  R"  value  of  approximately  390  which  was 
close  to  those  values  of  the  physical  experiment  at  the  low  Reynolds  number  of  5,000  but 
some  3  times  larger  than  the  value,  Rm  =  125  of  the  computational  experiment.  Although 
qualitatively  both  experimental  data  presented  here  agree  reasonably  well  with  that  reported 
in  the  literature,  either  the  computational  experiment  demonstrates  a  Reynolds  number  effect 
not  shown  before,  or  has  not  yet  reached  full  development. 

The  results  presented  in  figures  4.10  and  4.11  support  this  argument.  As  mentioned 
previously,  scaling  the  fluctuations  with  outer  and  local  variables  should  show  a  peak  value 
increasing  as  the  Reynolds  number  decreases.  Relative  to  the  physical  experiment  and  the  data 
of  Laufer,  Clark  and  Kreplin,  the  peak  in  the  computational  experiment  is  somewhat  lower  than 
the  trends  would  seem  to  indicate.  The  computational  grid  may  not  have  been  of  a  sufficient 
length  nor  the  code  run-time  of  sufficient  duration  to  fully  develop  this  profile.  The  assumption 
of  statistical  independence  mentioned  in  the  previous  section  may  also  be  contributing  to  this 
effect.  Thi6  becomes  an  important  factor  because  the  source  terms  contributing  the  most  to 
the  wall  pressure  field  occur  closest  to  the  wall  at  Y+  values  less  than  40.  One  could  well  argue 
that  the  simulated  experiment  has  not  reached  its  fully  developed  state  in  a  spatial  sense,  or  is 
still  experiencing  transition.  These  ideas  will  be  addressed  further  in  chapter  5  where  the  data 
for  the  wall  pressure  field  will  be  presented. 


Hussain  (1983)  discussed  the  important  role  which  the  initial  condition  plays  in  the  de¬ 
velopment  of  a  turbulent  boundary  layer  and  the  necessary  conditions  for  a  fully  developed 
turbulent  profile.  For  channel  flow,  these  conditions  can  be  represented  as;  1)  a  mean  profile 
of  the  proper  shape,  2)  a  mean  profile  indicating  a  log-law  and  defect  region  of  proper  extent 
when  plotted  with  respect  to  inner  variables,  3)  a  peak  in  the  turbulence  intensity  level  of 
tt/ti*  =  2.5  ±  10%  at  y+  =  15,  and  4)  an  intensity  profile  which  monotonically  decreases  to  the 
centerline  value.  Based  on  these  characteristics,  figures  4.1  thru  4.12  clearly  show  the  physi¬ 
cal  experiment  to  be  fully  developed.  The  intensity  values,  although  somewhat  high  and  low 
for  the  physical  and  computational  experiments  respectively,  are  certainly  within  the  ranges 
outlined  above  by  Hussain.  A  more  detailed  critique  of  the  computational  data  was  given  by 
Handler  el  al. (1989)  and  will  be  addressed  in  chapter  5  in  relation  to  the  wall  pressure. 


4.4  Shear  Velocity 


A  parameter  of  critical  importance  in  evaluating  the  development  of  turbulence  profiles  is 
the  shear  velocity  ‘u“\  An  additional  criterion  to  those  outlined  by  Hussain  (1983)  in  relation  to 
a  channel  flow,  is  a  linearly  decreasing  static  pressure  throughout  the  fully  developed  flow  region 
thereby  exhibiting  a  constant  pressure  gradient.  In  most  external  boundary  layer  experiments, 
the  shear  velocity  is  not  measured  but  rather  inferred  by  forcing  the  mean  velocity  to  conform  to 
the  “law  of  the  wall”,  given  by  equation  (4.1).  In  addition,  recognizing  that  the  viscous  sublayer 
close  to  the  wall  must  result  in  a  linear  profile,  the  shear  velocity  is  sometimes  calculated  from 
the  slope  at  the  wall  from  the  relation, 


»=o 


(4.2) 


Each  parameter  inferred  or  otherwise  determined  by  assumptions  about  the  flow,  incurs  the 
risk  of  introducing  additional  errors.  For  the  present  physical  experiment,  the  pressure  gradient 
was  measured  directly  as  outlined  in  chapter  3.  For  each  Reynolds  number,  the  static  pressure 
gradient  differential  was  recorded  in  time  during  the  full  experimental  run  and  averaged  with  the 
use  of  equation  (2.2).  The  value  of  shear  velocity  for  each  Reynolds  number  is  plotted  in  figure 
4.13  in  terms  of  R*  =  umh/v  versus  the  channel  center-line  Reynolds  number  Rh  =  U0h/i>. 
The  numerical  experimental  value  of  125  is  also  shown  and  the  solid  line  represents  the  data 
of  Hussain  (1975)  which  was  found  to  correlate  with  the  relation, 


Rm  =  0.1079  (Ra).0  9U 


(4.3) 
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The  results  of  both  the  computational  as  well  as  the  physical  experiments  are  in  excellent 
agreement  with  this  correlation. 

In  summary,  the  results  presented  from  both  experiments  suggest  the  following  conclusions: 

o  The  mean  velocity  profiles  of  the  physical  experiment  demonstrate  a  fully 

developed  two-dimensional  channel  flow  throughout  the  core  region  of  the  duct. 

o  The  turbulence  intensity  profiles  with  respect  to  outer  and  local  variables, 
show  a  Reynolds  number  dependence  consistent  with  other  investigations  for 
the  physical  experiment. 

o  The  peak  in  the  rms  level  u/u*  is  approximately  2.8  to  3.3  and  occurs  at 
Y+  =  15,  in  good  agreement  with  accepted  values  in  the  literature. 

o  The  location  of  the  channel  test  section  shows  an  ideal  constant  pressure 
gradient  resulting  in  constant  values  for  shear  velocity  for  each  Reynolds 
number  of  the  physical  experiment. 

o  The  computational  experiment  demonstrates  similar  agreement  for  the  mean 
velocity,  but  shows  somewhat  retarded  development  with  respect  to  the 
fluctuating  component  u(y). 
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Station  Across  the  Ch  anne I  . 


Station  Rcross  the  Channel. 
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FIGURE  4.4  Velocity  Profile  at  Rh  =  20,000  at  Each 

Station  Rcross  the  Channel. 
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FIGURE  4.8  Universal  Velocity  Pr 
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FIGURE  4.10  RMS  Profile  Versus  Outer  Variables 
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FIGURE  4.13  Friction  Velocity  Versus  Reynolds  Number. 


CHAPTER  5 


EXPERIMENTAL  RESULTS  FOR  PRESSURE 


5.1  Introduction 

This  chapter  presents  a  comparison  of  the  results  obtained  for  the  wall  pressure  field  for 
both  the  physical  and  computational  experiments  which  is  the  primary  focus  of  this  study. 
There  are  a  number  of  valid  reasons  for  examining  duct  flows  at  low  Reynolds  numbers.  First 
(and  foremost),  the  recent  theoretical  work  related  to  Direct  Numerical  Simulation  (DNS)  ex¬ 
periments  has  generated  excitement  in  the  research  community  in  that  a  great  deal  of  informa¬ 
tion  appears  to  be  available  without  the  overhead  of  maintaining  large,  complicated  laboratory 
facilities.  However,  this  excitement  has  also  introduced  a  cloudy  issue,  in  that  successful  eval¬ 
uation  of  a  computational  experiment  requires  accurate  comparison  with  an  identical  physical 
experim mt  with  all  parameters  matching,  before  addressing  issues  not  yet  resolved.  This  means 
that  not  only  should  the  kinematic  properties  of  the  flow,  such  as  mean  velocity,  turbulence 
intensity  and  Reynolds  stress  profiles  be  in  agreement,  but  both  the  spatially  averaged  “scalar” 
properties  such  as  the  total  wall  shear  stress  t^u,  prmi/Twaii,  PrmtlQ,  and  the  more  detailed 
spectral  characteristics  of  the  wall  pressure  field  should  also  be  in  agreement. 

As  first  postulated  by  Kraichnan  (1956),  the  pressure  generated  at  the  wall  is  the  summa¬ 
tion  of  kinematic  turbulent  source  terms  throughout  the  boundary  layer  with  over  80%  of  this 
contribution  taking  place  within  the  near  wall  buffer  region  (Y+  <  50).  The  Poisson  integral 
equation  for  the  wall  pressure  represents  an  “averaging”  process  which  tends  to  smooth  out 
small  discrepancies  not  exactly  matched  in  the  velocity  field.  Therefore,  if  the  spectral  content 
of  the  wall  pressure  field  (including  convection  velocity,  scales  of  motion  and  decay  rates)  does 
not  match  between  the  two  comparative  experiments,  it  could  be  argued  that  the  numerical 
simulation  does  not  represent  a  generic  flow  in  entirety.  Whether  the  fault  lies  with  assumed 
initial  and/or  boundary  conditions,  or  that  the  simulation  is  not  considered  “mature”  enough 
in  spatial  or  temporal  space,  are  issues  to  be  resolved. 

Secondly,  because  of  the  computer  limitations  due  to  spatial  resolution  and  complex  bound¬ 
ary  conditions,  computational  experiments  are  limited  at  this  time  to  internal  (duct  or  channel) 
flows  at  very  low  Reynolds  numbers.  As  mentioned  in  chapter  1,  virtually  all  past  investiga¬ 
tions  of  the  wall  pressure  field  have  been  conducted  in  external  flow  environments  simulating 
flat  plate  boundary  layers.  The  few  studies  done  within  internal  configurations,  Clinch  (1969), 
Wambsganss  and  Zaleski  (1970),  Wilson  et  al. (1979),  and  Horne  and  Hansen  (1981),  have  been 
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confined  to  pipe  or  annulus  geometries  at  Reynolds  numbers  2  to  3  orders  of  magnitude  larger 
than  the  computational  experiment.  Therefore,  the  generation  of  a  detailed  data  base  in  a 
physical  configuration  specifically  tailored  to  the  computational  scheme  provides  for  the  first 
time,  an  informational  basis  for  comparison.  In  addition,  the  results  presented  also  show  dis¬ 
tinct  features  of  duct  flows  such  as  length  scales,  time  scales  and  decay  rates  which  differentiate 
them  from  external  boundary  layers. 

The  wall  pressure  signature  was  recorded  for  the  physical  experiment  at  five  locations  on 
the  bottom  wall  of  the  channel  for  Reynolds  numbers,  =  10000,15000,20000  and  25000. 
Transducers  2  to  5  were  located  in-line  with  the  flow  direction,  approximately  5.1cm  off  the 
channel  centerline.  Transducer  1  was  located  5.1cm  on  the  opposite  side  of  the  channel  cen¬ 
terline  and  aligned  with  transducer  2  in  the  same  cross-sectional  plane.  The  distance  between 
transducers  1  and  2  resulted  in  a  transverse  separation  of  £z/h  =  4  which  was  considered 
sufficient  to  ensure  zero  correlation  of  turbulent  energy  between  these  transducers  and  hence 
transducer  1  was  used  for  noise  cancellation  as  outlined  in  chapter  three.  The  separation  dis¬ 
tances  in  the  flow  direction  for  transducers  2  thru  5  were  in  multiples  of  (x/h  =  0.16,  with  the 
largest  distance  possible  at  £x//i  =  0.48.  During  each  run,  all  five  transducers  were  recorded 
simultaneously  together  with  the  pressure  gradient,  and  digitally  stored  on  magnetic  tape.  This 
data  was  then  transferred  to  a  VAX  780  computer  to  be  processed  by  software  developed  for 
data  analysis. 

The  digital/software  approach  for  data  analysis,  in  contrast  to  utilizing  a  spectral  analyzer, 
was  superior  for  a  number  of  reasons.  First,  for  comparison  purposes,  both  experiments  could 
be  analyzed  via  the  same  method.  Secondly,  a  number  of  combinations  of  transform  size  in 
Fourier  space  could  be  chosen  when  analyzing  the  physical  data,  in  order  to  optimize  the  best 
choice  combination  between  spectral  accuracy  and  resolution.  According  to  Bendat  and  Piersol 
(1980),  and  outlined  in  chapter  3,  the  accuracy  in  measuring  all  spectral  statistical  functions  is 
inversely  proportional  to  N,  the  number  of  ensemble  averages  obtained.  However,  due  to  the 
limited  number  of  points  recorded  in  time,  higher  resolution  in  spectral  space  requires  a  larger 
transform  size  (larger  set  of  time  points),  resulting  in  an  overall  lower  number  of  ensemble 
averages.  Performing  this  analysis  on  a  computer  with  almost  immediate  graphical  output 
provides  the  researcher  with  enormous  convenience.  A  third  reason  for  this  approach  was  also 
outlined  in  chapter  3  and  is  related  to  cancellation.  The  largest  potential  source  of  error  in 
cancellation  schemes  involving  additions  and  subtractions  in  either  “real"  or  “spectral”  6pace,  is 
that  representative  voltages  are  used  rather  than  the  actual  pressure  signals.  Hence,  errors  are 
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introduced  proportional  to  the  sensitivity  ratios  of  the  channels  being  analyzed.  In  a  realistic 
sense,  it  is  virtually  impossible  to  apply  the  actual  sensitivity  constants  to  each  channel  being 
investigated  within  the  laboratory  environment.  On  the  other  hand,  having  the  entire  data 
sets  stored  on  a  computer  where  the  analysis  is  going  to  be  performed,  provides  a  convenient 
means  for  dealing  with  the  actual  pressure  signals  by  applying  the  appropriate  sensitivities  to 
each  channel  as  the  data  is  being  read  back  into  the  analysis  program. 

Table  5.1  lists  the  various  combinations  of  record  size  and  resulting  ensemble  averages  that 
were  investigated  for  the  physical  experiment  together  with  that  possible  for  the  computational 
experiment.  For  instance,  6et  1  indicates  512  time  points  per  record  size,  with  a  resolution  of 
7.81  Hz.  The  total  number  of  points  recorded  in  time  during  am  experimental  run  at  any 
Reynolds  number,  divided  by  the  512  points  required  per  record,  resulted  in  3500  ensemble 
averages.  Sets  1  thru  6  represent  the  possibilities  investigated  for  all  Reynolds  numbers. 


Table  5.1 

Possible  Record  Size  and  Ensemble  Averages 

Set 

Aj, me 

A  frequency 

A average 

Resolution 

(time  pts) 

(frequency  pts) 

(ensemble  averages) 

(Hz) 

1 

512 

257 

3500 

7.81 

2 

1024 

513 

1750 

3.91 

3 

2048 

1025 

875 

1.95 

4 

4096 

2049 

450 

0.98 

5 

16384 

8193 

110 

0.24 

6* 

( . ) 

1061 

110/875 

C. 24/1. 95 

Numerical 

512 

257 

128/2048 

~  7.81 

*  Set  #  6 

was  a  combination  of  sets  3  and  5 

5.2  Correction  Analysis  and  rms  Pressure  Levels 

5.2.1  Auto-Spectrum  Levels  and  Magnitude  Squared  Coherence 

Representative  auto-spectra  from  transducer  2  for  all  four  Reynolds  numbers  in  the  phys¬ 
ical  experiment  are  presented  in  figures  5.1  thru  5.4.  These  figures  indicate  the  direct  output 
(uncorrected)  from  the  transducer,  the  corrected  version  using  the  technique  outlined  in  chap¬ 
ter  3,  and  the  background  noise  level  with  all  turbulence  removed.  These  results  were  taken 
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from  data  set  2  in  table  5.1  because  the  number  of  time  points  and  ensemble  averages  gave  an 
expected  accuracy  in  the  auto-spectrum  on  the  order  of  2.4  %  and  the  resolution  obtained  was 
sufficient  to  show  effective  cancellation.  A  number  of  interesting  features  are  indicated.  All 
four  Reynolds  numbers  show  that  although  the  channel  attempts  to  minimize  acoustic  and/or 
background  facility  generated  noise  by  operating  in  a  “blowdown”  mode,  a  very  intense  con¬ 
taminating  noise  level  exists  for  frequencies  less  than  approximately  10#  z.  An  analysis  was 
conducted  to  determine  if  organ  pipe  resonances  were  present  in  the  system,  but  the  frequencies 
computed  were  at  least  an  order  of  magnitude  larger  than  calculated. 

In  order  to  better  illuminate  the  lower  frequency  region,  separate  runs  were  made  at  each 
Reynolds  number  over  a  bandwidth  of  100#z,  with  a  resolution  of  approximately  0.1  Hz.  The 
results  showed  three  distinct  frequencies  located  at  approximately  1,  5  and  10#z  respectively, 
with  their  intensities  decreasing  for  each  successively  larger  frequency  component.  Despite  ex¬ 
tensive  analysis,  no  sources  for  these  contaminating  frequencies  could  be  isolated.  Accelerome¬ 
ters  were  placed  on  the  test  section  at  a  number  of  locations  to  determine  if  structure  vibration 
was  being  transmitted  to  the  fluid  within  the  channel.  There  was  virtually  zero  coherence  with 
the  accelerometer  output  and  that  of  the  pressure  transducers,  and  no  peaks  were  indicated  at 
these  three  distinct  frequencies.  It  was  surmised  that  the  large  upstream  reservoir  might  be  the 
source.  In  order  to  maintain  constant  velocity  during  each  run,  high  pressure  air  is  discharged 
in  the  top  of  the  tank  to  ensure  a  constant  overall  dynamic  pressure.  One  possible  explanation 
is  that  a  very  low  frequency,  slow  moving  pressure  bubble  is  created  in  the  upstream  tank,  due 
to  the  interaction  of  the  high  pressure  air  with  the  free  surface  of  the  contained  water.  This 
low  frequency  bulk  wave  is  then  transmitted  acoustically  to  the  channel  test  section  within  the 
flowing  fluid. 

Another  interesting  feature  indicated  in  these  figures  is  the  limited  dynamic  range  of  the 
transducers  used.  Figure  5.1  6hows  that  at  a  Reynolds  number  of  10,000,  the  high  frequency 
end  of  the  spectrum  is  dominated  by  electronic  noise,  and  the  correction  scheme  tails  to  ade¬ 
quately  cancel  the  three  distinct  contaminating  frequencies  at  the  low  end.  It  should  be  noted 
however,  that  further  examination  of  figure  5.1  shows  that  cancellation  works  rather  well  at 
high  frequencies  where  turbulence  is  present,  but  lower  in  intensity  than  the  electronic  noise 
evident.  This  supports  the  hypothesis  that  the  “Wiener”  cancellation  scheme  presented  in 
chapter  3  is  not  only  “not”  bandwidth  limited  but  that  it  produces  excellent  results,  provided 
the  level  of  the  contaminating  noise  is  within  20 dB  of  the  desired  signal.  If  one  accepts  the 
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hypotheses  of  Panton  (1980)  and  Chase  (1980),  that  the  low  frequency  end  of  the  wall  pres¬ 
sure  signature  should  drop  off  as  u>2,  then  the  ratio  of  amplitudes  over  a  decade  in  frequency 
would  result  in  an  approximate  drop  of  20 dB  in  amplitude.  Figure  5.1  has  a  maximum  around 
—  10 dB  at  around  10#*  and  therefore  should  drop  to  about  -30 dB  at  approximately  1#*. 
The  intensity  of  the  contaminating  noise  is  on  the  order  of  50  to  60 dB  higher  in  this  range  and 
therefore  requires  extreme  cancellation.  The  lack  of  effectiveness  at  this  low’  Reynolds  number 
is  not  considered  to  be  a  failure  of  the  correction  scheme.  In  chapter  3  the  assumptions  of 
“homogeneity”  and  “signal-free"  reference  were  critical  to  obtaining  accurate  cancellation.  For 
very  low  Reynolds  numbers,  the  length  scales  and  decay  rates  of  turbulence,  particularly  the 
large  scale  components  at  low  frequency,  might  be  large  and  slow  enough  to  extend  transversely 
past  the  limited  separation  of  transducers  1  and  2.  Hence,  part  of  the  turbulence  itself  is  be¬ 
ing  cancelled  or  otherwise  modified.  Examining  the  other  three  Reynolds  numbers  shows  that 
cancellation  is  adequate  and  the  overall  background  noise  level  in  the  channel  is  approximately 
10  to  20 dB  below  the  turbulent  energy  spectrum  levels. 

In  order  to  show  the  full  range  of  the  cancellation  scheme  and  to  better  illustrate  the  three 
contaminating  frequencies,  the  composite  data  set  6  from  table  5.1  is  shown  in  figure  5.5  for  all 
four  Reynolds  numbers.  This  data  set  combines  two  analyses,  one  with  a  resolution  of  1.95 Hz 
and  one  with  a  resolution  of  0.244#*  with  875  and  110  ensemble  averages  respectively.  The 
merging  of  these  two  sets  is  apparent  in  figure  5.5  by  the  obvious  waviness  starting  at  10#* 
and  lower.  The  most  interesting  information  depicted  in  figure  5.5.  is  the  apparent  universal 
drop-off  in  spectrum  level  proportional  to  frequency  with  a  slope  of  1.7  (/1-7).  Although  an  w2 
dependence  resulting  in  2 OdB  difference  in  amplitudes  is  not  exactly  matched,  it  is  more  likely 
a  result  of  the  lower  spectral  accuracy,  due  to  an  ensemble  average  that  at  most,  only  gives 
an  expected  error  good  to  approximately  10  percent.  In  either  case,  figure  5.5  presents  “new” 
information  for  internal  flow  geometries  and  validates  past  and  present  theoretical  suppositions 
for  similar  frequency  dependence  for  external  flows. 

Finally,  figures  5.6  thru  5.9  show  representative  results  for  the  magnitude  squared  coher¬ 
ence  7J,  for  the  combinations  of  transducers  1 :  2, 2  :  3, 2  :  4  and  2  :  5  for  all  Reynolds  numbers. 
These  figures  are  obtained  from  data  set  2  in  table  5.1  for  illustrative  purposes  and  serve  to 
indicate  two  important  points.  First,  the  decay  in  intensity  levels  as  well  as  loss  in  coherence 
with  frequency  at  increasing  separation  for  transducers  2  thru  5  in  the  flow  direction  is  distinct 
for  all  Reynolds  numbers.  Secondly,  the  solid  line  represents  the  coherence  between  transducer 
1,  which  was  used  as  a  cancellation  source,  and  transducer  2  which  was  the  lead  transducer  in 
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the  linear  array.  The  separation  of  these  two  transducers  by  £2/h  =  4  was  sufficient  to  ensure 
zero  correlation  of  turbulence  in  a  cross-sectional  plane  for  the  largest  three  Reynolds  numbers. 
Examination  of  ail  four  figures  demonstrates  that  in  reality,  correction  need  only  be  applied 
for  frequencies  below  approximately  10 Hz.  In  addition,  as  pointed  out  in  Horne  and  Handler 
(1989),  the  potential  error  resulting  from  a  non-signal  free  reference  is  specifically  determined 
by  the  noise  coherence  level.  The  maximum  error  possible  is  equivalent  to  3 dB,  when  the 
magnitude  squared  coherence  is  equal  to  unity.  This  is  certainly  indicated  in  figures  5.6  to  5.9. 
within  experimental  and  computational  error. 


5.2.2  Root-Mean-Square  Pressure  Levels  (rmj) 

As  pointed  out  in  chapter  3,  there  are  three  methods  of  computing  the  rms  level  of  the 
wall  pressure  depending  on  the  method  of  analysis  chosen.  The  primary  method  is  to  utilize 
equation  (2.18)  and  average  the  time  history  data  directly,  provided  there  is  no  low  frequency 
contaminating  noise.  The  second  and  third  methods  depend  on  whether  the  correlation  co- 
variance  is  solved  for,  as  represented  by  equation  (2.32),  or  whether  the  power  spectral  density 
is  computed  directly  by  Fourier  transforming  the  pressure  signals  and  integrating  the  auto- 
spectrum  via  equation  (2.36).  All  three  methods  must  give  the  same  results,  and  in  reality,  the 
last  two  are  identical  as  the  correlation  covariance  and  spectral  density  functions  are  Fourier 
transform  partners. 


For  the  physical  experiment,  data  set  6  was  analyzed  in  order  to  take  advantage  of  the 
cancellation  scheme  to  obtain  an  accurate  representation  of  a  generic  channel  flow  as  is  possible 
to  date.  As  shown  in  the  previous  figures  5.1  to  5.9  it  was  not  possible  to  average  the  time 
history  records  to  obtain  prm *,  due  to  the  high  intensity  low  frequency  noise  present  in  each 
signal.  Therefore,  the  corrected  spectra  were  integrated  throughout  the  frequency  range  using 
equation  (2.36),  but  modified  for  frequency,  as  shown  by, 


P 


2 

rmi 


(5.1) 


(Note  that  $pp(u>)  differs  from  $pp(/)  by  the  factor  4 jt).  As  a  check  of  the  method  of  analysis, 
the  raw,  uncorrected  data  time  histories  were  averaged  and  compared  with  the  integrations  of 
the  uncorrected  auto-spectrums  for  which  identical  results  were  obtained.  The  rms  levels  are 
presented  in  table  5.2  as  ratios  Prmi/Q  and  pTmtlTwan  together  with  other  related  character¬ 
istics  describing  all  flow  configurations  investigated. 
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Table  5.2  Parameter?  Governing  Both  Experiments 


Type 

V avg 

Prm  1 

Tws|{ 

Pr  mi 

Q 

ti" 

d+ 

R" 

Physical 

2.46 

2.84 

0.0174 

0.056 

23 

572 

3.87 

2.12 

0.0106 

0.050 

32 

806 

4.80 

2.36 

0.0108 

0.048 

38 

955 

6.18 

2.26 

0.0092 

0.045 

46 

1156 

Numerical 

1.32 

0.0092 

0.059 

~  40 

125 

It  should  be  noted,  as  pointed  out  in  chapters  3  and  4  that  the  physical  experiment 
was  really  a  composite  of  two  experiments;  one  in  which  the  kinematic  (velocity)  data  was 
documented  and  the  second  for  which  the  wall  pressure  results  were  obtained.  In  the  wall 
pressure  experiment,  no  values  for  the  centerline  velocity  were  possible  as  the  value  for  Vavg 
was  the  only  variable  recorded.  Therefore,  the  average  velocity  was  used  to  normalize  the 
rms  pressure  by  calculating  a  dynamic  pressure  term  as  represented  by  Q  —  \pl7lvg.  Because 
there  is  no  ratio  of  Uavg/Uoc  generally  accepted  within  the  literature  for  channel  flows,  actual 
measurable  values  were  used  to  determine  these  ratios. 

Figure  5.10  6hows  the  values  of  prmi/Q  taken  from  table  5.2,  plotted  versus  transducer 
diameter  normalized  by  the  viscous  length  scale  vju' .  Included  are  the  measurements  of 
Blake  (1970),  Schewe  (1983),  and  Farabee  (1986).  It  is  not  clear  what  value  to  use  for  d+  for 
the  computational  results,  even  though  the  value  of  d+  ~  40  is  indicated  in  table  5.2.  The 
separation  distance  between  computational  grid  points  in  the  flow  direction  is  £xu* /v  =  40.  and 
limits  the  smallest  scales  resolvable.  As  shown  in  figure  5.10,  both  experiments  agree  quite  well 
with  the  excellent  data  of  Blake  and  Schewe,  and  the  more  recent  data  of  Farabee.  It  should 
be  noted  however,  that  Blake's,  Schewe’s  and  Farabee’s  data  were  obtained  from  experiments 
conducted  within  flat  plate  boundary  layers. 


5.3  Spectral  Features  of  Channel  Flow  Wall  Pressure 
5.3.1  Single  Point  Frequency  Spectra 

The  single  point  frequency  spectra  were  measured  and  noise  cancellation  techniques  applied 
at  all  five  locations  in  the  channel  test  section.  Except  for  the  lowest  Reynolds  number  at 
10, 000,  all  other  results  demonstrated  turbulent  homogeneity  and  produced  similar  results  at  all 
Reynolds  numbers.  The  single  point  frequency  spectra  at  the  Reynolds  numbers  Rh  =  15,000 
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and  25,000  were  chosen  to  compare  with  both  the  computational  experiment  and  flat  plate 
data  from  the  literature.  These  Reynolds  numbers  were  approximately  7  and  12  times  larger 
than  the  computational  experiment  respectively,  and  table  5.2  indicates  their  values  for  R‘  as 
806  and  1156.  Figures  5.2  and  5.4  show  that  the  intensity  level  of  the  pressure  fluctuations 
was  sufficiently  above  the  background  noise  and  the  corrected  results  show  good  cancellation 
of  the  low  frequency  acoustic  contamination. 

The  velocity  fluctuations  throughout  the  boundary  layer  comprise  the  source  terms  for 
the  wall  pressure  fluctuations.  If  one  examines  the  length  and  time  scales  which  appear  to 
collapse  the  pressure  spectra  to  those  scales  characteristic  of  the  turbulent  boundary  layer,  it 
may  be  possible  to  determine  the  specific  origins  of  these  source  terms.  Detailed  information 
which  quantify  these  regions  more  explicitly  must  be  obtained  from  cross-spectral  information. 
However,  at  least  in  a  global  sense,  single  point  data  may  be  examined  to  determine  orders  of 
magnitude. 


Figures  5.11  thru  5.16  show  the  single  point  pressure  spectra  normalized  by  three  sets 
of  variables  commonly  used  in  the  literature.  For  these  and  subsequent  plotted  results,  the 
following  relations  were  used. 


Uavg/U  oo  —  0.84 

(5-2) 

* 

II 

^7 

(5.3) 

6/h  =  1 

(5.4) 

physical  hydraulic  diameter  ,  =  1.89 

2  h 

(5.5) 

computational  hydraulic  diameter  ,  ^  =  2.00 

(5.6) 

It  must  be  noted,  that  the  above  relations  are  used  only  to  compare  the  experimental  results 
with  typical  data  in  the  literature  which  represent  flat  plate  boundary  layers  and  for  which 
these  parameters  are  usually  measured  and  have  a  certain  significance.  Although  there  has 
been  no  acceptable  value  for  Uavg/U »  reported  in  the  literature,  examination  of  the  velocity- 
data  for  the  physical  experiment  presented  in  chapter  3  indicates  that  the  above  value  given 
by  equation  (5.2)  is  reasonable.  Figure  5.11  shows  the  frequency  spectra  scaled  on  the  outer 
flow  length  and  velocity  scales  6m  and  Uoo  respectively,  in  the  form, 


— c -  vs,  - . 

Q2f  i/= o 


(5.7) 
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This  is  usually  termed  “outer-variable”  scaling  and  uses  a  time  scale  S' /Ux  associated  with 
the  outer  part  of  the  boundary  layer.  Examination  of  the  results  from  the  physical  experiment 
reconfirm  the  u2  dependence  at  low  frequencies  mentioned  in  section  5.2.1.  The  time  scale 
could  be  thought  of  as  a  combination  of  a  velocity  scale  V and  a  length  scale  S'  (or  6), 
which  are  characteristic  of  the  outer  flow,  large  scale,  slower  moving  structure  (low  frequency). 
Despite  both  the  physical  and  numerical  data  indicating  an  J1  dependence,  the  normalization 
used  doesnot  collapse  both  experimental  sets  of  data  exactly  for  very  low  frequencies.  In  a 
channel,  the  outer  flow  consists  of  that  structure  which  exists  and/or  occupies  the  entire  channel 
cross-section,  hence  the  centerline  velocity  and  channel  half-height  Uoo  and  h  respectively, 
certainly  constitute  proper  velocity  and  length  scales  for  normalization  of  the  pressure  spectra 
amplitudes.  For  the  flat  plate  flows  of  Blake  (1970)  and  Farabee  (1986),  Uoo  and  S’  represent 
proper  outer  flow  variables,  but  again,  no  collapse  of  the  data  is  readily  apparent  between  the 
two  internal/external  flow  regimes. 


Figure  5.12  shows  the  frequency  spectra  scaled  on  a  combination  of  inner  and  outer  flow- 
variables  in  the  form. 


vs. 


1  wall* 


U )S‘ 

uZ 


(5.8) 


A  time  scab  associated  with  the  outer  flow  is  used,  but  the  wall  shear  stress  is  now  incorporated 
to  normalize  the  intensity  of  the  wall  pressure  fluctuations.  Corcos  (1963),  Blake  (1970)  and 
Farabee  (1986)  have  all  shown  that  the  maximum  in  the  wall  pressure  spectrum  should  occur 
at  u>Sm /Ua o  approximately  equal  to  0.3.  The  flat  plate  boundary  layer  data  indicated  in  figure 
5.12  substantiate  this  value,  but  the  channel  flow  data  from  both  the  physical  and  numerical 
experiments  do  not  agree  with  this  result.  Figure  5.13  modifies  the  data  in  figure  5.12  by 
computing  a  hydraulic  diameter  given  by, 


H ydraulic  diameter  ,  d*  =  — — : - — .  (5.9) 

Perimeter  ’ 

The  value  of  d^/h  equals  4.0  and  3.8  for  the  computational  and  physical  experiments  re¬ 
spectively.  As  shown  in  figure  5.13,  the  maximum  in  the  channel  flow  data  does  peak  at 
uS’/U oo  =  0.3.  However  there  is  a  marked  difference  between  internal  and  external  flows  as 
indicated  by  the  flat  plate  data  of  Blake  (1970)  and  Farabee  (1986).  If  one  were  to  use  the 
hydraulic  diameter  only  on  the  x-axis  in  figure  5.13,  the  agreement  at  uS’/Uoo  =  0.3  is  met  as 
shown  in  figure  5.14,  but  major  differences  still  exist  within  the  low  frequency  region. 

A  number  of  points  can  be  made  from  the  collective  examination  of  figures  5.11  thru  5.14. 
First,  scaling  only  on  outer  flow  variables  does  appear  to  collapse  the  data  between  the  two 
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flow  configurations,  interna]  and  externa]  provided  hydraulic  diameter  is  used.  This  is  some¬ 
what  subjective,  as  there  is  virtually  no  information  in  the  literature  at  low  frequency.  Second, 
normalization  of  the  wall  pressure  fluctuations  with  Twan  in  mixed  (inner-outer)  variable  scal¬ 
ing  does  collapse  the  spectrum  amplitudes  for  internal  and  external  flows  separately,  and  the 
maximum  in  the  spectrums  does  appear  to  be  universal  at  uS'/Uoc  —  0-3,  provided  internal 
flows  can  be  characterized  by  hydraulic  diameter.  Third,  the  most  interesting  feature  with 
mixed  variable  scaling  is  the  important  role  that  the  wall  shear  plays  in  establishing  the  mag¬ 
nitudes  of  the  low*  frequency  pressure  fluctuations.  This  implies  that  inner  flow  length  scales  are 
much  different  between  the  two  regimes.  Fourth,  the  computational  experiment  demonstrates 
marked  degradation  at  low  frequencies.  The  u>2  dependence  appears  to  be  indicated,  but  the 
breadth  of  the  spectrum  at  low  frequencies  is  not  of  sufficient  bandwidth  to  collapse  this  data 
with  the  physical  experiment.  This  is  a  result  of  insufficient  time  history  data  and  the  periodic 
boundary  conditions  interfering  with  the  flow.  Finally,  neither  of  the  two  sets  of  variable  nor¬ 
malization  hold  for  the  high  frequency  portion  of  the  spectra.  This  is  not  unexpected  because 
smaller  scale  structures  corresponding  to  higher  frequencies  exist  within  the  inner  layer  of  the 
boundary  layer  and  have  a  time  scale  of  j//u*2.  Collectively,  these  results  show  that  external 
and  internal  flow  configurations  constitute  two  physically  different  flow  regimes  whose  outer 
and  inner  flows  are  “not"  similar. 


Figure  5.15  scales  the  high  frequency  portion  of  the  wall  pressure  fluctuations  on  inner 
variables.  The  velocity  scale  is  u*  and  the  length  scale  is  represented  by  v/uu.  Both  the 
physical  and  computational  data  sets  collapse  in  the  region  0.4  <  wi//u*2  <  1.  This  region 
is  represented  by  an  w-7/3  dependence,  as  shown  by  Schewe  (1984).  There  is  a  discrepancy 
in  the  computational  data  at  ui//um2  ~  1.  The  shape  of  the  numerical  spectrum  appears  to 
change  slope  sharply  beyond  this  frequency,  which  is  typical  of  reduced  resolution  due  to  finite 
transducer  size  in  physical  experiments.  One  of  the  critical  aspects  of  DNS  calculations  as 
6hown  by  Reynolds  (1989),  is  that  the  grid  resolution  size  must  be  small  enough  to  resolve 
both  the  largest  and  smallest  eddies  that  are  dynamically  significant  to  the  flow.  The  grid 
6ize  for  the  computational  experiment  is  represented  by  z*  =  £ru*/*/  =  40.  From  table  5.2, 
u'/t/oo  =  .059.  As  this  region  of  the  boundary  layer  is  very  dose  to  the  wall,  an  assumption  of 
convection  velodty  Ue/Uoo  =  0.5  is  not  unreasonable.  Using  these  values,  a  cut-off  frequency 
due  to  finite  ‘computational’  resolution  could  be  determined  as, 


f  uv  i  _  E  2?r 
lu.2  leutojf  ~  •  x.  u' I U,*, 


1.33. 


(5.10) 
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Apparently,  the  discrepancy  between  the  physical  and  computational  experiments  at  very  high 
frequencies,  is  a  result  of  the  relatively  poor  resolution  capability  resulting  from  the  computa¬ 
tional  grid  size. 

Figure  5.16  shows  the  computational  spectrum  compared  with  the  physical  experiment 
Reynolds  number  Rh  =  15,000  from  the  composite  data  set  6  and  set  1  which  also  had  512 
time  points  per  ensemble  record.  Clearly,  the  only  discrepancies  in  the  software  approach  for 
the  physical  data  is  in  that  region  where  the  higher  resolution  had  only  100  averages.  This  is 
common  in  power  spectra  that  do  not  contain  enough  time  points  to  achieve  sufficient  accuracy. 
The  amplitudes  exhibit  a  waviness  and  have  not  settled  down  to  their  steady  state  values.  As 
a  final  note,  the  wall  pressure  spectrum  from  the  computational  experiment  exhibits  a  number 
of  peaks  in  the  low  frequency  region.  Similar  to  the  physical  experiment,  these  peaks  represent 
noise  and  are  a  result  of  the  periodic  boundary  conditions  in  the  DNS  calculations.  This 
phenomenon  will  be  addressed  further  in  section  5.4. 


5.3.2  Cross- Spectra]  (two-point)  Pressure  Measurements 

Measurements  of  the  two-point  space-time  statistics  of  the  wall  pressure  field  aid  in  quan¬ 
tifying  the  convective  properties  more  fully  and  provide  insight  into  the  size  of  the  eddies 
representative  of  the  inner  and  outer  flow.  The  magnitude  of  the  cross-spectrum  is  expressed 
in  normalized  form  as  the  square- root  of  the  magnitude-squared  coherence  function,  as  out¬ 
lined  in  chapter  2,  equation  (2.37).  For  each  successive  streamwise  separation,  the  phase  of  the 
cross-spectrum  represents  a  narrow-band  interpretation  of  the  convection  velocity  and  has  the 
form, 


Uc  = 


o(w.€r)' 


(5.11) 


For  flow  in  the  streamwise  direction,  Taylor  proposed  that  the  cross  spectrum  phase  can 
assume  a  physical  meaning  when  expressed  as, 
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(5.12a) 


or 


—  u^x  —  u  c 

■a  ~  V'  -  Mr. 


(5.125) 
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where  kx  =  2x/Xx.  As  shown,  the  phase  can  be  interpreted  as  a  non-dimensional  wave  number. 
By  plotting  the  magnitude  of  the  coherence  versus  the  phase,  an  indication  of  the  decay  of  the 
dominant  pressure  producing  eddies  of  a  specific  size  Xx  may  be  obtained. 

Two-point  measurements  were  made  for  the  physical  experiment  at  three  streamwise  sep¬ 
arations  characterized  by  £x/h  =  0.16,  0.32  and  0.48.  The  computational  experiment’s  grid 

spacing  provided  ratios  of  (£x//i)j  =  j(0.32)  with  j  =  1,2, . 15.  Only  128  ensemble  averages 

were  possible  for  the  computational  experiment  when  calculating  the  cross-spectrum.  In  addi¬ 
tion,  the  short  length  of  the  computational  box  resulted  in  the  periodic  boundary  conditions 
affecting  adversely  the  phase  information,  therefore,  only  the  physical  experimental  data  for 
coherence  will  be  displayed.  Figures  5.17  to  5.19  present  the  streamwise  coherence  plotted 
versus  the  phase  a(u>,(x)  for  the  largest  three  Reynolds  numbers.  Bull  (1967)  proposed  that 
the  streamwise  coherence  exhibited  an  exponential  decay  of  the  form, 

I  7(«,&)  I  =  exp  -C0(2-£)  ,  (5.13) 

and  this  relation  is  indicated  in  each  figure.  Bull  found  a  value  for  C0  equal  to  approximately 
0.12.  For  other  flat  plate  flows  in  the  literature,  C0  has  taken  on  values  from  0.1  to  0.2  and 
appear  to  decrease  as  the  Reynolds  number  increases.  The  more  recent  study  of  Farabee  (1986) 
found  values  of  C0  =  0.145  and  0.125.  In  figures  5.17  to  5.19,  the  values  indicated  are  larger, 
ranging  from  0.27  to  0.18  and  do  exhibit  this  same  Reynolds  number  dependence.  Bull  found 
that  the  streamwise  coherence  decayed  to  a  value  of  0.05  when  a(u!,{x)  equaled  24.5.  Using 
equations  (5.12a)  and  (5.13),  he  obtained  the  relation, 

ix  ~  4AX  .  (5.14) 

This  suggests  that  the  largest  eddies  of  size  Xx  lose  coherence  upon  being  convected  4  times 
their  size.  Similar  results  were  found  by  Willmarth  and  Wooldridge  (1962)  and  Farabee  (1986). 
However,  the  data  of  the  physical  experiment  presented  here  show  that  eddies  of  a  specific  size 
Xx  lose  all  their  coherence,  upon  being  convected  a  distance  2  to  2|  times  their  respective  sizes. 
For  example,  taking  the  value  C0  =  0.22  from  figure  5.18  and  evaluating  the  term  in  brackets 
from  equation  (5.13)  where  |  7(w,(z)  j=  0.05  results  in, 

tx  =  2.2  Ar  .  (5.15) 

This  implies  a  faster  decay  rate  for  internal  flow  geometries  versus  external  flat  plate  boundary 
layers  which  have  no  constraining  wall  to  limit  the  size  of  the  boundary  layer  thickness  6.  This 
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suggests  that  channel  flows  which  are  characterized  by  a  fixed  length  scale  (either  2/i  for  ducts, 
or  D  the  diameter  in  pipes)  force  eddies  of  a  specific  size  to  decay  in  half  the  distance  for 
internal  flow's  versus  external.  It  is  not  clear  whether  this  has  been  previously  observed  since 
the  pipe  flow  experiments  of  Bakewell  (1963)  and  Clinch  (1969)  were  broadband  analyses  for 
which  transducer  resolution  was  affected  by  cancellation  effects. 

In  duct  flows,  to  a  large  extent,  it  has  been  proposed  that  turbulence  is  characterized 
primarily  by  what  goes  on  close  to  the  wall  rather  than  in  regions  closer  to  the  centerline.  The 
outer  intermittent  region  of  external  flows  represented  by  boundary  layers  have  been  hypothe¬ 
sized  to  be  different  than  the  region  existing  at  the  channel  centerlines.  This  is  substantiated 
somewhat  by  the  kinematic  results  from  chapter  4,  in  that  the  viscous  sub-layer,  buffer,  and 
log-law  regions  are  very  similar  to  those  results  reported  for  external  flows.  However,  the 
shorter  decay  lengths  and  lower  intensity  values  for  all  frequencies  in  the  pressure  spectra  for 
channel  flows  gives  evidence  that  both  internal  and  external  boundary  layers  are  fundamentally 
different. 


5.3.3  Convection  Velocity  (narrow-band) 

A  narrow  band  analysis  of  the  convection  velocity  is  performed  by  solving  equation  5.11 
for  each  separation  distance  between  transducers.  For  both  experiments,  these  results  are 
presented  in  figures  5.20  to  5.23  scaled  on  outer  variables  h  and  V «>,  and  in  figures  5.24  to  5.27 
scaled  on  the  inner  variable  time  scale  v/u'2. 

The  most  distinguishing  feature  displayed  by  these  results  is  the  Reynolds  number  depen¬ 
dence.  First,  at  all  Reynolds  numbers  (including  the  computational  experiment),  the  convection 
velocity  decreases  to  zero  with  decreasing  frequency.  This  is  a  trend  first  measured  by  Blake 
(1970)  and  Wills  (1970)  and  predicted  by  Chase  (1980)  and  recently  confirmed  by  Farabee 
(1986).  Because  this  trend  is  indicated  no  matter  the  size  of  transducer  separation,  the  im¬ 
plication  is  that  a  finite  lower  frequency  limit  exists  for  which  there  is  measurable  coherence. 
Referring  to  figures  5.17  to  5.19  this  is  shown  dearly  by  the  drop  in  coherence  at  low  frequency. 
A6  Farabee  showed,  equation  (5.13)  can  achieve  unity,  either  by  u>  going  to  zero  or  going  to 
zero.  The  coherence  going  to  1,  as  goes  to  zero,  is  a  required  result,  since  the  cross-spectrum 
essentially  becomes  the  auto-spectrum  for  which  7  =  1  for  all  frequencies.  However,  if  7  goes 
to  1  as  w  goes  to  zero,  this  implies  that  low  frequency  components  remain  coherent  over  very 
long  distances,  without  being  distorted  by  the  flow,  which  is  not  physically  realizable.  Thus, 
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as  figures  5.17  to  5.19  show,  there  is  a  loss  of  coherence  with  decreasing  frequency,  and  as  the 
convection  velocities  support,  there  is  a  finite  frequency  for  which  structure  at  all  wavenumbers 
exist  within  the  channel  and  are  convected  downstream.  Secondly,  for  all  Reynolds  numbers, 
the  convection  velocity  assumes  higher  values  for  larger  separation  distances,  implying  that  as 
the  transducer  separation  increases,  smaller  scale  structure  has  decayed  and  all  that  remains 
is  the  larger  size  eddies  which  physically  must  exist  further  out  into  the  flow  and  are  subject 
to  a  higher  velocity.  Thirdly,  at  a  fixed  frequency,  say  uh/Uoo  =  0.3,  as  the  Reynolds  number 
increases,  even  smaller  scale  structure  exists  within  the  flow  and  is  convected  downstream.  This 
i6  sensed  by  the  transducers,  by  registering  an  even  higher  convection  velocity,  from  a  value  of 
UJUoc  equal  to  0.65  for  the  low  Reynolds  number,  to  a  value  of  0.98  at  a  Reynolds  number 
of  25, 000.  The  convection  velocity  ratio  at  the  smaller  separation  seems  to  asymptote  at  large 
frequencies  to  a  value  of  approximately  0.55  and  for  the  larger  separations,  a  value  closer  to 
0.75. 

The  computational  experiment  exhibits  the  same  Reynolds  number  dependence  as  shown 
by  the  physical  experiment,  except  for  the  relative  peak  values  and  magnitude  dependence 
on  separation  distance.  As  mentioned  previously,  the  cross-spectrum  phase  was  very  difficult 
to  measure,  due  to  the  low  number  of  ensemble  averages  available  and  interference  from  the 
periodic  boundary  conditions.  Therefore,  the  only  expectation  was  to  confirm  the  Reynolds 
number  dependence  from  this  set  of  data.  Figures  5.24  to  5.27  show  the  same  information  just 
discussed,  except  they  are  scaled  on  the  inner  variable  time  scale,  i //u“2. 

In  section  5.3.2,  an  assumption  was  made  that  kx  =  u/Uc(w)  which  is  Taylor’s  hypothesis 
that  a  unique  convection  velocity  exists  that  relates  the  wavenumber  of  a  certain  fluctuation 
to  its  frequency.  Figures  5.20  to  5.27  show  that  the  convection  velocity  varies  with  £r/h.  This 
supports  the  belief  that  the  pressure  field  is  dispersive,  in  which  fluctuations  of  approximately 
the  same  Bize  may  convect  at  different  velocities  depending  on  their  location  within  the  flow. 
This  idea  was  explicity  6hown  by  Wills  (1970)  who  spatially  transformed  his  cross-spectral 
measurements  into  contour  plots  showing  both  the  dispersive  nature  of  the  wall  pressure  field 
and  that  for  decreasing  frequency,  UefU^o  should  tend  to  zero. 
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5.4  Temporal  Features  of  Channel  Flow 


5.4.1  Correlation  Analysis  (broadband) 

As  6hown  in  the  previous  sections,  a  narrow-band  analysis  consists  of  evaluating  the  cross- 
spectral  density  functions  in  narrow  frequency  bandwidths.  Assuming  that  Taylor’s  hypothesis 
holds,  this  provides  a  means  for  evaluating  the  properties  of  the  wall  pressure  field  for  ail  wave 
numbers  contained  within  a  narrow  frequency  band  centered  at  u.  The  6pace-time  correlation 
function  is  the  Fourier  transform  of  the  cross-spectra  for  which  all  phase  information  is  lost. 
The  retained  information  consists  of  the  relative  measure  of  correlation  of  the  large  scale  com¬ 
ponents,  and  the  distances  and  limits  in  time  for  which  they  remain  coherent.  In  addition,  a 
broad  band  convection  speed  can  be  determined  for  all  scales  and  frequencies  by  the  proper 
interpretation  of  this  function. 

For  the  most  part,  a  broadband  anaysis  does  not  provide  very  much  detail  about  the 
turbulence  present  and/or  its  relationship  to  the  wall  pressure.  On  the  other  hand,  one  of 
the  goals  of  this  experimental  endeavor  was  to  evaluate  the  wall  pressure  results  from  a  DNS 
calculation  by  comparing  with  those  of  the  physical  experiment.  All  information  extracted 
from  both  experiments  should  be  considered  valuable.  In  addition,  one  of  the  assets  of  the 
software  approach  to  analysis,  was  to  evaluate  both  experiments  via  the  same  approach.  As 
mentioned,  the  computational  experiment  had  indications  of  low  frequency  noise  contamination 
and  the  possibility  of  resolution  problems  due  to  insufficient  grid  size.  The  temporal  features,  as 
obtained  through  the  Fourier  transform  of  the  cross- spectral  density  function,  will  substantiate 
these  observations  more  clearly. 

The  space-time  correlation  functions  for  the  physical  experiment  are  presented  in  figures 
5.28  to  5.31  for  all  four  Reynolds  numbers  as  a  function  of  r,  the  separation  in  time.  Figure 
5.32  contains  similar  results  for  the  computational  experiment  as  a  function  of  7'*,  a  non- 
dimensional  time.  As  mentioned  previously,  there  was  a  low  frequency  contamination  centered 
at  approximately  1,  5  and  10  hz,  respectively  for  all  four  Reynolds  numbers.  This  low  freqency 
noise  had  such  large  amplitudes,  that  it  completely  swamped  the  time  signals.  The  single  point 
auto-spectra  presented  in  section  5.3.1  showed  a  contamination  over  the  low  range  of  frequencies 
which  were  30  to  50 dB  higher  then  the  rest  of  the  spectrum  levels.  When  Fourier  transforming 
frequency  spectra,  large  amplitudes,  particularly  at  low  frequencies,  tend  to  mask  the  overall 
output.  This  can  be  observed  in  figures  5.28  to  5.31,  by  the  almost  solid  horizontal  lines  at  the 
top  of  each  figure.  These  lines  represent  the  Fourier  transforms  of  the  “un-corrected”  spectrums. 
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The  other  four  curves  represent  the  fourier  transforms  of  the  “corrected"  cross-spectra]  density 
functions  from  section  5.3.2.  It  was  speculated  in  section  5.2.1  that  the  separation  between 
the  two  lead  transducers  1  and  2  for  the  low  Reynolds  number  at  10,000  was  insufficient  for 
total  cancellation  of  the  contaminating  noise.  This  is  evident  in  figure  5.28  by  the  somewhat 
arbitrary  coefficients  obtained  at  each  maximum.  The  higher  Reynolds  numbers,  for  which 
cancellation  produced  excellent  results,  clearly  show  the  dominant  characteristic  of  the  wall- 
pressure  field  as  a  convecting  pattern  which  rapidly  loses  coherence.  Although  there  were  not 
sufficient  transducers  to  obtain  a  much  broader  separation  in  the  streamwise  direction,  the 
three  separations  presented  do  indicate  slight  asymmetries  for  each  successive  correlation  as 
represented  by  a  dispersive  field. 

The  results  from  the  computational  experimental  results  presented  in  figure  5.32  show 
similar  behavior  for  the  first  nine  separational  grid  points.  Asymmetry  in  each  successive  cor¬ 
relation  is  somewhat  subjective  to  isolate,  however,  this  is  consistent  with  a  very  low  Reynolds 
number  flow  which  would  only  be  weakly  dispersive.  The  most  interesting  feature  in  the  com¬ 
putational  results  is  the  indication  of  periodicities  at  large  time  separations  for  each  correlation. 
This  is  a  direct  result  of  the  periodic  boundary  conditions  used  within  DNS  calculations.  If 
the  computational  box  over  which  the  flow  is  being  computed  is  too  small  for  the  Reynolds 
number  being  investigated,  the  flow  senses  or  feels  this  periodicity.  This  same  periodicity  was 
observed  in  the  single  point  spectra  of  the  computational  results  in  figures  5.11  thru  5.16  as 
the  individual  peaks  at  low  frequency.  This  periodicity  represents  noise,  and  is  very  similar  to 
that  observed  in  the  physical  experiment.  The  only  method  of  eliminating  it,  is  to  increase  the 
box  size  such  that  it  is  much  larger  than  the  size  of  the  largest  structure  remaining  correlated. 
In  section  5.3.2  it  was  found  that  the  structure  at  all  wavenumbers  (wavelengths)  decay  over 
distances  between  2  and  2|  times  their  size.  For  the  computational  experiment  investigated, 
the  box  size  was  5,  2  and  5  units  in  the  x,  y,  z  directions,  respectively.  Because  y  =  2 h,  then 
x  =  5  corresponds  to  a  computational  box  only  bh  units  long  in  the  x  direction.  In  order  to 
eliminate  the  periodicity  in  the  boundary  conditions  from  interfering  with  the  results,  this  box 
size  should  be  increased  bv  a  factor  two  to  four  times  its  present  size. 

In  order  to  more  explicitly  compare  the  two  experiments,  a  comparison  was  made  based  on 
outer  and  inner  unit  time  scales.  Figures  5.33  and  5.34  use  a  time  scale  h/Uoc  associated  with 
the  outer  flow,  and  figures  5.35  and  5.36  use  a  time  scale  u/um7  associated  with  flow  close  to  the 
wall.  Figure  5.33  shows  the  largest  three  Reynolds  numbers  for  the  physical  experiment  at  their 
respective  values  of  r\/h  separation.  Two  distinct  features  are  present  as  a  function  of  Reynolds 
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number.  First,  aJ!  three  separations  show  a  frozen  convection  pattern.  In  other  words,  each 
successive  correlation  is  coincident  about  its  value  of  rj/h  for  all  Reynolds  numbers.  Secondly, 
the  correlation  coefficient  has  the  highest  value  and  broadest  curve  for  the  smallest  Reynolds 
number  at  the  first  separation,  but  as  the  flow  convects  and  decays,  all  curves  approach  a  similar 
shape  and  intensity.  Figure  5.34  shows  a  comparison  between  the  computational  experiment 
and  the  lowest  Reynolds  number  at  Rh  =  15,000  for  the  same  separation,  r^/h  =  0.32.  The 
computational  experiment  indicates  a  higher  intensity  and  broader  shape  at  this  separation, 
consistent  with  the  previous  information.  It  may  also  be  noted,  that  the  physical  experiment 
supports  the  decay  properties  shown  in  section  5.3.2  via  a  narrow  band  analysis  with  virtually 
all  turbulence  decaying  within  two  non-dimensional  length  units. 

Figures  5.35  and  5.36  show  similar  trends  when  non-dimensionalized  on  an  inner  time  scale. 
The  difference  is  that  a  convecting  pattern  is  observed  based  on  Reynolds  number.  The  same 
trend  indicated  earlier  is  observed  relative  to  intensity  levels,  but  the  largest  Reynolds  numbers 
include  much  smaller  length  scales  with  resulting  smaller  time  scales  and  show  a  pseudo  moving 
convection  pattern.  Again,  the  similar  trend  is  observed  in  figure  5.36  for  the  computational 
experiment,  with  its  time  scale  being  much  larger  at  its  lower  Reynolds  number  than  for  the 
larger  Reynolds  numbers  of  the  physical  experiment.  The  moving  convecting  pattern  can  be 
explained  by  the  differing  values  of  of  each  experiment.  In  figures  5.33  and  5.34,  each  succes¬ 
sive  correlation  hump  is  at  the  same  rj/h  value  which  does  not  change  with  Reynolds  number. 
For  figures  5.35  and  5.36,  each  experiment  incurs  a  different  value  for  the  shear  velocity  vm , 
hence  =  rju'/x/  has  ever  increasing  values  from  the  computational  experiment  through  the 
physical  experiment.  To  summarize,  the  computational  experiment  represents  similar  results 
to  the  physical  experiment  which  are  consistent  with  its  lower  Reynolds  number.  However,  its 
limited  computational  box  size,  contributed  to  a  much  larger  decay  rate  than  expected  with 
the  periodic  boundary  conditions  contributing  noise  at  the  larger  time  separations,  hence  lower 
frequencies. 

Figure  5.37  presents  the  convection  velocity  based  on  a  broadband  approach.  In  other 
words,  the  time  for  the  maximum  to  occur  at  each  successive  separation  in  the  curves  of  figures 
5.33  thru  5.36  is  combined  with  the  separation  distance  to  obtain  a  convection  velocity  according 
to  equation  (2.43).  Similar  to  the  results  of  Bull  (1967)  and  Blake  (1970),  the  convection 
velocity  increases  with  separation  for  both  experiments.  Both  Bull  and  Blake  showed  that  the 
broadband  convection  velocity  should  begin  about  Ue/Uco  =  0.5  and  increase  with  separation 
to  an  asymptotic  value  of  about  0.75.  The  computational  results  are  somewhat  influenced  by 
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the  periodicity  interfering  with  their  intensity  levels,  and  the  physical  experiment  did  not  have 
enough  separations  to  reach  an  asymptotic  value.  On  the  whole,  the  results  of  both  present 
experiments  are  consistent  with  those  in  the  literature. 

5.4.2  Length  and  Time  Scales 

As  postulated  by  Kraichnan  (1956),  the  integration  of  the  auto-covariance  function  for 
pressure  over  a  flat  plate  to  include  the  wall  of  a  channel,  should  be  equal  to  zero  such  that  no 
net  force  per  unit  area  would  be  present.  However,  Hansen  et  a/.(l 987)  argued  that  channel 
flows  can  support  very  large  scale  low  frequency  structure  such  that  this  integration  is  finite. 
Because  the  spatial  correlation  as  defined  by  equation  (2.47)  was  not  obtained,  a  macro-length 
scale  could  not  be  defined  via  this  analysis.  Therefore,  the  evaluation  of  the  auto  correlation 
functions  at  the  first  ‘zero’  crossing  for  each  Reynolds  number  of  the  physical  experiment  were 
obtained  as  an  approximation  for  macro-time  scales  and  are  indicated  as, 


Rh  =  15,000  — * 

To,  =  6.25  (msec)  , 

Rh  =  20,000  — ► 

To,  =  8.50  (msec)  , 

Rh  ~  25,000  — * 

T0 3  =  6.75  (msec)  . 

From  Taylor’s  hypothesis,  one  can  define  a  macro-length  scale  Lx  from  equation  (2.49)  by 
assuming  that  the  largest  structure  existing  in  the  flow  is  proportional  to  the  maximum  velocity 
and  the  maximum  time  scale.  The  scales  for  the  three  highest  Reynolds  numbers  are, 

Rh  =  15,000  — *  LzJh  =  0.66, 

Rh  =  20,000  — *  LxJh  =  1.10  , 

Rh  =  25,000  — ♦  LxJh  =  1.14  . 

These  scales  give  a  relative  measure  of  the  size  of  the  largest  eddies  which  exist  in  the  flow 
and  contain  most  of  the  energy.  They  are  also  indicative  of  the  longest  distance  over  which 
large  eddies  are  correlated.  From  section  5.3.2  and  5.4.1,  it  was  shown  that  eddies  of  a  specific 
6ize  Ar  decay  in  a  distance  approximately  2  \  times  their  size  in  channel  flows.  The  macro¬ 
scales  indicated  above  suggest  that  it  is  the  limiting  dimension  of  a  channel  (2 h)  which  controls 
the  size  and  extent  over  which  structure  exists  and  remains  coherent.  For  the  computational 
experiment,  the  total  box  size  of  At  =  bh  would  seem  to  be  only  twice  as  long  as  the  correlation 


length  scale  for  the  largest  eddies.  Thus,  only  two  statistically  independent  points  in  the  flow 
direction  would  be  appropriate  for  ensemble  averaging  purposes  with  16  in  the  ‘z’  direction. 
Therefore,  instead  of  an  accurate  statistical  average  of  2048  ensembles  for  the  auto  spectral 
density  function,  one  only  obtains  approximately  16  which,  according  to  equation  (2.52)  gives 
an  error  of, 

«1*pp(/)]  =  0.125.  (5.16) 

Finally,  using  Taylor’s  analysis  presented  in  equations  (2.50)  and  (2.51),  the  time  and  length 
micro  scales  were  computed  and  are  listed  below. 

Rh  =  15,000  — ♦  rmini  =  1.21  (msec)  , 

Rh  =  20,000  — *  rmin2  =  0.93  (msec)  , 

Rh  =  25,000  — «  rmin:s  =  0.56  (msec)  ; 

Rh  =  15,000  — .  lxJh  =  0.128  , 

Rh  =  20,000  —  lxJh  =  0.122, 

Rh  =  25,000  — ♦  lxJh  =  0.094. 

The  definition  of  a  micro-scale  is  an  average  and  certainly  can  be  considered  an  upper  bound 

on  the  size  of  the  smallest  eddies.  The  scale  size  analysis  presented  above,  is  general  and 
requires  further  investigation  of  similar  computations  from  the  velocity  field  which  the  present 
experimental  investigation  was  unable  to  provide. 

5.4.3  Probability  and  Time  Series  Analysis 

Schewe  (1983)  conducted  one  of  the  more  definitive  analyses  on  external  boundary  layer 
flows  for  which  the  transducer  dimension  was  as  small  as  nineteen  viscous  unit6  in  size,  one 
of  the  smallest  amongst  all  investigations.  Schewe  showed  that  the  wall  pressure  intensity 
exceeded  values  as  high  as  pt  —  3 prmi  for  only  one  percent  of  the  time  but  still  contributed 
significantly  to  the  overall  pressure  level  of  intensity.  One  of  the  primary  concerns  of  this 
experimental  investigation  was  that  the  resolution  of  data  obtained  for  the  physical  experiment 
was  sufficient  to  capture  all  relevant  information  for  the  small  scales.  Figures  5.38  and  5.39 
present  the  probability  distribution  functions  (PDF),  as  defined  in  chapter  2,  normalized  by 
Prmi  and  the  dynamic  pressure  Q,  for  both  the  physical  and  computational  experiments.  The 
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data  of  Schewe  is  shown  as  well  as  the  typical  curve  representing  a  Gaussian  random  function. 
The  data  from  both  the  physical  and  computational  experiments  are  in  excellent  agreement 
with  Schewe's  data.  An  analysis  was  conducted  for  both  experiments  to  determine  the  percent 
the  wall  pressure  exceeded  3 prm4.  The  results  of  the  computations  showed  the  same  one 
percent  threshold  reported  by  Schewe  which  is  consistent  with  the  PDF  results  presented  in 
these  figures. 

Representive  time  records  are  plotted  in  figures  5.40  and  5.41  for  the  computational  exper¬ 
iment  and  figures  5.42  and  5.43  for  the  physical  experiment  at  the  Reynolds  numbers  of  15000 
and  25000,  respectively.  Although  a  visual  analysis  of  time  records  is  somewhat  -subjective, 
certain  observations  can  be  noted.  A  low  frequency  periodicity  is  clearly  indicated  in  figures 
5.41  and  5.42  for  the  computational  experiment  which  shows  the  interference  caused  by  the 
periodic  boundary  conditions  used  in  the  computation.  Wall  pressure  fluctuations  typify  a 
very  intermittent  process  in  which  short  time  segments  with  large  fluctuations  appear  to  be 
followed  by  long  time  segments  with  smaller  levels  of  fluctuations.  Events  exceeding  3prm, 
have  short  time  histories  and  decay  much  faster  than  those  of  larger  structures.  Overall,  these 
comments  seem  to  be  substantiated  by  the  the  temporal  results.  If  one  observes  the  temporal 
data  somewhat  subjectively,  a  frozen  pattern  of  turbulence  moving  at  a  convection  velocity  of 
L'c  =  0.70  can  be  observed,  as  shown  by  the  broad  band  analysis  discussed  in  section  5.4.1. 
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Figure  5.1  Power  Spectral  Density  at  R_h  =  10,000 
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Figure  5.2  Power  Spectral  Density  at  R_h  -  15,000 
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Figure  5.4  Power  Spectral  Density  at  R_h  =  25,000 
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Figure  5.5  Highest  Resolution  Power  Spectral  Density  (corrected) 
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Figure  5.6  Magnitude  Squared  Coherence,  R_h  =  10,000. 
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Figure  5.7  Magnitude  Squared  Coherence,  R_h  =  15,000. 


Figure  5.8  Magnitude  Squared  Coherence,  R_h  =  20,000. 
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Figure  5.9  Magnitude  Squared  Coherence,  R_h  =  25,000. 
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Figure  5.10  Comparison  of  RMS  Levels 
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Figure  5.11  Outer  Variable  Similarity 
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Figure  5.12  Mixed  Variable  Similarity 
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Figure  5.13  Mixed  Variable  Similarity  With  Hydraulic  Diameter  on  Both  Axes. 
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Figure  5.14  Mixed  Variable  Similarity  With  Hydraulic  Diameter  on  X-axis  Only. 
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Figure  5.17  Streamwise  Coherence,  R_h  =  15,000 


Figure  5.18  Streamwise  Coherence,  R_h  =  20,000 
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Figure  5.19  Streamwise  Coherence,  R_h  =  25,000 


Figure  5.20  Narrow  Band  Convection  Velocity  of  Computational 
Experiment,  (Outer  Units) 


Figure  5.21  Narrow  Band  Convection  Velocity  of  Physical 
Experiment  At  R_h  =  15,000  (Outer  Units) 
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Figure  5.23  Narrow  Band  Convection  Velocity  of  Physical 
Experiment  At  R_h  =  25,000  (Outer  Units) 
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Figure  5.34  Narrow  Band  Convection  Velocity  of  Computational 
Experiment,  (Inner  Units) 
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Figure  5.26  Narrow  Band  Convection  Velocity  of  Physical 
Experiment  At  R_h  =  20,000  (Inner  Units) 
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Figure  5.27  Narrow  Band  Convection  Velocity  of  Physical 
Experiment  At  R_h  =  25,000  (Inner  Units) 
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Figure  5.28  Time  Correlation  for  Physical  Experiment 
At  R_h  =  10,000. 


Figure  5.29  Time  Correlation  for  Physical  Experiment 
At  R_h  =  15,000. 


Figure  5.31  Time  Correlation  for  Physical  Experiment 
At  R_h  =  25,000. 
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Figure  5.32  Time  Correlation  for  Computational  Experiment 
At  R_h  =  2,129. 


Figure  5.33  Outer  Unit  Correlation  Versus  Reynolds  Number 
for  the  Physical  Experiment. 
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Figure  5  34  Outer  Unit  Correlation;  Computational  and 
Physical  Experiments. 


Figure  5.35  Inner  Unit  Correlation  Versus  Reynolds  Number 
for  the  Physical  Experiment. 
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Figure  5.37  Broadband  Convection  Velocity  for  Both  Experiments. 
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Figure  5.38  Probability  Distribution  Normalized  with  rms  Pressure. 
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Figure  5.39  Probability  Distribution  Normalized  with  Dynamic  Head. 
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5.40  Computational 
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Figure  5.41  Computational  Temporal  Data  at  Station  4. 
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Figure  5.42  Temporal  Data  at  R_h  =  15,000 
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CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 


6.1  Summary  of  Accomplishments 

Throughout  this  experimental  study,  a  number  of  issues  of  scientific  merit  have  been  raised 
and  effort  has  been  focused  on  addressing  these  issues.  Wall  pressure  experiments  in  various 
flow  geometries  have  been  conducted  for  more  than  the  last  three  decades  and  one  might  well 
ask  “What  has  been  learned?”.  Why  does  a  characterization  of  “low  Reynolds  number”  appear 
to  be  so  important?  In  relation  to  pressure  in  an  acoustic  sense,  has  the  knowledge  gained  in 
recent  experiments  increased  our  ability  to  discriminate  against  “flow  noise”?  Questions  also 
have  been  raised  from  a  more  generic  standpoint,  about  “internal  flow”  versus  “external  flow". 
For  instance,  do  unique  length,  velocity  and/or  time  scales  exist  for  both  flow  geometries 
regardless  of  Reynolds  number-and -do  both  flow  regimes  demonstrate  similarity  scaling?  All 
of  the  above  issues  have  played  significant  roles  in  the  design  of  the  present  experiment  to  in 
some  way  further  a  basic  understanding  of  turbulent  flows. 

Three  major  reasons  can  be  identified  from  a  practical  standpoint  to  justify  the  exami¬ 
nation  of  wall  pressure  fluctuations.  First,  a  wall  pressure  fluctuation  experiment  provides  a 
non-intrusive  means  for  examining  the  fundamental  physics  of  turbulent  flows.  By  making  the 
measurement  at  the  wall,  no  flow  disturbance  is  created  by  the  introduction  of  pressure  probes 
within  the  flow  area.  Distinct  regions  of  the  pressure  spectrum  can  be  investigated  for  scaling 
laws  and/or  decay  characteristics,  as  represented  by  similarity  scaling.  Certain  properties  of 
the  space-time  correlation  function  may  give  insights  into  characteristic  length  and  time  scales 
of  the  underlying  turbulence.  Second,  wall  pressure  fluctuation  spectra  provide  information  on 
fluid/structure  interaction.  Vibration  and  sound  are  generated  by  coupling  phenomena  which 
may  lead  to  damaging  fatigue  if  a  resonance  condition  exists.  Third,  “flow  noise”  effects  can 
be  examined  in  relation  to  surface  and  sub-surface  sonar  platforms  and  towed  arrays  from 
the  standpoint  of  contamination  and  detection  liabilities.  These  issues  provided  the  major 
motivation  for  conducting  this  investigation. 

A  recent  study  by  Kim  (1989)  has  raised  interesting  questions  about  the  origins  (source 

♦ 

terms)  of  the  wall  pressure  fluctuations.  Kim  analyzed  the  results  from  a  DNS  calculation 
at  a  larger  Reynolds  number  ( Rm  =  180)  than  that  of  the  present  experiment.  Decomposing 
the  velocity  into  mean  and  fluctuating  components,  he  split  the  pressure  into  a  “rapid”  part 
containing  the  mean  velocity  gradient  and  a  “slow”  part  containing  a  product  of  fluctuating 
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component  gradients.  The  rapid  and  slow  terms  come  about  by  assuming  that  only  the  rapid 
part  responds  immediately  to  a  change  imposed  on  the  mean  field,  whereas  the  slow  part  feels 
the  change  through  non-linear  interactions.  Sometimes,  the  slow  part  is  termed  the  “return" 
term  because  a  number  of  researchers  believe  this  term  to  be  responsible  for  returning  the 
Reynolds-stress  tensor  toward  local  isotropy.  Kim  found  that,  contrary  to  prior  experimental 
work,  the  slower  part  was  much  larger  in  intensity  than  the  rapid  part  except  very  close  to 
the  wall  (Y+  <  15).  Both  parts  were  found  to  contribute  approximately  75  percent  of  the 
total  rms  pressure  at  the  wall.  In  terms  of  single  point  frequency  spectra,  Kim  found  no 
distinct  differences  in  length  scales  associated  with  the  two  pressures.  Because  the  rapid  part 
is  composed  of  the  mean  velocity  gradient,  one  would  expect  large  contributions  from  the  low 
wavenumbers  and  so  this  finding  is  somewhat  surprising.  This  recent  work  of  Kim,  together 
with  the  results  presented  here,  suggest  the  viability  of  using  DNS  calculations  to  further  our 
understanding  and  interpretation  of  the  fundamental  physics  of  turbulent  flows. 

The  overall  objective  of  this  research  was  to  obtain  experimentally  for  the  first  time, 
the  statistical  properties  of  the  wall  pressure  fluctuations  for  a  generic  turbulent  flow  within 
a  rectangular  channel  and  to  directly  compare  these  results  with  those  obtained  from  a  full 
solution  to  the  Navier-Stokes  equations  via  Direct  Numerical  Simulation  (DNS).  For  both  the 
physical  and  computational  experiments,  serious  limitations  on  the  quality  of  the  results  had 
to  be  identified.  For  instance,  facility  generated  noise,  transducer  resolution/sensitivity  and 
Reynolds  number  effects  played  primary  roles  in  the  determination  of  whether  the  physical 
experiment  represented  true  equilibrium  flow  conditions.  Equilibrium  conditions  also  had  to 
be  met  in  the  computational  experiment  but  there  were  additional  constraints  inherent  in 
the  computational  technique  employed.  In  other  words,  computational  accuracy  was  directly 
related  to  the  grid  resolution,  the  size  of  the  computational  box  and  the  extent  in  time  over 
which  the  computation  was  conducted. 

An  extensive  data  base  was  obtained  from  the  present  physical  experiment.  Profiles  for  the 
mean  and  fluctuating  components  of  the  streamwise  velocity  were  obtained  at  three  locations 
in  the  spanwise  direction  for  a  range  of  Reynolds  numbers  J?/,  =  5000,  10000, 15000  and  20000. 
Utilizing  a  new  cancellation  technique,  a  complete  set  of  cross-spectral  density  measurements 
was  obtained  for  the  wall  pressure  field  for  Reynolds  numbers  ranging  from  Rh  =  10,000  to 
25,000  and  at  streamwise  separations  represented  by  £r//i  =0.16  to  0.48.  All  data  was  stored 
on  magnetic  tape,  then  transferred  to  and  subsequently  analyzed  on  a  CRAY  super-computer. 
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The  original  time-dependent  data  together  with  the  processed  results  are  stored  in  CRAY  tape 
files  for  easy  access. 


6.2  Major  Conclusions 

A  number  of  issues  were  investigated  in  the  process  of  conducting  this  research  with  the 
following  conclusions.  First,  a  careful  characterization  of  what  constitutes  a  generic  equilibrium 
experiment  for  wall  pressure  was  obtained.  It  was  found  that  the  physical  experiment  repre¬ 
sented  a  generic  channel  flow  with  fully-developed  profiles  of  streamwise  mean  and  fluctuating 
components  of  velocity.  It  was  further  shown  that  the  small  transducers  utilized  in  obtain¬ 
ing  the  wall  pressure  spectral  properties  provided  sufficient  resolution  to  adequately  resolve 
all  relevant  information  at  the  small  scales  and  at  high  frequencies.  Facility-generated  noise 
was  effectively  eliminated  by  the  application  of  a  new  cancellation  technique  which  thereby 
produced  new  experimental  results  at  low  frequencies  within  internal  flow  geometries. 

A  second  issue  of  importance  was  related  to  the  physics  of  the  flow  through  the  examination 
of  the  fundamental  nature  of  the  wall  pressure  spectrum.  A  broadband  continuous  pressure 
spectrum  with  no  bulges  was  established  for  the  physical  experiment  thereby  demonstrating 
an  equilibrium  channel  flow.  Scaling  laws  were  shown  to  collapse  the  data  over  specific  regions 
related  to  the  predominant  source  terms  at  those  respective  frequencies.  The  single  point 
spectra  of  a  channel  flow  can  be  related  to  similar  spectra  obtained  in  fully  developed  pipe 
flow,  provided  the  length  scale  used  to  collapse  the  data  is  based  on  the  hydraulic  diameter. 
Both  spectra  demonstrate  maximum  spectral  levels  at  fdh/Uoo  =  1.0,  where  /  is  the  frequency 
(Hertz),  dh  represents  the  hydraulic  diameter,  and  Uoc  is  the  centerline  mean  velocity.  An  w,2‘ 
dependence  in  the  pressure  spectrum  amplitude  was  shown  to  exist  at  low  frequencies,  as  has 
been  predicted  by  Patton  et  al. (1980)  and  Chase  (1980).  An  examination  of  the  streamwise  two- 
point  cross-spectral  results  show  that  internal  flows  do  exhibit  similarity  scaling  in  a  manner 
similar  to  external  flows,  wherein  an  eddy  of  a  specific  size  decays  in  a  distance  proportional 
to  its  size.  The  narrow  band  convection  velocity  was  shown  to  approach  zero  at  decreasing 
frequencies,  as  first  shown  by  Blake  (1970)  and  more  recently  substantiated  by  the  extensive 
experimental  work  of  Farabee  (1986). 

These  two  issues  represent  fundamental  concerns  to  experimentalists,  because  low 
Reynolds  number  data  are  not  well  documented  and  are  virtually  non-existent  for  Direct  Nu¬ 
merical  Simulation  channel  flow  experiments.  A  third  issue  representing  the  primary  goal  of  this 
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research  was  to  compare  directly  the  results  from  both  experiments.  Overall,  the  results  from 
the  computational  experiment  showed  good  qualitative  agreement  with  those  obtained  from  the 
physical  experiment.  However,  the  results  from  the  present  research  demonstrated  that  specific 
limitations  can  be  characterized  for  the  results  obtained  from  the  DNS  experiment.  First,  the 
computational  results  are  somewhat  underdeveloped  for  higher  order  statistical  functions.  Two 
possible  explanations  are  that  the  flow  was  not  developed  fully  enough  in  time  before  recording 
the  data,  and  that  statistical  errors  were  introduced  by  calculating  ensemble  averages  via  a 
spatial  average  rather  than  time  averages  based  on  a  stationary  ergodic  process.  Second,  there 
were  resolution  problems  at  both  ends  of  the  frequency  spectrum  related  to  the  computational 
technique  and  the  related  governing  parameters.  The  small  scales  at  high  frequencies  were 
not  adequately  resolved  due  to  the  insufficient  size  of  the  computational  grid.  This  cutoff  at 
high  frequencies  was  denoted  by  the  attenuation  of  spectrum  levels  for  the  fine-scale  structure, 
much  as  has  been  demonstrated  in  physical  experiments  which  use  pressure  transducers  with 
too  large  a  diameter.  However,  physically  large  transducers  do  provide  some  information  (al¬ 
though  attenuated),  whereas  the  computational  technique  acts  as  a  low-pass  filter.  The  large 
grid  spacing  apparently  results  in  the  rejection  of  ‘all’  smaller  scale  information.  Hence,  the 
cut-off  in  frequency  spectra,  as  demonstrated  via  inner  variable  scaling  in  figure  5.15,  is  much 
sharper  than  for  a  physical  transducer  with  the  same  value  of  d+.  The  size  of  the  computational 
box  in  which  the  simulation  was  carried  out  also  affects  the  spectrum  levels  at  low  frequen¬ 
cies.  As  shown  in  chapter  5,  this  box  size  appeared  to  be  too  small.  The  periodic  boundary 
conditions  were  sensed  by  the  flow  and  thus  represented  “noise”  peaks  in  the  wall  pressure 
spectrum  levels  at  low  frequencies.  This  interference  prevented  adequate  resolution  of  the  flow 
at  low  frequencies  and,  furthermore,  it  demonstrated  a  retardation  in  the  decay  characteristics 
as  demonstrated  by  the  results  obtained  from  the  physical  experiment. 

Finally,  an  issue  of  lesser  importance  but  of  potentially  significant  scientific  merit  was  the 
comparison  of  the  wall  pressure  results  from  internal  flows  (as  represented  by  generic  channel 
flows)  with  the  results  of  typical  outer  flow  regimes,  such  as  flat  plate  boundary  layers.  The 
preponderance  of  data  in  the  literature  is  represented  by  the  latter  set  of  experimental  flows, 
at  Reynolds  numbers  more  than  two  orders  of  magnitude  larger  than  those  of  the  present 
investigation.  It  can  be  argued  that  these  data  represent  an  incorrect  data  base,  when  used  to 
critique  the  accuracy  of  the  results  from  DNS  calculations.  The  results  of  the  present  research 
show  two  fundamental  characteristics  which  differentiate  external  flows  from  internal  flows. 
Channel  flows  were  shown  to  decay  at  approximately  twice  the  decay  rate  of  external  flows. 
Inner/mixed /outer  variable  scaling  was  found  to  collapse  the  data  from  both  regimes  separately, 
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but  no  “universal"  law  was  found  to  exist  which  brings  both  internal  and  external  flow  regimes 
into  coincidence.  This  was  most  noticeable  by  observing  the  difference  in  intensity  levels 
between  the  spectra  of  internal  flows  versus  external  flows.  The  pressure  gradient  represents 
an  “edge"  or  boundary  condition  of  a  flat  plate  boundary  layer  flow  and,  thus  is  controlled  by 
the  outer  intermittent  flow.  There  is  no  derivable  relationship  which,  relates  this  outer  flow 
phenomenon  to  shear  stress  at  the  wall,  whereas  within  internal  flows  a  direct  relationship 
does  exist.  Thus,  the  shear  velocity,  which  is  used  to  define  time  scales  of  the  smaller  structure 
(locally  isotropic)  at  high  frequencies,  appears  to  be  coupled  to  the  outer  flow  region  for  internal 
flows.  This  refutes  or  supports  the  hypothesis,  depending  on  one’s  perspective,  that  the  outer 
flow  of  the  two  flow  regimes  does  not  represent  similar  phenomena,  as  has  been  thought  in  the 
past. 


6.3  Original  Contributions  and  Concluding  Remarks 

This  investigation  has  provided  new  information  not  previously  recorded  in  the  literature 
regarding  wall  pressure  physics  related  to  internal  flow  geometries.  In  particular,  an  extensive 
data  base  has  been  obtained  at  low  Reynolds  numbers  and  may  be  used  to  assess  the  quality 
of  typical  Direct  Numerical  Simulation  (DNS)  experiments.  It  must  be  noted  that  the  data 
from  the  computational  experiment  used  in  the  present  study  is  not  the  only  data  base  which 
exists,  and  certainly  not  the  last  which  will  be  obtained.  However,  a  careful  examination  of  the 
results  obtained  from  the  physical  experiment  should  prove  useful  in  determing  those  specific 
governing  parameters  in  future  DNS  experiments  which  will  provide  a  more  accurate  data  base 
for  comparison. 


As  a  final  remark,  the  original  contributions  to  the  literature  as  obtained  in  the  present 
research  are  summarized  by  the  following: 

o  A  new  cancellation  technique  was  used  for  the  first  time  to  effectively  eliminate 
facility-generated  noise  and  thereby  substantiate  an  w2  dependence  at  low 
frequencies,  as  6hown  by  theoretical  predictions. 
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o  The  full  digital  analysis  technique  used  to  process  the  data  from  the  physical 
experiment  demonstrates  an  innovative  methodology,  not  used  previously,  in  the 
analysis  of  wall  pressure  data. 

o  Baseline  wall  pressure  data  at  low  Reynolds  numbers  has  been  presented  for 
the  first  time. 

o  The  comparison  of  the  results  from  a  Direct  Numerical  Simulation  (DNS)  experiment 
with  the  results  obtained  from  a  generic  channel  flow  experiment  has  been  accomplished 
for  the  first  time. 

o  Unique  features  of  internal  flows  were  examined  in  contrast  to  external  flows  and 
distinct  differences  were  shown  to  exist. 
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